if (FALSE) {
# Load and prepare data
data <- Exam
classes <- c(1, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.7, 8.5, Inf)
data$examsc.class <- cut(data$examsc, classes)
# Run model with random intercept and default settings
model1 <- semLme(
formula = examsc.class ~ standLRT + schavg + (1 | school),
data = data, classes = classes
)
summary(model1)
}
# \dontshow{
# Load and prepare data
data <- Exam
classes <- c(1, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.7, 8.5, Inf)
data$examsc.class <- cut(data$examsc, classes)
# Run model with random intercept and default settings
model1 <- semLme(
formula = examsc.class ~ standLRT + schavg + (1 | school),
data = data, classes = classes, burnin = 4, samples = 10
)
summary(model1)
# }
Run the code above in your browser using DataLab