I have never seen 'break' used in an 'ifelse'; you probably meant to
use an 'if' statement there.
On Thu, Mar 1, 2012 at 8:24 AM, Tsai, Pei-Chien
<pei-chien.t...@kcl.ac.uk> wrote:
> Dear R helpers,
>
> I have some difficulties in using 'break' function with loop, and the
> followings are my script. What I try to do is (1) permute 'or' first; (2)
> doing t-test if this 'or' pass criteria 1 (k=1); (3) end the loop when I get
> 10 permutations; (4) redo everything again but this time use criteria 2 (k=2)
> (I have more criteria 1:n).
>
> Somehow using my script, the final dataset (results1) only contains the
> result from criteria 1 (twice!) but not the result from criteria 2. I guess
> probably I put the break function under the wrong loop but I cannot fix it.
> Sorry if the whole script looks quite messy, I will be very appreciate if
> someone can help me fix this problem or probably give me some advices to
> write it in a smart way.
>
>
> My attempt:
>
> controlall <- rbeta(10000,1.5,6)
> caseall <- rbeta(10000,1.6,6)
>
> results <- NULL
> results1 <- NULL
> or <-vector("list",length=10)
> criteria <- matrix(data=c(1.05,1.15,1.15,1.25),ncol=2,nrow=2)
>
> for (k in 1:2) {
>
> for (i in 1:1000)
> {
> control <- sample(controlall,100)
> case <- sample(caseall,100)
> or[i] <-
> round(mean(case)*(1-mean(control))/(mean(control)*(1-mean(case))),digit=2)
>
> if (or[i]<criteria[k,2]&or[i]>=criteria[k,1])
> {
> group <- c(rep(1,100),rep(0,100))
> value <- c(case,control)
> ttest <- (t.test(case,control,alternative="two.sided",paired=F))$p.value
> all <- c(or[i],ttest)
> results <- rbind(results,all)
> }
> ifelse(nrow(results)==5,break,1)
> }
> results1 <- rbind(results1,results)
> }
>
>
> Thank you so much in advance,
> Amber
>
> [[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.
--
Jim Holtman
Data Munger Guru
What is the problem that you are trying to solve?
Tell me what you want to do, not how you want to do it.
______________________________________________
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.
