search for: currentp

Displaying 1 result from an estimated 1 matches for "currentp".

Did you mean: current
2008 Jun 20
1
Unexpected Behavior (potentially) in t.test
...ep="")) #feedback for (lcv in 1:t1Nsim) #Run simulation t1Nsim times... { if (lcv %% 20000 == 0) {print(paste((lcv/t1Nsim)*100,"%",sep=""))} #feedback Cgroup <- rnorm(t1N) #Initial random draw for Group1 Tgroup <- rnorm(t1N) #Initial random draw for Group 2 currentp <- t.test(Cgroup,Tgroup)[["p.value"]] #Get t1 p value t1pvals[lcv] <- currentp #Store t1 p value #If p >= .05 or <= continue threshold then run more subjects if ((currentp <= contthreshold) & (currentp >= .05)) { Cgroup <- c(Cgroup,rnorm(t2N)) #Add t2N subject...