Skip to content
Prev 87281 / 398502 Next

help with step()

Smith, Phil wrote:
Could be your "data" argument. ?all is a base function. Try re-naming to 
something else. Here's an example:

set.seed(1)
z <- all <- data.frame(x = rnorm(10), y = rnorm(10))

glm.all <- glm(y ~ 1, data = all)
step(glm.all, list(upper = y ~ x, lower = y ~ 1))

# Start:  AIC= 32.67
# y ~ 1
#
# Error in eval(expr, envir, enclos) : invalid 'envir' argument

glm.z <- glm(y ~ 1, data = z)
step(glm.z, list(upper = y ~ x, lower = y ~ 1))

# Start:  AIC= 32.67
#  y ~ 1
#
#        Df Deviance    AIC
# <none>      10.295 32.669
# + x     1    8.834 33.139
#
# Call:  glm(formula = y ~ 1, data = z)
#
# Coefficients:
# (Intercept)
#      0.2488
#
# Degrees of Freedom: 9 Total (i.e. Null);  9 Residual
# Null Deviance:      10.29
# Residual Deviance: 10.29        AIC: 32.67