用汽车包拟合二项分布答案

【问题标题】：Fitting Binomial distribution with car package用汽车包拟合二项分布
【发布时间】：2017-01-27 16:40:53
【问题描述】：

为了拟合我的数据集的多个概率分布，我尝试了使用 car 和 MASS 包的多个分布，例如泊松分布：

require(car); require (MASS)

#Data set
X1 <- rnorm(100, 175, 7)
Y1 <- rbinom(100,1, 0.35)
##Fitting Poisson:
poisson <- fitdistr(Y1, "Poisson")
qqp(Y1, "pois", poisson$estimate)

但是，如果我的选项是二项分布，MASS 包不支持这种分布，那么我使用fitdistrplus 包：

require(fitdistrplus)
fitBinom=fitdist(data=Y1, dist="binom", fix.arg=list(size=100), start=list(prob=0.35))

当我尝试时：

qqp(Y1, "Binomial", fitBinom$estimate)

不适用于car 包的qqp 功能。还有其他适合二项分布的选项，例如qqp 函数？

谢谢，

【问题讨论】：

标签： r r-car probability-distribution

【解决方案1】：

qqp 实际上适用于 fitdistrplus。您只需添加二项分布的所有参数。您还需要将“Binomial”替换为“binom”

fitBinom <- fitdist(data=Y1, dist="binom", fix.arg=list(size=100), start=list(prob=0.35))
qqp(Y1, "binom", size=100, fitBinom$estimate)

【讨论】：