>
> Hi All:
>
> I am using R to calculate exact 95% confidence interval using Clopper
> Pearson method. I am using the following code but it seems to get into a
> loop and not get out of it, it goes on forever although I am looping it
> only 10 times across 63 sites with 10 observations per site. I was hoping
> to get some help.
>
> Thanks
>
Anamika

> set.seed(111)  #so that you can regenerate the same values
> # another time if you need to
> k<-63
> n<-10
> x<-NULL
> eta<-rnorm(63,-0.085,.990)  # generating data using Binomial Logit Normal
> p<-exp(eta)/(1+exp(eta))
> for(i in 1:k)
> x<-cbind(x,rbinom(10,n,p[i]))
>
> rate<-t(t(x)/n)
>
> # Exact Confidence Interval
>
> for (i in 1:10){
>
> for (j in 1:63)
> {
>
> if (x==0)
> {
> l_cl_exact<-0
> u_cl_exact<-qbeta(.975,x+1,n-x)
> }
> else if (x==n)
> {
> l_cl_exact<-qbeta(.025,x,n-x+1)
> u_cl_exact<-1
> }
> else
> l_cl_exact<-qbeta(.025,x,n-x+1)
> u_cl_exact<-qbeta(.975,x+1,n-x)
>
> }
> }_____________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html>
> and provide commented, minimal, self-contained, reproducible code.
>

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to