# NOT RUN {
### get allometry model for each species
CSP <- plyr::dlply(allometry, c("spp"), bCM)
CSP.coef <- plyr::ldply(CSP, stats::coef)
### add lambda value
CSP.coef$lam <- plyr::ddply(allometry, c("spp"), function(df)
bCM(df, lam.only = TRUE))[, "V1"]
# }
Run the code above in your browser using DataLab