# NOT RUN {
data <- data.frame(y=ChickWeight$weight,g=ChickWeight$Diet)
aggregate(data$y,by=list(data$g),mean)
#create matrices HR: g4>(g1,g2,g3)
HR <- create_matrices(varnames = c("g1","g2","g3","g4"),
hyp = "g4>(g1,g2,g3)")
Fbar.ineq(data,Amat=HR$Amat)
# }
Run the code above in your browser using DataLab