Trouble plotting a maximum likelihood model(s) to rank abundance curves in fitsad
Dear Cory, Please see section 5 of the package vignette to add lines of model predictions to rad and octave plots. In a nutshell you need objects with the results of each fit and then you can use these objects to get predicted values with 'radpred' and 'octavpred'. Please also note that the data input to 'fitsad' and 'fitrad' are numeric vectors, not matrices. Best wishes Paulo 2018-09-24 13:22 GMT-03:00 Cory Redman <corymredman at gmail.com>:
I have been trying to create a plot showing maximum likelihood models fitted to a rank abundance curve in the 'fitsad' package in R. See attached pic for what I am aiming for. Y-axis can be in log scale or not. I can't even create a plot that shows the results of the model. I have looked over the online literature, but can't figure it out. Please help. Cheers, Cory My code:
library(sads) Unit2 = matrix(c(529,230,129,108,98,62,40,23,18,13,9,5,3)) plot(Unit2) Unit2.ln <- fitsad(Unit2, 'lnorm') plot(Unit2.ln)
Error in (function (classes, fdef, mtable) : unable to find an inherited method for function ?octav? for signature ?"matrix"?
*Paulo In?cio de Knegt L?pez de Prado *Depto. de Ecologia - Instituto de Bioci?ncias - USP Associate Editor, Ecological Complexity <https://www.journals.elsevier.com/ecological-complexity/editorial-board> Recommender, Peer Community in Ecology <https://ecology.peercommunityin.org/> Address: Rua do Mat?o, travessa 14, n? 321 Cid. Universit?ria, S?o Paulo - SP - Brasil CEP 05508-900 +55-11-30910963 (sala) +55-11-30917600 (Secretaria) [[alternative HTML version deleted]]