while I'm struggling with associating labels to lm() terms, I found out this cool thing (well, you know what I mean):
Code: Select all
afac<-replicate(100,sample(2,1))
afac<-factor(afac)
afac2<-rnorm(100,0,1)
y<-rnorm(100,0,1)
model<-lm(y~afac+afac2)
### Regression ####
summary(model)
### ANOVA ####
drop1(model,.~.,test="F")