library("PlackettLuce")
example("beans", package = "PlackettLuce")
G = group(R, rep(seq_len(nrow(beans)), 4))
d = cbind(G, beans)
split = ifelse(d$maxTN < 18.7175, TRUE, FALSE)
likelihood_ratio(G, split)
mod = PlackettLuce(G)
anova(mod)
Run the code above in your browser using DataLab