statsmodels.stats.rates.test_poisson

statsmodels.stats.rates.test_poisson(count, nobs, value, method=None, alternative='two-sided', dispersion=1)[source]

Test for one sample poisson mean or rate

Parameters:
count : array_like

Observed count, number of events.

nobs : arrat_like

Currently this is total exposure time of the count variable. This will likely change.

value : float, array_like

This is the value of poisson rate under the null hypothesis.

method : str

Method to use for confidence interval. This is required, there is currently no default method. See Notes for available methods.

alternative : {'two-sided', 'smaller', 'larger'}

alternative hypothesis, which can be two-sided or either one of the one-sided tests.

dispersion : float

Dispersion scale coefficient for Poisson QMLE. Default is that the data follows a Poisson distribution. Dispersion different from 1 correspond to excess-dispersion in Poisson quasi-likelihood (GLM). Dispersion coeffficient different from one is currently only used in wald and score method.

Return type:

HolderTuple instance with test statistic, pvalue and other attributes.

Notes

The implementatio of the hypothesis test is mainly based on the references for the confidence interval, see confint_poisson.

Available methods are:

  • “score” : based on score test, uses variance under null value

  • “wald” : based on wald test, uses variance base on estimated rate.

  • “waldccv” : based on wald test with 0.5 count added to variance computation. This does not use continuity correction for the center of the confidence interval.

  • “exact-c” central confidence interval based on gamma distribution

  • “midp-c” : based on midp correction of central exact confidence interval. this uses numerical inversion of the test function. not vectorized.

  • “sqrt” : based on square root transformed counts

  • “sqrt-a” based on Anscombe square root transformation of counts + 3/8.

See also

confint_poisson