# 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])Run the code above in your browser using DataLab