拟合泊松分布

Question

我无法计算泊松分布的最大似然估计和BIC。我能够得到直方图，但无法在其上叠加核密度估计。

你能告诉我哪里出错了吗？

x.pois<-rpois(Y1, 20)
hist(x.pois, breaks=100,freq=FALSE)
lines(density(Y1, bw=0.8), col="red")
library(MASS)
fitdistr(Y1,densfun="pois")
my.mle<-fitdistr(Y1, densfun="poison")
print(my.mle)
BIC(my.mle)

Answer 1

你需要（1）正确拼写“泊松”; （2）使用x.pois（泊松样本），而不是Y1（根据您的代码示例，它应该是您尝试采样的点数）。

请注意，离散分布的核密度估计和直方图不一定非常有意义。

Y1 <- 100
set.seed(101) ## for reproducibility
x.pois<-rpois(Y1, 20)
hist(x.pois, breaks=100,freq=FALSE)
lines(density(x.pois, bw=0.8), col="red")
library(MASS)
(my.mle<-fitdistr(x.pois, densfun="poisson"))
##      lambda  
##   20.6700000 
##  ( 0.4546427)
BIC(my.mle)
## [1] 572.7861

更新：your other question清楚地表明Y1确实是您的样本，在这种情况下，整个rpois() - 抽样事物只是一个红色的鲱鱼。在这种情况下，您应该省略前三行，并在上面的代码中将Y1替换为x.pois。