[R] Stepwise

arun smartpink111 at yahoo.com
Fri Nov 16 00:44:16 CET 2012


Hi,
May be this helps:
set.seed(232)
mat1<-matrix(sample(1:100,80,replace=TRUE),ncol=8) #with 8 columns
dat1<-data.frame(mat1)
 names(dat1)[1]<-"Y"
 form1<-as.formula(paste("Y~",paste(names(dat1)[2:8],collapse="+"))) #should change "8" to "400"

fit.final<-lm(form1,data=dat1)

fit.final
#Call:
#lm(formula = form1, data = dat1)
#
#Coefficients:
#(Intercept)           X2           X3           X4           X5           X6  
# 143.763081    -0.009030    -1.844810     1.436647    -0.005407    -0.050101  
#         X7           X8  
#   0.563285    -1.901918  
A.K.






________________________________
From: farnoosh sheikhi <farnoosh_81 at yahoo.com>
To: arun <smartpink111 at yahoo.com> 
Sent: Thursday, November 15, 2012 6:14 PM
Subject: Stepwise


Hi there,

I want to apply the following code fo my data with 400 predictors.
I was wondering if there is an alternative way instead of typing 400 predictors for the following code.
I really appreciate your help.

fit0<-lm(Y~1, data= mydata)
fit.final<- lm(Y~X1+X2+X3+.....+X400, data=mydata) ???
step(fit0, scope=list(lower=fit0, upper=fit.final), data=mydata, direction="forward")
step(fit.final, scope=list(lower=fit.final, upper=fit0), data=mydata, direction="backward")

Best,Farnoosh Sheikhi 



More information about the R-help mailing list