Re: [R] Confidence Intervals for Poisson

On Jan 15, 2009, at 1:10 AM, John Kerpel wrote: Hi folks! I run the following code to get a CI for a Poisson with lambda=12.73 library(MASS) set.seed(125) x <- rpois(100,12.73) lambda_hat<-fitdistr(x, dpois, list(lambda=12))$estimate #Confidence Intervals - Normal Approx. alpha<-c(.05,.02

[R] Confidence Intervals for Poisson

Hi folks! I run the following code to get a CI for a Poisson with lambda=12.73 library(MASS) set.seed(125) x <- rpois(100,12.73) lambda_hat<-fitdistr(x, dpois, list(lambda=12))$estimate #Confidence Intervals - Normal Approx. alpha<-c(.05,.025,.01) for(n in 1:length(alpha)) { LowerCI<-mean(x