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...