Displaying 1 result from an estimated 1 matches for "tp3005183p3005183".
2010 Oct 21
1
All other variables in upper scope arg for stepAIC
...a<-cbind(y1,y2,y3)
library(MASS)
for (i in 1:3)
{mod<-lm(get(paste("y", i, sep=""))~1)
stepAIC(mod, scope=list(lower=~1, upper=~the_other_y))}
Thanks
Paul
--
View this message in context: http://r.789695.n4.nabble.com/All-other-variables-in-upper-scope-arg-for-stepAIC-tp3005183p3005183.html
Sent from the R help mailing list archive at Nabble.com.