我试图自动化我的一段代码,这样编程就不那么乏味了.
基本上,我试着用rms包中的fastbw()
个逐步 Select 变量.我想将fastbw()
Select 的变量列表传递到一个公式中,即y ~ x1+x2+x3
,"x1""x2""x3"是fastbw()
Select 的变量列表
这是我try 过但没有成功的代码
olsOAW0.r060 <- ols(roll_pct~byoy+trans_YoY+change18m,
subset= helper=="POPNOAW0_r060",
na.action = na.exclude,
data = modelready)
OAW0 <- fastbw(olsOAW0.r060, rule="p", type="residual", sls= 0.05)
vec <- as.vector(OAW0$names.kept, mode="any")
b <- paste(vec, sep ="+") ##I even tried b <- paste(OAW0$names.kept, sep="+")
bestp.OAW0.r060 <- lm(roll_pct ~ b ,
data = modelready,
subset = helper =="POPNOAW0_r060",
na.action = na.exclude)
我是R新手,还没有跟上陡峭的学习曲线,所以为明显的编程错误道歉.