Displaying 1 result from an estimated 1 matches for "frdm".
Did you mean:
fdm
2004 Apr 17
0
nlme - sum of squares - permutation test
...ble F
# the permutation number is kept in the ouput attributes and thus
# can be called from there
an<-anova(lme.model)
Fobs<-an[[3]]
n<-rep(0,length(an[[3]]))
for (i in 1:pn){
lm1<-update(lme.model, sample(.)~.)
an<-anova(lm1)
Frdm<-an[[3]]
for (j in 1:(length(an[[3]]))){
if (Frdm[j]>=Fobs[j]) n[j]<-n[j]+1
}
}
n<-n/pn
names(n)<-row.names(an)
n<-as.data.frame(n)
names(n)<-"p.value"
attributes(n)<-c(attributes(n),permutations=p...