[R] Problem with stepAIC, "invalid second argument"

Klaus Thul klaus_thul at yahoo.com
Sun Dec 26 16:14:42 CET 2004


Dear all,

am I doing something wrong, or is this a bug (R 2.0.1 & 1.9.1,
linux&windows):

    > library(MASS)
    > df <- data.frame(x = rnorm(100), y = rnorm(100), z = rnorm(100))
    > f <- lm(x ~ y + z, data = df)
    > stepAIC(f, scope = list(upper = ~y * z, lower = ~1))
    Start:  AIC= -16.56
     x ~ y + z
    Error in eval(expr, envir, enclos) : invalid second argument
    >


I can avoid the problem doing the following, but in my application I do
not want to use attach / detach because of possible name conflicts:

    > attach(df)
    > f <- lm(x ~ y + z)
    > stepAIC(f, scope = list(upper = ~y * z, lower = ~1))
    Start:  AIC= 23.55
     x ~ y + z

    ...

    Call:
    lm(formula = x ~ 1)

    Coefficients:
    (Intercept)
        -0.1972
    > detach(df)


Thanks for any help,
Klaus Thul




More information about the R-help mailing list