SimSST Examples

Mohsen Soltanifar

2023-01-09

SimSST

The goal of SimSST is to simulate stop signal task data based on fixed ssd method and the tracking method.

Installation

You can install the development version of descfarspkg with:

library(dplyr)
#> Warning: package 'dplyr' was built under R version 4.1.3
#> 
#> Attaching package: 'dplyr'
#> The following objects are masked from 'package:stats':
#> 
#>     filter, lag
#> The following objects are masked from 'package:base':
#> 
#>     intersect, setdiff, setequal, union
library(gamlss.dist)
#> Warning: package 'gamlss.dist' was built under R version 4.1.3
#> Loading required package: MASS
#> Warning: package 'MASS' was built under R version 4.1.3
#> 
#> Attaching package: 'MASS'
#> The following object is masked from 'package:dplyr':
#> 
#>     select
library(MASS)
library(SimSST)
set.seed(1)

Example: Fixed SSD method-based Simulation

mySSTdata1 <- 
  simssfixed(
    pid = c("FNLN1","FNLN1"),
    block = c(1,2), 
    n = c(10,10), 
    m = c(4,4), 
    SSD.b = c(220,240),
    dist.go = c("ExG","ExG"),
    theta.go = as.matrix(rbind(c(440,90,90),c(440,90,90))),
    dist.stop = c("ExG","ExG"),
    theta.stop = as.matrix(rbind(c(120,80,70),c(120,80,70)))
  )
mySSTdata1 
#>       Participant.id Block Trial  Inhibition GORT    SSRT    SRRT    SSD   
#>  [1,] "FNLN1"        "1"   "Go"   "-999"     "558.2" "-999"  "-999"  "-999"
#>  [2,] "FNLN1"        "1"   "Go"   "-999"     "426.8" "-999"  "-999"  "-999"
#>  [3,] "FNLN1"        "1"   "Go"   "-999"     "354.3" "-999"  "-999"  "-999"
#>  [4,] "FNLN1"        "1"   "Stop" "1"        "-999"  "104.6" "-999"  "220" 
#>  [5,] "FNLN1"        "1"   "Go"   "-999"     "745.3" "-999"  "-999"  "-999"
#>  [6,] "FNLN1"        "1"   "Go"   "-999"     "568.7" "-999"  "-999"  "-999"
#>  [7,] "FNLN1"        "1"   "Go"   "-999"     "689.8" "-999"  "-999"  "-999"
#>  [8,] "FNLN1"        "1"   "Stop" "1"        "-999"  "94.7"  "-999"  "220" 
#>  [9,] "FNLN1"        "1"   "Stop" "1"        "-999"  "232.3" "-999"  "220" 
#> [10,] "FNLN1"        "1"   "Stop" "1"        "-999"  "153.8" "-999"  "220" 
#> [11,] "FNLN1"        "2"   "Go"   "-999"     "493.2" "-999"  "-999"  "-999"
#> [12,] "FNLN1"        "2"   "Go"   "-999"     "636.1" "-999"  "-999"  "-999"
#> [13,] "FNLN1"        "2"   "Go"   "-999"     "591.7" "-999"  "-999"  "-999"
#> [14,] "FNLN1"        "2"   "Go"   "-999"     "385.1" "-999"  "-999"  "-999"
#> [15,] "FNLN1"        "2"   "Stop" "0"        "-999"  "195.7" "420.9" "240" 
#> [16,] "FNLN1"        "2"   "Stop" "1"        "-999"  "189.8" "-999"  "240" 
#> [17,] "FNLN1"        "2"   "Go"   "-999"     "564.1" "-999"  "-999"  "-999"
#> [18,] "FNLN1"        "2"   "Stop" "1"        "-999"  "267.2" "-999"  "240" 
#> [19,] "FNLN1"        "2"   "Go"   "-999"     "622.3" "-999"  "-999"  "-999"
#> [20,] "FNLN1"        "2"   "Stop" "1"        "-999"  "264.7" "-999"  "240"

Example: Tracking Method-based Simulation

mySSTdata2 <- 
  simsstrack(
    pid = c("FNLN1","FNLN1"),
    block = c(1,2), 
    n = c(10,10), 
    m = c(4,4), 
    SSD.b = c(220,240), 
    dist.go = c("ExG","ExG" ),
    theta.go = as.matrix(rbind(c(440,90,90),c(440,90,90))),
    dist.stop = c("ExG","ExG" ),
    theta.stop = as.matrix(rbind(c(120,80,70),c(120,80,70)))
  )
