search for: vardep

Hello, as far as I see, R reports type I sums of squares. I'd like to get R to print out type III sums of squares. e.g. I have the following model: vardep~factor1*factor2 to get the type III sum of squares for factor1 I've tried anova(lm(vardep~factor2+factor1:factor2),lm(vardep~factor1*factor2)) but that didn't yield the desired result. Could anyone give me a hint how to proceed? thanks in advance Josef

...se it in spite of hard heteroscedasticity ? 4/ A last question concerning SAS. My model appears to not converge in SAS, indicating a "structure" in the variance. Is it implying something in lmer or lme ? Many Thanks Here is the function for the permutation : permut<-function(data,vardep,var1,var2,var.random,nboot=10){ # initialisation ms.var1<-numeric(length=nboot) ms.var2<-numeric(length=nboot) ms.var3<-numeric(length=nboot) var1.permut<-numeric(length=length(var1)) var2.permut<-numeric(length=length(var2)) var3.permut<-numeric(length=length(var2)) var1<-fac...