# NOT RUN {
library(MPDiR) # contains the Vernier data
fit <- quickpsy(Vernier, Phaseshift, NumUpward, N,
grouping = .(Direction, WaveForm, TempFreq), B = 20)
loglikssaturated(fit)
# }
Run the code above in your browser using DataLab