data(trichoptera)
trichoptera <- prepare_data(trichoptera$Abundance, trichoptera$Covariate)
myPLNLDA <- PLNLDA(Abundance ~ 1, grouping = Group, data = trichoptera)
if (FALSE) {
plot(myPLNLDA, map = "individual", nb_axes = 2)
}
Run the code above in your browser using DataLab