r - stepAIC handling of multinom models -
i seeing weird behavior stepaic
function in mass
package when dealing multinomial logistic models. here sample code:
library(nnet) library(mass) example("birthwt") race.model <- multinom(race ~ smoke, bwt) race.model2 <- stepaic(race.model, k = 2)
in case race.model
, race.model2
have identical terms; stepaic
did not prune anything. however, need query attributes of models, , error race.model2
:
formula(race.model)[2]
returns race()
but
formula(race.model2)[2]
gives error:
error in terms.formula(newformula, specials = names(attr(termobj, "specials"))) : invalid model formula in extractvars
this behavior seems occur when stepaic
not remove terms model. in following code, terms removed stepaic
, , both models can queried:
race.big <- multinom(race ~ ., bwt) race.big2 <- stepaic(race.big, k = 2) formula(race.big)[2] formula(race.big2)[2]
any ideas going wrong here?
Comments
Post a Comment