ポアソン定数の信頼区間

こちらから。
　

poisson.conf <- function(x, conf=0.95){
	N2 <- length(x)*2
	df <- 2*sum(x)+2
	alpha2 <- (1-conf)/2
	return(qchisq(c(alpha2, 1-alpha2), df)/N2)
}
n<- 100
poisCI<- matrix(0,n+1,2)
dimnames(poisCI)<- list(0:n,c("lower","upper"))
for(i in 0:n){
	poisCI[i+1,]<- poisson.conf(i)
}