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 <corymred...@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]] _______________________________________________ R-sig-ecology mailing list R-sig-ecology@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology