data(Fires)
names(Fires)
#Fit full model and compare with anova and best AIC subset
#Fit full model
bestglm(Fires, IC="BICq", t=1)
#ANOVA for xyarea is significant at 1.1%.
summary(aov(lburned~xyarea, data=Fires))
bestglm(Fires, IC="AIC")
Run the code above in your browser using DataLab