Displaying 1 result from an estimated 1 matches for "road_d".
Did you mean:
road_de
2009 Jan 26
1
glm StepAIC with all interactions and update to remove a term vs. glm specifying all but a few terms and stepAIC
...scope argument is
scope=list(upper=~.^2,lower=NULL). This means there are 78 predictor terms
considered, the 12 primary terms and 66 interactions [n(n+1)/2]. I see this
with trace=T also. Here is the code used:
>glm1<-glm(formula = PRESENCE == "1" ~ SNOW + I(SNOW^2) + POP_DEN + ROAD_DE
+ ADJELEV + I(ADJELEV^2) + TRI + I(TRI^2) + EDGE + I(EDGE^2) + TREECOV +
I(TREECOV^2),family = binomial, data = wolv)
summary(glm1)
>library(MASS)
>stepglm2<-stepAIC(glm1,scope=list(upper=~.^2,lower=NULL),
trace=T,k=log(4828),direction="both")
> summary(stepglm2)...