[R-sig-eco] Fit log-series to species-rank abundance

[Paulo Prado]

Dear Guillem,

May 1975 (in Patterns of Species Abundance and Diversity. In M. L. Cody & J.
M. Diamond (Eds.)) gives the recipe to fit expected curves of a SAD model in
rank-abundance plots. The trick is to find the inverse function of the CDF
of the SAD model. A simple solution in R is to calculate the cumulative
number of species predicted by the logseries and to use this values as
independente variables for the function splinefun, to create a function that
plots the expected curve of the model in the rank-abundance plot. Please
find attached a code that seems to do the job at least for visual comparison
of models. There are also the PDF function  for the logseries, and a
function to fit it that returns the log-likelihood and the AIC of the model.

Hope it helps

Best Wishes

Paulo Prado

http://r-sig-ecology.471788.n2.nabble.com/file/n5734514/logseries.r
logseries.r 
-- 
View this message in context: http://r-sig-ecology.471788.n2.nabble.com/Fit-log-series-to-species-rank-abundance-tp5729089p5734514.html
Sent from the r-sig-ecology mailing list archive at Nabble.com.