## symmetric CoCA
data(beetles)
data(plants)
## log transform the bettle data
beetles <- log(beetles + 1)
## fit the model
bp.sym <- coca(beetles ~ ., data = plants, method = "symmetric")
## draw a screeplot of the eignevalues
screeplot(bp.sym)
Run the code above in your browser using DataLab