Test for Poisson distribution

20 次查看（过去 30 天）
SBinary
SBinary 2021-4-9
编辑： Abdolkarim Mohammadi 2021-4-9
I have a vector of observations (length: 1195).
Previous studies have found the data to have a Poisson distribution, and when I plot the data it closely resembles what I would expect to see from a Poisson distribution (see image below).
And the fit (and quality metrics) looks like the following,
Distribution: Poisson
Log likelihood: -5094.64
Domain: 0 <= y < Inf
Mean: 5.63264
Variance: 5.63264
Parameter Estimate Std. Err.
lambda 5.63264 0.068655
Estimated covariance of parameter estimates:
lambda
lambda 0.0047135
However, when I test the hypothesis using a Cho-squared goodness of fit test (chi2gof) as detailed here, Chi-square goodness-of-fit test - MATLAB chi2gof - MathWorks Nordic, the results indicate that the hypothesis should be rejected (at alpha 0.05).
I am not confident I have not made a mistake, and would really appreciate some input. The following is the code for the chi2gof test,
n = sum(x);
bins = 0:1194;
pd = fitdist(bins','Poisson','Frequency',x');
expCounts = n * pdf(pd,bins);
[h,p,st] = chi2gof(bins,'Ctrs',bins,...
'Frequency',x, ...
'Expected',expCounts,...
'NParams',1)
Here is a sample of x (1:100),
x = [3,1,3,4,0,0,0,4,10,2,0,4,6,2,2,2,5,2,0,4,4,2,1,2,3,3,6,0,3,2,0,0,2,4,4,3,0,0,2,3,0,0,0,0,0,0,0,0,1,1,3,2,2,2,1,8,2,4,6,0,3,1,3,1,3,1,1,4,1,2,9,6,6,0,11,3,0,1,2,2,3,0,6,0,3,0,8,5,8,6,7,5,5,0,8,1,3,0,0,8]

请先登录，再回答此问题。

回答（1 个）

Abdolkarim Mohammadi
编辑：Abdolkarim Mohammadi 2021-4-9
Maybe the problem is that you have created too many bins; one bin for each observation. That would not yield good results. Consider lowering the number of bins.
In addition, x is your data and NOT the frequency. You can use the histcounts() function to get the desired bins.

类别

在 Help Center 和 File Exchange 中查找有关 Hypothesis Tests 的更多信息

标签

产品


版本

R2018a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by