I am curious about performing poisson test where I have uncertainty about my count. For example, I expect to see 15 bunny rabbits per hike. On a given hike, I positively identify 19 bunny rabbits, and I think I see 4 more bunny rabbits, but I can't be sure. How do I include my maybe bunny rabbits into my poisson test? Is there a more difficult to pronounce test I need to use in this case?
I'm not sure this is the right approach, and I'm not sure what to do with the resulting standard error.
> install.packages(Amelia); > library(Amelia); > q <- data.frame("bunny-count"=c(20,19,23,22,19,22)) > se <- data.frame("se"=c(0.0625,0.0625,0.0625,0.0625,0.0625,0.0625)) > mi.meld(q, se) $q.mi bunny.count [1,] 20.83333 $se.mi se [1,] 1.861456