Constructing a ZSM SAD in R
Dear Roy, I haven't done this, but I would start with the function zsm in the package untb by Robin Hankin. See also the function etienne. Yours, P?ter P?ter S?lymos, PhD Department of Mathematical and Statistical Sciences University of Alberta Edmonton, Alberta, T6G 2G1 Canada
On Sun, Jan 11, 2009 at 12:33 PM, Roy Martin <royworth at gmail.com> wrote:
Hi all,
I was hoping someone on the board would have experience constructing the
zero-sum multinomial (ZSM) species abundance distribution (SAD) in R.
Unfortunately, my expertise in probability theory and R programming appear
limited, which has left me largely confused about how to construct the
expected SAD dataset in R using the maximum likelihood estimates of 'theta'
and 'm' that I have obtained for my dataset using the methods and PARI/GP
code provided in Etienne (2005). If anyone has any experience with this in
R, I would be extremely grateful for any suggestions.
Sincerely,
Roy Martin
PhD Candidate
Wildlife and Fisheries Resources
West Virginia University
Morgantown, WV 26506
rmarti10 at mix.wvu.edu
royworth at gmail.com
Etienne, R.S. 2005. A new sampling formula for neutral biodiversity. Ecology
Letters 8: 253-260.
[[alternative HTML version deleted]]
_______________________________________________ R-sig-ecology mailing list R-sig-ecology at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology