# NOT RUN {
# Loading of the verbal data
data(verbal)
attach(verbal)
# Excluding the "Anger" variable
verbal <- verbal[colnames(verbal)!="Anger"]
# Keeping the first 5 items and the first 50 subjects
# (this is an artificial simplification to reduce the computational time)
# Sixth column holds the group membership
verbal <- verbal[1:50, c(1:5, 25)]
# Likelihood-ratio statistics
LRT(verbal[,1:5], verbal[,6])
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab