if (FALSE) {
data(Skabholmen)
Y <- Skabbholmen$Y
X <- Skabbholmen$X
model <- gllvm(y = Y, X = X, studyDesign = X,
lv.formula = ~scale(Elevation) + I(Year-1978),
num.RR = 2,
family = "ordinal",
zeta.struc="common",
row.eff=~(1|transectID))
}
Run the code above in your browser using DataLab