# NOT RUN {
## 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")
## fitted values
bp.fit <- fitted(bp.sym)
# }
# NOT RUN {
bp.fit
# }
# NOT RUN {
## fitted values for beetles only
beetle.fit <- fitted(bp.sym, which = "y1")
# }
Run the code above in your browser using DataLab