simdata1 <- sim.ny(par = c(30, 200, 250, 350), esterrorpar = c(-2.67, 0.0094,
0.0299, 0.0020, 1.13), c = 4, k = 12, N = 10000, pr = TRUE, seed = 0)
simdata2 <- sim.ny(par = c(30, 200, 250, 350), esterrorpar = c(-2.67, 0.0094,
0.0299, 0.0020, 1.13), c = 4, k = 12, N = 10000, pr = TRUE, seed =
12345)
simcorrectRT1 <- subset(simdata1, simdata1[,2] == TRUE, select = 1)
simcorrectRT2 <- subset(simdata2, simdata2[,2] == TRUE, select = 1)
LqDist(simcorrectRT1, simcorrectRT2, q=1)
Run the code above in your browser using DataLab