使用聚合计数数据拟合R中的泊松GLM

时间:2015-12-15 01:24:25

标签: r glm poisson

我有一个关于澳大利亚昆士兰沿岸各个地方报告的搁浅龟数量的数据集。我想知道的是在这些地点没有报告的搁浅龟的数量。为了估计这个数字,我收集了关于将乌龟报告到搁浅地点的频率的数据;也就是说,沿海岸约20个点,单个乌龟搁浅的次数是多少次?所以我有计数数据,表明一次,两次或三次或更多次报告到绞线位置的海龟数量。最后,我想将这些数据与协变量联系起来,例如当地人口密度和距离最近的道路的距离,以便预测零报告"其他沿海地区的发病率也是如此。

我的数据看起来应该是这样的,然后:

loc<-c("A","B","C")
rep1<-c(51,24,10)
rep2<-c(4,8,3)
rep3ormore<-c(2,1,0)
pop<-c(50,1000,100)
turtle- cbind.data.frame(loc, rep1, rep2, rep3ormore, pop)

还有其他可能的协变量,但我现在要保持简单!我认为应该可以使用Poisson分布来完成,但是我无法绕过如何做到这一点。

此外,在某些情况下,我没有确切报告过的海龟的数字,但我有类别; 4-6,7-10,&gt; 10等。如果有可能对这种可能性进行建模,那也会很棒!

0 个答案:

没有答案