## analyze in respect to amplification
## efficiency outliers
ml1 <- modlist(reps)
res1 <- KOD(ml1, check = "uni2")
## which runs are outliers?
outl <- is.outlier(res1)
outl
which(outl)
## test for sigmoidal outliers
## with the 'testdat' dataset
ml2 <- modlist(testdat, model = l5, check = "uni2")
is.outlier(ml2)
Run the code above in your browser using DataLab