Thank you very much. This was, was i needed. Unfortunately I have one futher
problem with this Code. I don't only need the coverage probability for one
but for a range of different odds ratios. (for example [1;30]). I tried it
with a loop but I get an error. I think again, that I'm almost there but
having a little mistake. The complete code is:
#setting values
n1 <- 10
n2 <- 10
y <- 100
alpha <- 1
z<-1.96
# creating 2x2 table
for (i in 1:30)
{
theta <- i
x1 <- exp(alpha +theta)/ (1+ exp(alpha +theta))
x2 <- exp(alpha)/ (1+ exp(alpha))
n11 <- rbinom(y, 10, x1)
n12 <- n1 - n11
n21 <- rbinom(y, 10, x2)
n22 <- n2 - n21
# upper and lower limit gart interval
gartu <-function(z,d,e, f, g){log(((d+.5)*(g+.5))/((e+.5)*(f+.5)))+
z*sqrt(1/(d+.5)+1/(e+.5)+1/(f+.5)+1/(g+.5))}
gartl <-function(z,d,e, f, g){log(((d+.5)*(g+.5))/((e+.5)*(f+.5)))-
z*sqrt(1/(d+.5)+1/(e+.5)+1/(f+.5)+1/(g+.5))}
u <- gartu(z, n11[i],n22[i],n12[i],n21[i])
l <- gartl(z, n11[i],n22[i],n12[i],n21[i])
foo <- function(theta, u, l) mean(theta >= l & theta <= u, na.rm = TRUE)
foo(theta, u, l)
}
--
View this message in context:
http://r.789695.n4.nabble.com/Coverage-Probability-tp4485511p4485865.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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.
