本笔记中原始数据及代码均来源于李东风先生的R语言教程,在此对李东风先生的无私分享表示感谢。
vcd包提供了一个goodfit函数, 可以用来拟合指定的某种理论分布(包括泊松、二项、负二项分布), 并检验服从该理论分布的零假设。 例如,我们生成一组速率参数为2的泊松随机数, 检验其分布是否泊松分布:
library(vcd)
## 载入需要的程辑包:grid## Warning message:## 程辑包‘vcd’是用R版本4.0.2 来建造的
set.seed(101)datax 100, summary(goodfit(datax, "poisson"))
## Goodness-of-fit test for poisson distribution## ## X^2 df P(> X^2)## Likelihood Ratio 4.289456 5 0.5085374
检验其是否服从二项分布,取二项分布试验数为10:
summary(goodfit(datax, "binomial", par = list(size = 10)))