search for: tp18342119p18342119

...group<-factor(rep(1:k,n))) eval.simu<-function(e) { y<-mu[group]+e fit<-aov(y~group) s.fit<-summary(fit)[[1]] s.fit$P[1] } p<-apply(ee,1,"eval.simu") #then count number of p?alpha to get power -- View this message in context: http://www.nabble.com/simulate-data-for-lme-tp18342119p18342119.html Sent from the R help mailing list archive at Nabble.com.