[R-sig-eco] Trouble plotting a maximum likelihood model(s) to rank abundance curves in fitsad

[Paulo Inácio Prado]

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 using 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]]