Human / S. meleagridis strain I data [1]
| Dose |
Infected |
Non-infected |
Total |
| 12000 |
3 |
3 |
6 |
| 24000 |
4 |
2 |
6 |
| 52000 |
3 |
3 |
6 |
| 96000 |
3 |
3 |
6 |
| 155000 |
5 |
1 |
6 |
| 3E+05 |
6 |
0 |
6 |
| 720000 |
4 |
1 |
5 |
| 1150000 |
6 |
0 |
6 |
| 5.5E+06 |
5 |
1 |
6 |
| 2.4E+07 |
5 |
0 |
5 |
| 5E+07 |
6 |
0 |
6 |
|
Goodness of fit and model selection
| Model |
Deviance |
Δ |
Degrees
of freedom |
χ20.95,1
p-value |
χ20.95,m-k
p-value |
| Exponential |
71.8 |
63.8 |
10 |
3.84
1.33e-15 |
18.3
1.98e-11 |
| Beta Poisson |
7.99 |
9 |
16.9
0.535 |
| Beta-Poisson fits better than exponential; cannot reject good fit for beta-Poisson. |
|
Optimized parameters for the beta-Poisson model, from 10000 bootstrap iterations
| Parameter |
MLE estimate |
Percentiles |
| 0.5% |
2.5% |
5% |
95% |
97.5% |
99.5% |
| α |
3.89E-01 |
1.23E-01 |
1.74E-01 |
2.01E-01 |
1.12E+00 |
1.91E+00 |
3.82E+02 |
| N50 |
1.68E+04 |
7.98E+01 |
1.08E+03 |
2.27E+03 |
4.78E+04 |
5.71E+04 |
7.49E+04 |
|
Parameter scatter plot for beta Poisson model ellipses signify the 0.9, 0.95 and 0.99 confidence of the parameters.
beta Poisson model plot, with confidence bounds around optimized model
References
QMRA