Hi hubinho,
You are almost there. Try this slightly modification of your function:
# theta, u and l are vectors of the same length
foo <- function(theta, u, l) mean(theta >= l & theta <= u, na.rm =
TRUE)
foo(theta, u, l)
HTH,
Jorge.-
On Mon, Mar 19, 2012 at 12:55 PM, hubinho <> wrote:
> Hello.
>
> I'm allready this far. I have a function which is calculating the lower
(l)
> and upper (u) limit for a confidence interval for the odds ratio.
>
> For example for 5 simulated 2x2 tables the upper and lower limits are:
>
> > u
> [1] 2.496141 7.436524 8.209161 4.313587 3.318612
> > l
> [1] -0.9718608 1.1000713 1.5715373 0.1135158 -0.2700517
>
> With (l[1]; u[1]) being the confidence interval for the odds ratio for the
> first simulated table and so on.
>
> Now I want to compute the coverage probability. For that I've created a
> function which is return 1 if the odds ratio is in the interval and 0 if it
> isn't.
>
> cover <- function(theta, u, l){
> if(theta >= l && theta <= u){z=1}
> if(theta < l || theta > u){z=0}; return(z)
> }
>
> This works but unfortunately not if I want to summarize the function and
> divide it with the sample size to get the coverage probability.
>
> I tried it this way
>
> for(for(x in 1:5) {a = (sum(cover(theta, u[x], l[x]))/5; return(a)}
>
> Maybe someone can help me. Thank you
>
>
>
> --
> View this message in context:
> http://r.789695.n4.nabble.com/Coverage-Probability-tp4485511p4485511.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.
>
[[alternative HTML version deleted]]