od <- options(digits = 4)
## 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")
## weights == R[0]
weights(bp.sym)
options(od)
Run the code above in your browser using DataLab