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@mix.wvu.edu
royworth@gmail.com

Etienne, R.S. 2005. A new sampling formula for neutral biodiversity. Ecology
Letters 8: 253-260.

	[[alternative HTML version deleted]]