mySSTdata2
#>       Participant.id Block Trial  Inhibition GORT    SSRT    SRRT    SSD   
#>  [1,] "FNLN1"        "1"   "Stop" "0"        "-999"  "319.4" "446.6" "220" 
#>  [2,] "FNLN1"        "1"   "Go"   "-999"     "670.7" "-999"  "-999"  "-999"
#>  [3,] "FNLN1"        "1"   "Stop" "1"        "-999"  "302.1" "-999"  "170" 
#>  [4,] "FNLN1"        "1"   "Go"   "-999"     "472"   "-999"  "-999"  "-999"
#>  [5,] "FNLN1"        "1"   "Stop" "1"        "-999"  "155.2" "-999"  "220" 
#>  [6,] "FNLN1"        "1"   "Go"   "-999"     "646.7" "-999"  "-999"  "-999"
#>  [7,] "FNLN1"        "1"   "Stop" "0"        "-999"  "262.5" "370.8" "270" 
#>  [8,] "FNLN1"        "1"   "Go"   "-999"     "474.3" "-999"  "-999"  "-999"
#>  [9,] "FNLN1"        "1"   "Go"   "-999"     "470.5" "-999"  "-999"  "-999"
#> [10,] "FNLN1"        "1"   "Go"   "-999"     "512"   "-999"  "-999"  "-999"
#> [11,] "FNLN1"        "2"   "Stop" "0"        "-999"  "270.3" "427.2" "240" 
#> [12,] "FNLN1"        "2"   "Go"   "-999"     "402.8" "-999"  "-999"  "-999"
#> [13,] "FNLN1"        "2"   "Stop" "1"        "-999"  "171.4" "-999"  "190" 
#> [14,] "FNLN1"        "2"   "Go"   "-999"     "440.4" "-999"  "-999"  "-999"
#> [15,] "FNLN1"        "2"   "Stop" "0"        "-999"  "160.2" "392.4" "240" 
#> [16,] "FNLN1"        "2"   "Go"   "-999"     "351.9" "-999"  "-999"  "-999"
#> [17,] "FNLN1"        "2"   "Stop" "0"        "-999"  "276"   "442.4" "190" 
#> [18,] "FNLN1"        "2"   "Go"   "-999"     "562.5" "-999"  "-999"  "-999"
#> [19,] "FNLN1"        "2"   "Go"   "-999"     "671.4" "-999"  "-999"  "-999"
#> [20,] "FNLN1"        "2"   "Go"   "-999"     "615"   "-999"  "-999"  "-999"

To convert tracking-based SimSST simulated data to BEESTS software input data, we have:

Example: Converstion to BEESTs software format

Datatemp2 <- mySSTdata2 
ss_presented <- recode(Datatemp2[,3], 'Stop' = "1", 'Go' = "0")
inhibited <- Datatemp2[,4]
ssd <- Datatemp2[,8]
rt <- Datatemp2[,5]
srrt <- Datatemp2[,7]
Data2 <- cbind.data.frame(ss_presented, inhibited, ssd, rt, srrt)
for(i in 1:20) if(Data2$inhibited[i]==0) Data2$rt[i] <- Data2$srrt[i]  
myBEESTSdata2 <- (Data2[,-5])[order(ss_presented),]
myBEESTSdata2
#>    ss_presented inhibited  ssd    rt
#> 2             0      -999 -999 670.7
#> 4             0      -999 -999   472
#> 6             0      -999 -999 646.7
#> 8             0      -999 -999 474.3
#> 9             0      -999 -999 470.5
#> 10            0      -999 -999   512
#> 12            0      -999 -999 402.8
#> 14            0      -999 -999 440.4
#> 16            0      -999 -999 351.9
#> 18            0      -999 -999 562.5
#> 19            0      -999 -999 671.4
#> 20            0      -999 -999   615
#> 1             1         0  220 446.6
#> 3             1         1  170  -999
#> 5             1         1  220  -999
#> 7             1         0  270 370.8
#> 11            1         0  240 427.2
#> 13            1         1  190  -999
#> 15            1         0  240 392.4
#> 17            1         0  190 442.4

Once users reformat the simulated datasets to BEESTS dataset format, they may feed the BEESTS data to the BEESTS software to compute their plausible parameters (mu, sigma, tau).

Example: Simulating correlated SST data using general tracking method

mySSTdata3 <- simssgen(
     pid = c("FNLN1","FNLN2","FNLN2"),
     block = c(1,1,2),
     n = c(50,100,150),
     m = c(10,20,30),
     SSD.b = c(200,220,240),
     dist.go = c("ExG","ExG","ExG"),
     theta.go = as.matrix(rbind(c(400,60,30),c(440,90,90),c(440,90,90))),
     dist.stop = c("ExG","ExG","ExG"),
     theta.stop = as.matrix(rbind(c(100,70,60),c(120,80,70),c(120,80,70))),
     rho = c(0.35,0.45,0.45),
     d = c(50,65,75))
mySSTdata3
#>        Participant.id Block Trial  Inhibition GORT    SSRT    SRRT    SSD   
#>   [1,] "FNLN1"        "1"   "Stop" "0"        "-999"  "411.7" "384"   "200" 
#>   [2,] "FNLN1"        "1"   "Go"   "-999"     "421.1" "-999"  "-999"  "-999"
#>   [3,] "FNLN1"        "1"   "Stop" "1"        "-999"  "265.1" "-999"  "150" 
#>   [4,] "FNLN1"        "1"   "Go"   "-999"     "446.5" "-999"  "-999"  "-999"
#>   [5,] "FNLN1"        "1"   "Stop" "1"        "-999"  "214.1" "-999"  "200" 
#>   [6,] "FNLN1"        "1"   "Go"   "-999"     "402.6" "-999"  "-999"  "-999"
#>   [7,] "FNLN1"        "1"   "Stop" "1"        "-999"  "77.8"  "-999"  "250" 
#>   [8,] "FNLN1"        "1"   "Go"   "-999"     "392.8" "-999"  "-999"  "-999"
#>   [9,] "FNLN1"        "1"   "Stop" "0"        "-999"  "240.9" "359.9" "300" 
#>  [10,] "FNLN1"        "1"   "Go"   "-999"     "435.8" "-999"  "-999"  "-999"
#>  [11,] "FNLN1"        "1"   "Stop" "0"        "-999"  "149.8" "328.5" "250" 
#>  [12,] "FNLN1"        "1"   "Go"   "-999"     "423.5" "-999"  "-999"  "-999"
#>  [13,] "FNLN1"        "1"   "Stop" "0"        "-999"  "104.2" "301.1" "200" 
#>  [14,] "FNLN1"        "1"   "Go"   "-999"     "426.8" "-999"  "-999"  "-999"
#>  [15,] "FNLN1"        "1"   "Stop" "1"        "-999"  "225.2" "-999"  "150" 
#>  [16,] "FNLN1"        "1"   "Go"   "-999"     "433.8" "-999"  "-999"  "-999"
#>  [17,] "FNLN1"        "1"   "Stop" "1"        "-999"  "121.4" "-999"  "200" 
#>  [18,] "FNLN1"        "1"   "Go"   "-999"     "435"   "-999"  "-999"  "-999"
#>  [19,] "FNLN1"        "1"   "Stop" "0"        "-999"  "313.9" "377"   "250" 
#>  [20,] "FNLN1"        "1"   "Go"   "-999"     "432.7" "-999"  "-999"  "-999"
#>  [21,] "FNLN1"        "1"   "Go"   "-999"     "451.9" "-999"  "-999"  "-999"
#>  [22,] "FNLN1"        "1"   "Go"   "-999"     "508.4" "-999"  "-999"  "-999"
#>  [23,] "FNLN1"        "1"   "Go"   "-999"     "429.6" "-999"  "-999"  "-999"
#>  [24,] "FNLN1"        "1"   "Go"   "-999"     "381.2" "-999"  "-999"  "-999"
#>  [25,] "FNLN1"        "1"   "Go"   "-999"     "412.3" "-999"  "-999"  "-999"
#>  [26,] "FNLN1"        "1"   "Go"   "-999"     "424.5" "-999"  "-999"  "-999"
#>  [27,] "FNLN1"        "1"   "Go"   "-999"     "351.9" "-999"  "-999"  "-999"
#>  [28,] "FNLN1"        "1"   "Go"   "-999"     "331.3" "-999"  "-999"  "-999"
#>  [29,] "FNLN1"        "1"   "Go"   "-999"     "517.5" "-999"  "-999"  "-999"
#>  [30,] "FNLN1"        "1"   "Go"   "-999"     "440.8" "-999"  "-999"  "-999"
#>  [31,] "FNLN1"        "1"   "Go"   "-999"     "452.8" "-999"  "-999"  "-999"
#>  [32,] "FNLN1"        "1"   "Go"   "-999"     "380.3" "-999"  "-999"  "-999"
#>  [33,] "FNLN1"        "1"   "Go"   "-999"     "441.1" "-999"  "-999"  "-999"
#>  [34,] "FNLN1"        "1"   "Go"   "-999"     "533.1" "-999"  "-999"  "-999"
#>  [35,] "FNLN1"        "1"   "Go"   "-999"     "449"   "-999"  "-999"  "-999"
#>  [36,] "FNLN1"        "1"   "Go"   "-999"     "533.8" "-999"  "-999"  "-999"
#>  [37,] "FNLN1"        "1"   "Go"   "-999"     "499.3" "-999"  "-999"  "-999"
#>  [38,] "FNLN1"        "1"   "Go"   "-999"     "472.2" "-999"  "-999"  "-999"
#>  [39,] "FNLN1"        "1"   "Go"   "-999"     "379"   "-999"  "-999"  "-999"
#>  [40,] "FNLN1"        "1"   "Go"   "-999"     "440.9" "-999"  "-999"  "-999"
#>  [41,] "FNLN1"        "1"   "Go"   "-999"     "420.8" "-999"  "-999"  "-999"
#>  [42,] "FNLN1"        "1"   "Go"   "-999"     "500.5" "-999"  "-999"  "-999"
#>  [43,] "FNLN1"        "1"   "Go"   "-999"     "463.8" "-999"  "-999"  "-999"
#>  [44,] "FNLN1"        "1"   "Go"   "-999"     "454.8" "-999"  "-999"  "-999"
#>  [45,] "FNLN1"        "1"   "Go"   "-999"     "528.5" "-999"  "-999"  "-999"
#>  [46,] "FNLN1"        "1"   "Go"   "-999"     "423.3" "-999"  "-999"  "-999"
#>  [47,] "FNLN1"        "1"   "Go"   "-999"     "531"   "-999"  "-999"  "-999"
#>  [48,] "FNLN1"        "1"   "Go"   "-999"     "416.4" "-999"  "-999"  "-999"
#>  [49,] "FNLN1"        "1"   "Go"   "-999"     "434.8" "-999"  "-999"  "-999"
#>  [50,] "FNLN1"        "1"   "Go"   "-999"     "398.7" "-999"  "-999"  "-999"
#>  [51,] "FNLN2"        "1"   "Stop" "1"        "-999"  "76.2"  "-999"  "220" 
#>  [52,] "FNLN2"        "1"   "Go"   "-999"     "423.1" "-999"  "-999"  "-999"
#>  [53,] "FNLN2"        "1"   "Stop" "1"        "-999"  "88.6"  "-999"  "285" 
#>  [54,] "FNLN2"        "1"   "Go"   "-999"     "434.1" "-999"  "-999"  "-999"
#>  [55,] "FNLN2"        "1"   "Stop" "1"        "-999"  "32.5"  "-999"  "350" 
#>  [56,] "FNLN2"        "1"   "Go"   "-999"     "554.1" "-999"  "-999"  "-999"
#>  [57,] "FNLN2"        "1"   "Stop" "0"        "-999"  "283.4" "515.1" "415" 
#>  [58,] "FNLN2"        "1"   "Go"   "-999"     "461.4" "-999"  "-999"  "-999"
#>  [59,] "FNLN2"        "1"   "Stop" "1"        "-999"  "85"    "-999"  "350" 
#>  [60,] "FNLN2"        "1"   "Go"   "-999"     "767.6" "-999"  "-999"  "-999"
#>  [61,] "FNLN2"        "1"   "Stop" "0"        "-999"  "171.1" "444.3" "415" 
#>  [62,] "FNLN2"        "1"   "Go"   "-999"     "541.6" "-999"  "-999"  "-999"
#>  [63,] "FNLN2"        "1"   "Stop" "0"        "-999"  "269.4" "596.8" "350" 
#>  [64,] "FNLN2"        "1"   "Go"   "-999"     "407.3" "-999"  "-999"  "-999"
#>  [65,] "FNLN2"        "1"   "Stop" "0"        "-999"  "343.1" "510.2" "285" 
#>  [66,] "FNLN2"        "1"   "Go"   "-999"     "630.8" "-999"  "-999"  "-999"
#>  [67,] "FNLN2"        "1"   "Stop" "0"        "-999"  "439.3" "558.9" "220" 
#>  [68,] "FNLN2"        "1"   "Go"   "-999"     "469.7" "-999"  "-999"  "-999"
#>  [69,] "FNLN2"        "1"   "Stop" "1"        "-999"  "87.4"  "-999"  "155" 
#>  [70,] "FNLN2"        "1"   "Go"   "-999"     "640.8" "-999"  "-999"  "-999"
#>  [71,] "FNLN2"        "1"   "Stop" "1"        "-999"  "348.7" "-999"  "220" 
#>  [72,] "FNLN2"        "1"   "Go"   "-999"     "639.5" "-999"  "-999"  "-999"
#>  [73,] "FNLN2"        "1"   "Stop" "1"        "-999"  "129.9" "-999"  "285" 
#>  [74,] "FNLN2"        "1"   "Go"   "-999"     "519.3" "-999"  "-999"  "-999"
#>  [75,] "FNLN2"        "1"   "Stop" "1"        "-999"  "145.4" "-999"  "350" 
#>  [76,] "FNLN2"        "1"   "Go"   "-999"     "677.1" "-999"  "-999"  "-999"
#>  [77,] "FNLN2"        "1"   "Stop" "0"        "-999"  "245.1" "340.2" "415" 
#>  [78,] "FNLN2"        "1"   "Go"   "-999"     "446.6" "-999"  "-999"  "-999"
#>  [79,] "FNLN2"        "1"   "Stop" "1"        "-999"  "121"   "-999"  "350" 
#>  [80,] "FNLN2"        "1"   "Go"   "-999"     "259.5" "-999"  "-999"  "-999"
#>  [81,] "FNLN2"        "1"   "Stop" "0"        "-999"  "192.7" "389.5" "415" 
#>  [82,] "FNLN2"        "1"   "Go"   "-999"     "536.1" "-999"  "-999"  "-999"
#>  [83,] "FNLN2"        "1"   "Stop" "0"        "-999"  "195.6" "425"   "350" 
#>  [84,] "FNLN2"        "1"   "Go"   "-999"     "457.1" "-999"  "-999"  "-999"
#>  [85,] "FNLN2"        "1"   "Stop" "1"        "-999"  "132.2" "-999"  "285" 
#>  [86,] "FNLN2"        "1"   "Go"   "-999"     "394.8" "-999"  "-999"  "-999"
#>  [87,] "FNLN2"        "1"   "Stop" "1"        "-999"  "280.5" "-999"  "350" 
#>  [88,] "FNLN2"        "1"   "Go"   "-999"     "664.4" "-999"  "-999"  "-999"
#>  [89,] "FNLN2"        "1"   "Stop" "0"        "-999"  "115.9" "511.3" "415" 
#>  [90,] "FNLN2"        "1"   "Go"   "-999"     "357.1" "-999"  "-999"  "-999"
#>  [91,] "FNLN2"        "1"   "Go"   "-999"     "556.9" "-999"  "-999"  "-999"
#>  [92,] "FNLN2"        "1"   "Go"   "-999"     "559.8" "-999"  "-999"  "-999"
#>  [93,] "FNLN2"        "1"   "Go"   "-999"     "701.3" "-999"  "-999"  "-999"
#>  [94,] "FNLN2"        "1"   "Go"   "-999"     "544.2" "-999"  "-999"  "-999"
#>  [95,] "FNLN2"        "1"   "Go"   "-999"     "585.1" "-999"  "-999"  "-999"
#>  [96,] "FNLN2"        "1"   "Go"   "-999"     "537.8" "-999"  "-999"  "-999"
#>  [97,] "FNLN2"        "1"   "Go"   "-999"     "515.2" "-999"  "-999"  "-999"
#>  [98,] "FNLN2"        "1"   "Go"   "-999"     "529.9" "-999"  "-999"  "-999"
#>  [99,] "FNLN2"        "1"   "Go"   "-999"     "868.7" "-999"  "-999"  "-999"
#> [100,] "FNLN2"        "1"   "Go"   "-999"     "546.9" "-999"  "-999"  "-999"
#> [101,] "FNLN2"        "1"   "Go"   "-999"     "717.4" "-999"  "-999"  "-999"
#> [102,] "FNLN2"        "1"   "Go"   "-999"     "402.9" "-999"  "-999"  "-999"
#> [103,] "FNLN2"        "1"   "Go"   "-999"     "506.4" "-999"  "-999"  "-999"
#> [104,] "FNLN2"        "1"   "Go"   "-999"     "497"   "-999"  "-999"  "-999"
#> [105,] "FNLN2"        "1"   "Go"   "-999"     "428.5" "-999"  "-999"  "-999"
#> [106,] "FNLN2"        "1"   "Go"   "-999"     "488.7" "-999"  "-999"  "-999"
#> [107,] "FNLN2"        "1"   "Go"   "-999"     "402.6" "-999"  "-999"  "-999"
#> [108,] "FNLN2"        "1"   "Go"   "-999"     "343.6" "-999"  "-999"  "-999"
#> [109,] "FNLN2"        "1"   "Go"   "-999"     "435.5" "-999"  "-999"  "-999"
#> [110,] "FNLN2"        "1"   "Go"   "-999"     "367.4" "-999"  "-999"  "-999"
#> [111,] "FNLN2"        "1"   "Go"   "-999"     "677.3" "-999"  "-999"  "-999"
#> [112,] "FNLN2"        "1"   "Go"   "-999"     "575.8" "-999"  "-999"  "-999"
#> [113,] "FNLN2"        "1"   "Go"   "-999"     "648.6" "-999"  "-999"  "-999"
#> [114,] "FNLN2"        "1"   "Go"   "-999"     "412.3" "-999"  "-999"  "-999"
#> [115,] "FNLN2"        "1"   "Go"   "-999"     "568.9" "-999"  "-999"  "-999"
#> [116,] "FNLN2"        "1"   "Go"   "-999"     "461.6" "-999"  "-999"  "-999"
#> [117,] "FNLN2"        "1"   "Go"   "-999"     "319.7" "-999"  "-999"  "-999"
#> [118,] "FNLN2"        "1"   "Go"   "-999"     "562.5" "-999"  "-999"  "-999"
#> [119,] "FNLN2"        "1"   "Go"   "-999"     "533.6" "-999"  "-999"  "-999"
#> [120,] "FNLN2"        "1"   "Go"   "-999"     "441"   "-999"  "-999"  "-999"
#> [121,] "FNLN2"        "1"   "Go"   "-999"     "455.9" "-999"  "-999"  "-999"
#> [122,] "FNLN2"        "1"   "Go"   "-999"     "611.4" "-999"  "-999"  "-999"
#> [123,] "FNLN2"        "1"   "Go"   "-999"     "735.9" "-999"  "-999"  "-999"
#> [124,] "FNLN2"        "1"   "Go"   "-999"     "582.8" "-999"  "-999"  "-999"
#> [125,] "FNLN2"        "1"   "Go"   "-999"     "378.9" "-999"  "-999"  "-999"
#> [126,] "FNLN2"        "1"   "Go"   "-999"     "496.6" "-999"  "-999"  "-999"
#> [127,] "FNLN2"        "1"   "Go"   "-999"     "506.1" "-999"  "-999"  "-999"
#> [128,] "FNLN2"        "1"   "Go"   "-999"     "698.7" "-999"  "-999"  "-999"
#> [129,] "FNLN2"        "1"   "Go"   "-999"     "654"   "-999"  "-999"  "-999"
#> [130,] "FNLN2"        "1"   "Go"   "-999"     "679.3" "-999"  "-999"  "-999"
#> [131,] "FNLN2"        "1"   "Go"   "-999"     "557.5" "-999"  "-999"  "-999"
#> [132,] "FNLN2"        "1"   "Go"   "-999"     "420.6" "-999"  "-999"  "-999"
#> [133,] "FNLN2"        "1"   "Go"   "-999"     "665.3" "-999"  "-999"  "-999"
#> [134,] "FNLN2"        "1"   "Go"   "-999"     "414.5" "-999"  "-999"  "-999"
#> [135,] "FNLN2"        "1"   "Go"   "-999"     "573.2" "-999"  "-999"  "-999"
#> [136,] "FNLN2"        "1"   "Go"   "-999"     "648.3" "-999"  "-999"  "-999"
#> [137,] "FNLN2"        "1"   "Go"   "-999"     "416.6" "-999"  "-999"  "-999"
#> [138,] "FNLN2"        "1"   "Go"   "-999"     "336.1" "-999"  "-999"  "-999"
#> [139,] "FNLN2"        "1"   "Go"   "-999"     "499.6" "-999"  "-999"  "-999"
#> [140,] "FNLN2"        "1"   "Go"   "-999"     "449.1" "-999"  "-999"  "-999"
#> [141,] "FNLN2"        "1"   "Go"   "-999"     "474.3" "-999"  "-999"  "-999"
#> [142,] "FNLN2"        "1"   "Go"   "-999"     "473.9" "-999"  "-999"  "-999"
#> [143,] "FNLN2"        "1"   "Go"   "-999"     "715.3" "-999"  "-999"  "-999"
#> [144,] "FNLN2"        "1"   "Go"   "-999"     "618.6" "-999"  "-999"  "-999"
#> [145,] "FNLN2"        "1"   "Go"   "-999"     "704.4" "-999"  "-999"  "-999"
#> [146,] "FNLN2"        "1"   "Go"   "-999"     "610.6" "-999"  "-999"  "-999"
#> [147,] "FNLN2"        "1"   "Go"   "-999"     "530"   "-999"  "-999"  "-999"
#> [148,] "FNLN2"        "1"   "Go"   "-999"     "573.6" "-999"  "-999"  "-999"
#> [149,] "FNLN2"        "1"   "Go"   "-999"     "609.2" "-999"  "-999"  "-999"
#> [150,] "FNLN2"        "1"   "Go"   "-999"     "626.2" "-999"  "-999"  "-999"
#> [151,] "FNLN2"        "2"   "Stop" "1"        "-999"  "133.1" "-999"  "240" 
#> [152,] "FNLN2"        "2"   "Go"   "-999"     "380.4" "-999"  "-999"  "-999"
#> [153,] "FNLN2"        "2"   "Stop" "0"        "-999"  "265.4" "553.9" "315" 
#> [154,] "FNLN2"        "2"   "Go"   "-999"     "476.7" "-999"  "-999"  "-999"
#> [155,] "FNLN2"        "2"   "Stop" "1"        "-999"  "87.9"  "-999"  "240" 
#> [156,] "FNLN2"        "2"   "Go"   "-999"     "577.3" "-999"  "-999"  "-999"
#> [157,] "FNLN2"        "2"   "Stop" "0"        "-999"  "207.9" "473.3" "315" 
#> [158,] "FNLN2"        "2"   "Go"   "-999"     "346.6" "-999"  "-999"  "-999"
#> [159,] "FNLN2"        "2"   "Stop" "1"        "-999"  "99.6"  "-999"  "240" 
#> [160,] "FNLN2"        "2"   "Go"   "-999"     "558.8" "-999"  "-999"  "-999"
#> [161,] "FNLN2"        "2"   "Stop" "1"        "-999"  "-33.6" "-999"  "315" 
#> [162,] "FNLN2"        "2"   "Go"   "-999"     "566"   "-999"  "-999"  "-999"
#> [163,] "FNLN2"        "2"   "Stop" "0"        "-999"  "288.7" "451.7" "390" 
#> [164,] "FNLN2"        "2"   "Go"   "-999"     "544"   "-999"  "-999"  "-999"
#> [165,] "FNLN2"        "2"   "Stop" "1"        "-999"  "49.6"  "-999"  "315" 
#> [166,] "FNLN2"        "2"   "Go"   "-999"     "619.1" "-999"  "-999"  "-999"
#> [167,] "FNLN2"        "2"   "Stop" "0"        "-999"  "486.7" "394.4" "390" 
#> [168,] "FNLN2"        "2"   "Go"   "-999"     "533"   "-999"  "-999"  "-999"
#> [169,] "FNLN2"        "2"   "Stop" "1"        "-999"  "122.4" "-999"  "315" 
#> [170,] "FNLN2"        "2"   "Go"   "-999"     "404.9" "-999"  "-999"  "-999"
#> [171,] "FNLN2"        "2"   "Stop" "1"        "-999"  "90.4"  "-999"  "390" 
#> [172,] "FNLN2"        "2"   "Go"   "-999"     "414.1" "-999"  "-999"  "-999"
#> [173,] "FNLN2"        "2"   "Stop" "0"        "-999"  "210"   "306.4" "465" 
#> [174,] "FNLN2"        "2"   "Go"   "-999"     "543.4" "-999"  "-999"  "-999"
#> [175,] "FNLN2"        "2"   "Stop" "0"        "-999"  "295.3" "574.3" "390" 
#> [176,] "FNLN2"        "2"   "Go"   "-999"     "416.4" "-999"  "-999"  "-999"
#> [177,] "FNLN2"        "2"   "Stop" "1"        "-999"  "116.8" "-999"  "315" 
#> [178,] "FNLN2"        "2"   "Go"   "-999"     "442.5" "-999"  "-999"  "-999"
#> [179,] "FNLN2"        "2"   "Stop" "1"        "-999"  "139.2" "-999"  "390" 
#> [180,] "FNLN2"        "2"   "Go"   "-999"     "791.7" "-999"  "-999"  "-999"
#> [181,] "FNLN2"        "2"   "Stop" "0"        "-999"  "397"   "543.7" "465" 
#> [182,] "FNLN2"        "2"   "Go"   "-999"     "506.2" "-999"  "-999"  "-999"
#> [183,] "FNLN2"        "2"   "Stop" "1"        "-999"  "34.7"  "-999"  "390" 
#> [184,] "FNLN2"        "2"   "Go"   "-999"     "593.6" "-999"  "-999"  "-999"
#> [185,] "FNLN2"        "2"   "Stop" "0"        "-999"  "231.3" "571.5" "465" 
#> [186,] "FNLN2"        "2"   "Go"   "-999"     "695.9" "-999"  "-999"  "-999"
#> [187,] "FNLN2"        "2"   "Stop" "0"        "-999"  "163.2" "392.4" "390" 
#> [188,] "FNLN2"        "2"   "Go"   "-999"     "526.6" "-999"  "-999"  "-999"
#> [189,] "FNLN2"        "2"   "Stop" "1"        "-999"  "51.5"  "-999"  "315" 
#> [190,] "FNLN2"        "2"   "Go"   "-999"     "558.2" "-999"  "-999"  "-999"
#> [191,] "FNLN2"        "2"   "Stop" "1"        "-999"  "4"     "-999"  "390" 
#> [192,] "FNLN2"        "2"   "Go"   "-999"     "684.1" "-999"  "-999"  "-999"
#> [193,] "FNLN2"        "2"   "Stop" "0"        "-999"  "276.9" "487.3" "465" 
#> [194,] "FNLN2"        "2"   "Go"   "-999"     "350"   "-999"  "-999"  "-999"
#> [195,] "FNLN2"        "2"   "Stop" "0"        "-999"  "156"   "478"   "390" 
#> [196,] "FNLN2"        "2"   "Go"   "-999"     "458.2" "-999"  "-999"  "-999"
#> [197,] "FNLN2"        "2"   "Stop" "1"        "-999"  "131.8" "-999"  "315" 
#> [198,] "FNLN2"        "2"   "Go"   "-999"     "544.7" "-999"  "-999"  "-999"
#> [199,] "FNLN2"        "2"   "Stop" "0"        "-999"  "141.1" "518.2" "390" 
#> [200,] "FNLN2"        "2"   "Go"   "-999"     "576.5" "-999"  "-999"  "-999"
#> [201,] "FNLN2"        "2"   "Stop" "1"        "-999"  "48.4"  "-999"  "315" 
#> [202,] "FNLN2"        "2"   "Go"   "-999"     "447.2" "-999"  "-999"  "-999"
#> [203,] "FNLN2"        "2"   "Stop" "0"        "-999"  "99.4"  "383.4" "390" 
#> [204,] "FNLN2"        "2"   "Go"   "-999"     "574.7" "-999"  "-999"  "-999"
#> [205,] "FNLN2"        "2"   "Stop" "1"        "-999"  "131.7" "-999"  "315" 
#> [206,] "FNLN2"        "2"   "Go"   "-999"     "483.9" "-999"  "-999"  "-999"
#> [207,] "FNLN2"        "2"   "Stop" "1"        "-999"  "107.5" "-999"  "390" 
#> [208,] "FNLN2"        "2"   "Go"   "-999"     "598.1" "-999"  "-999"  "-999"
#> [209,] "FNLN2"        "2"   "Stop" "0"        "-999"  "209.7" "363.8" "465" 
#> [210,] "FNLN2"        "2"   "Go"   "-999"     "298.9" "-999"  "-999"  "-999"
#> [211,] "FNLN2"        "2"   "Go"   "-999"     "328.9" "-999"  "-999"  "-999"
#> [212,] "FNLN2"        "2"   "Go"   "-999"     "416.3" "-999"  "-999"  "-999"
#> [213,] "FNLN2"        "2"   "Go"   "-999"     "553.7" "-999"  "-999"  "-999"
#> [214,] "FNLN2"        "2"   "Go"   "-999"     "755.1" "-999"  "-999"  "-999"
#> [215,] "FNLN2"        "2"   "Go"   "-999"     "461"   "-999"  "-999"  "-999"
#> [216,] "FNLN2"        "2"   "Go"   "-999"     "472.4" "-999"  "-999"  "-999"
#> [217,] "FNLN2"        "2"   "Go"   "-999"     "691.3" "-999"  "-999"  "-999"
#> [218,] "FNLN2"        "2"   "Go"   "-999"     "574"   "-999"  "-999"  "-999"
#> [219,] "FNLN2"        "2"   "Go"   "-999"     "381.6" "-999"  "-999"  "-999"
#> [220,] "FNLN2"        "2"   "Go"   "-999"     "581.1" "-999"  "-999"  "-999"
#> [221,] "FNLN2"        "2"   "Go"   "-999"     "538.9" "-999"  "-999"  "-999"
#> [222,] "FNLN2"        "2"   "Go"   "-999"     "405.3" "-999"  "-999"  "-999"
#> [223,] "FNLN2"        "2"   "Go"   "-999"     "594.2" "-999"  "-999"  "-999"
#> [224,] "FNLN2"        "2"   "Go"   "-999"     "436.3" "-999"  "-999"  "-999"
#> [225,] "FNLN2"        "2"   "Go"   "-999"     "618.3" "-999"  "-999"  "-999"
#> [226,] "FNLN2"        "2"   "Go"   "-999"     "507.2" "-999"  "-999"  "-999"
#> [227,] "FNLN2"        "2"   "Go"   "-999"     "344.1" "-999"  "-999"  "-999"
#> [228,] "FNLN2"        "2"   "Go"   "-999"     "535.6" "-999"  "-999"  "-999"
#> [229,] "FNLN2"        "2"   "Go"   "-999"     "506.7" "-999"  "-999"  "-999"
#> [230,] "FNLN2"        "2"   "Go"   "-999"     "475.2" "-999"  "-999"  "-999"
#> [231,] "FNLN2"        "2"   "Go"   "-999"     "456.1" "-999"  "-999"  "-999"
#> [232,] "FNLN2"        "2"   "Go"   "-999"     "598.3" "-999"  "-999"  "-999"
#> [233,] "FNLN2"        "2"   "Go"   "-999"     "431.2" "-999"  "-999"  "-999"
#> [234,] "FNLN2"        "2"   "Go"   "-999"     "644.7" "-999"  "-999"  "-999"
#> [235,] "FNLN2"        "2"   "Go"   "-999"     "461.5" "-999"  "-999"  "-999"
#> [236,] "FNLN2"        "2"   "Go"   "-999"     "591.9" "-999"  "-999"  "-999"
#> [237,] "FNLN2"        "2"   "Go"   "-999"     "587.3" "-999"  "-999"  "-999"
#> [238,] "FNLN2"        "2"   "Go"   "-999"     "633"   "-999"  "-999"  "-999"
#> [239,] "FNLN2"        "2"   "Go"   "-999"     "395.3" "-999"  "-999"  "-999"
#> [240,] "FNLN2"        "2"   "Go"   "-999"     "516.8" "-999"  "-999"  "-999"
#> [241,] "FNLN2"        "2"   "Go"   "-999"     "580.7" "-999"  "-999"  "-999"
#> [242,] "FNLN2"        "2"   "Go"   "-999"     "638.7" "-999"  "-999"  "-999"
#> [243,] "FNLN2"        "2"   "Go"   "-999"     "540.9" "-999"  "-999"  "-999"
#> [244,] "FNLN2"        "2"   "Go"   "-999"     "406.2" "-999"  "-999"  "-999"
#> [245,] "FNLN2"        "2"   "Go"   "-999"     "535.9" "-999"  "-999"  "-999"
#> [246,] "FNLN2"        "2"   "Go"   "-999"     "566.4" "-999"  "-999"  "-999"
#> [247,] "FNLN2"        "2"   "Go"   "-999"     "726"   "-999"  "-999"  "-999"
#> [248,] "FNLN2"        "2"   "Go"   "-999"     "546"   "-999"  "-999"  "-999"
#> [249,] "FNLN2"        "2"   "Go"   "-999"     "674.1" "-999"  "-999"  "-999"
#> [250,] "FNLN2"        "2"   "Go"   "-999"     "522.2" "-999"  "-999"  "-999"
#> [251,] "FNLN2"        "2"   "Go"   "-999"     "492.1" "-999"  "-999"  "-999"
#> [252,] "FNLN2"        "2"   "Go"   "-999"     "463.4" "-999"  "-999"  "-999"
#> [253,] "FNLN2"        "2"   "Go"   "-999"     "569.2" "-999"  "-999"  "-999"
#> [254,] "FNLN2"        "2"   "Go"   "-999"     "607.9" "-999"  "-999"  "-999"
#> [255,] "FNLN2"        "2"   "Go"   "-999"     "766.2" "-999"  "-999"  "-999"
#> [256,] "FNLN2"        "2"   "Go"   "-999"     "552.5" "-999"  "-999"  "-999"
#> [257,] "FNLN2"        "2"   "Go"   "-999"     "354"   "-999"  "-999"  "-999"
#> [258,] "FNLN2"        "2"   "Go"   "-999"     "481.3" "-999"  "-999"  "-999"
#> [259,] "FNLN2"        "2"   "Go"   "-999"     "573.1" "-999"  "-999"  "-999"
#> [260,] "FNLN2"        "2"   "Go"   "-999"     "739.4" "-999"  "-999"  "-999"
#> [261,] "FNLN2"        "2"   "Go"   "-999"     "610.3" "-999"  "-999"  "-999"
#> [262,] "FNLN2"        "2"   "Go"   "-999"     "644.2" "-999"  "-999"  "-999"
#> [263,] "FNLN2"        "2"   "Go"   "-999"     "682.6" "-999"  "-999"  "-999"
#> [264,] "FNLN2"        "2"   "Go"   "-999"     "385.3" "-999"  "-999"  "-999"
#> [265,] "FNLN2"        "2"   "Go"   "-999"     "335.8" "-999"  "-999"  "-999"
#> [266,] "FNLN2"        "2"   "Go"   "-999"     "426.8" "-999"  "-999"  "-999"
#> [267,] "FNLN2"        "2"   "Go"   "-999"     "538"   "-999"  "-999"  "-999"
#> [268,] "FNLN2"        "2"   "Go"   "-999"     "615"   "-999"  "-999"  "-999"
#> [269,] "FNLN2"        "2"   "Go"   "-999"     "219.5" "-999"  "-999"  "-999"
#> [270,] "FNLN2"        "2"   "Go"   "-999"     "528"   "-999"  "-999"  "-999"
#> [271,] "FNLN2"        "2"   "Go"   "-999"     "591.4" "-999"  "-999"  "-999"
#> [272,] "FNLN2"        "2"   "Go"   "-999"     "514.3" "-999"  "-999"  "-999"
#> [273,] "FNLN2"        "2"   "Go"   "-999"     "619.7" "-999"  "-999"  "-999"
#> [274,] "FNLN2"        "2"   "Go"   "-999"     "462.4" "-999"  "-999"  "-999"
#> [275,] "FNLN2"        "2"   "Go"   "-999"     "393.3" "-999"  "-999"  "-999"
#> [276,] "FNLN2"        "2"   "Go"   "-999"     "552.5" "-999"  "-999"  "-999"
#> [277,] "FNLN2"        "2"   "Go"   "-999"     "597.9" "-999"  "-999"  "-999"
#> [278,] "FNLN2"        "2"   "Go"   "-999"     "558.3" "-999"  "-999"  "-999"
#> [279,] "FNLN2"        "2"   "Go"   "-999"     "741.3" "-999"  "-999"  "-999"
#> [280,] "FNLN2"        "2"   "Go"   "-999"     "336.7" "-999"  "-999"  "-999"
#> [281,] "FNLN2"        "2"   "Go"   "-999"     "687.7" "-999"  "-999"  "-999"
#> [282,] "FNLN2"        "2"   "Go"   "-999"     "657.4" "-999"  "-999"  "-999"
#> [283,] "FNLN2"        "2"   "Go"   "-999"     "470.6" "-999"  "-999"  "-999"
#> [284,] "FNLN2"        "2"   "Go"   "-999"     "533.5" "-999"  "-999"  "-999"
#> [285,] "FNLN2"        "2"   "Go"   "-999"     "564.5" "-999"  "-999"  "-999"
#> [286,] "FNLN2"        "2"   "Go"   "-999"     "865.7" "-999"  "-999"  "-999"
#> [287,] "FNLN2"        "2"   "Go"   "-999"     "451.9" "-999"  "-999"  "-999"
#> [288,] "FNLN2"        "2"   "Go"   "-999"     "617.5" "-999"  "-999"  "-999"
#> [289,] "FNLN2"        "2"   "Go"   "-999"     "609.3" "-999"  "-999"  "-999"
#> [290,] "FNLN2"        "2"   "Go"   "-999"     "645.8" "-999"  "-999"  "-999"
#> [291,] "FNLN2"        "2"   "Go"   "-999"     "592"   "-999"  "-999"  "-999"
#> [292,] "FNLN2"        "2"   "Go"   "-999"     "565.8" "-999"  "-999"  "-999"
#> [293,] "FNLN2"        "2"   "Go"   "-999"     "583.2" "-999"  "-999"  "-999"
#> [294,] "FNLN2"        "2"   "Go"   "-999"     "632"   "-999"  "-999"  "-999"
#> [295,] "FNLN2"        "2"   "Go"   "-999"     "482.4" "-999"  "-999"  "-999"
#> [296,] "FNLN2"        "2"   "Go"   "-999"     "368.1" "-999"  "-999"  "-999"
#> [297,] "FNLN2"        "2"   "Go"   "-999"     "573.5" "-999"  "-999"  "-999"
#> [298,] "FNLN2"        "2"   "Go"   "-999"     "395.1" "-999"  "-999"  "-999"
#> [299,] "FNLN2"        "2"   "Go"   "-999"     "443"   "-999"  "-999"  "-999"
#> [300,] "FNLN2"        "2"   "Go"   "-999"     "586.4" "-999"  "-999"  "-999"