Many thanks to Dr Byrne, the explicit expression of the mean
hinted at the correct direction, the error I made was confusing
natural logarithm and base 10 logarithm!
I apologize that I didn't post the complete example data,
which follows:
No. of species Theoretical frequency Observed frequency
1 135.05 118
2 67.33 74
3 44.75 44
4 33.46 24
5 26.69 29
6 22.17 22
7 18.95 20
8 16.53 19
9 14.65 20
10 13.14 15
11 11.91 12
12 10.89 14
13 10.02 6
14 9.28 12
15 8.63 6
16 8.07 9
17 7.57 9
18 7.13 6
19 6.74 10
20 6.38 10
21 6.06 11
22 5.77 5
23 5.5 3
24 5.25 3
(from Fisher 1943, Kendall-Stuart 1977)
With thanks,
Vincent Vinh-Hung
Oncologisch Centrum, AZ-VUB
101 Laarbeeklaan
B-1090 Jette
Graeme Byrne wrote:
(...)
>
> You will need a value of q first. This will either be estimated from the raw
> data or assumed by some hypothesis. Once you have this just plug in the
> value of r you want and multiply the resulting probability by the sum of
> the observed frequencies.
>
> You might also be able to use the theorectical mean q/((q - 1 )*Log[1 - q])
> to estimate q by equating it to the sample mean and solving for q.
>
>
> --
> Dr Graeme Byrne
> La Trobe University, Bendigo
> PO Box 199, Bendigo, 3552
> Phone: 61 3 5444 7263
> Fax: 61 3 5444 7998
> [EMAIL PROTECTED]
===========================================================================
This list is open to everyone. Occasionally, less thoughtful
people send inappropriate messages. Please DO NOT COMPLAIN TO
THE POSTMASTER about these messages because the postmaster has no
way of controlling them, and excessive complaints will result in
termination of the list.
For information about this list, including information about the
problem of inappropriate messages and information about how to
unsubscribe, please see the web page at
http://jse.stat.ncsu.edu/
===========================================================================
