Great, that helps. With all those variables, your function worked for
me. I did slightly edit it for testing, mostly just adding spaces to
make it easier for me to read.
N <- 10; q <- 8; k <- q-2; L <- 5; n <- 2 * k + L + 1
p_0 <- c(0.2, 0.4, 0.4)
f_1 <- qf(0.05, 2, 2)
f_2 <- qf(0.95, 2, 2)
p_11 <- seq(0, 0.3, 0.1)
p_12 <- seq(0.1, 0.4, 0.1)
guete <- function(p_11, p_12) {
set.seed(1000)
S_vek <- matrix(0, nrow = N, ncol = 1)
zz <- cbind(p_11, p_12, (1 - p_11 - p_12))
for(i in 1:N) {
X_0 <- rmultinom(q - 1, size = 1, prob = p_0)
X_1 <- rmultinom(n - q + 1, size = 1, prob = zz)
N_0 <- apply(X_0[ ,(n - 2 * k - L + 1):(n - k - L)], 1, sum)
N_1 <- apply(X_1[ ,(n - q - k + 2):(n - q + 1)], 1, sum)
S_vek[i] <- ( (sum( ((N_1 - k * zz)^2) / k * zz)) /
(sum( ((N_0 - k * p_0)^2) / k * p_0))
) - 1
}
1 - mean(f_1 <= S_vek & S_vek <= f_2)
}
guete(p_11 = p_11, p_12 = p_12)
For your second problem:
p_11=seq(0,1,0.1)
p_12=seq(0,1,0.1)
then i get also this error message:
Error in rmultinom(n - q + 1, size = 1, prob = rbind(p_11, p_12, (1 - :
non-positive probability
the probabilities for rmultinom need to be non-negative. But when you
do seq(0, 1, .1), you get some 1s in p_11 and p_12. This means part
of your definition for the probs argument is (1 - 1 - 1) = -1. So you
need to make sure that that maximum value of p_11 + p_12 is <= 1 (if
you are going to subtract both from 1 and use that as a probability).
Hope that helps,
Josh
On Wed, Sep 8, 2010 at 3:24 PM, tuggi <tugbaguec...@web.de> wrote:
>
> hello,
>
> thank you for your answer. here the parameter
>
> N=10
>
>
> q=8
>
>
>
> k=q-2
>
>
>
> L=5
>
>
>
> n=2*k+L+1
>
>
> p_0=c(0.2,0.4,0.4)
>
> f_1=qf(0.05,2,2)
>
>
> f_2=qf(0.95,2,2)
>
> tuggi
> --
> View this message in context:
> http://r.789695.n4.nabble.com/problem-with-outer-tp2532074p2532132.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.
>
--
Joshua Wiley
Ph.D. Student, Health Psychology
University of California, Los Angeles
http://www.joshuawiley.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.
