我有以下数据:
data= c(0,1,0,11,2,0,3,0,0,2, 1,3,1, 0,1,0,0,0,2,3,0,0,0,8,1,1,1,0,1,1,
2,7,0,0,0,5,2,3, 6, 1, 1, 5 , 2 , 9 , 0, 0, 1, 21, 16, 2, 9, 6,
25, 2, 1, 12, 16, 14, 15, 15, 6, 1 ,12 ,12, 13 , 5 ,5 , 6 , 4 , 7 ,11 ,
8 , 4 , 5 , 8, 3 , 8 , 4 , 7 , 4 , 7 , 2 , 5, 6 , 4 , 5 , 1 , 0, 8 , 5, 6 ,
8 , 9, 8, 9 , 7 , 7 , 9 , 8 , 9 ,4 ,4 , 7 ,13 , 9 ,13 ,12, 10 , 9 , 8, 7, 11,
5 , 5 , 0, 1, 33 , 4 ,22, 19, 22, 9 , 5,4 ,17 , 7 , 7 , 4, 5, 3 , 0 , 0 , 9,
3, 0 , 0 ,36 ,40, 5 , 4, 0, 11 , 0 ,7 ,5 ,25, 39, 26 , 4 ,20, 12, 4, 17 , 3, 22,
12 ,14 , 8 , 9, 11, 7 ,11, 10 , 9 ,16, 6, 24, 8 , 5 , 6, 14 , 3 ,9 , 4, 1 ,20 , 0 ,
1 , 7 ,9 ,0 ,12, 2, 29 ,56, 16 , 8 ,28,0 ,19, 25 ,35, 87, 56, 66, 60 ,58, 14, 10 ,12,
13 ,13 ,34, 26, 18 ,13, 22, 13 ,12 ,15 ,41 ,11 ,11 ,11 , 5 , 6, 7 , 8 , 8 ,17, 16,
12, 21, 38 ,34 ,10 ,77 ,41 , 7 ,12 , 1, 16, 20, 8 , 5, 2, 20 , 7 ,16, 12 , 6, 10, 31, 12)
我使用了 fitdistrplus 包通过泊松分布拟合了这些数据。
x=fitdist(data,"pois")
一旦我在 x 上使用绘图函数:
plot(x)
默认情况下构建两个并排图。一个有经验和理论分布图,另一个有 CDF。
数据集有 (29/244)=12% 零值,因此我关心的是我想使用零截断泊松分布拟合数据,并希望拥有使用 fitdistrplus 生成的相同类型的图形包 [ 正如我在前面部分提到的: plot(x) ]
感谢任何帮助。谢谢
答案 0 :(得分:1)
actuar
包包含一个零截断的泊松分布 dztpois
。您需要指定一个起始值。
library(actuar)
d2 <- data[data>0]
xz=fitdist(d2,"ztpois", start=list(lambda=mean(d2)))
plot(xz)
如评论中所述,在泊松拟合和零截断泊松拟合之间对已丢弃零的数据进行统计比较是没有意义的。