我試圖使我的代碼片段自動(dòng)化，以使編程變得不再那么乏味。基本上，我試圖fastbw()在rms包中使用逐步選擇變量。我想將選擇的變量列表傳遞fastbw()到公式中y ~ x1+x2+x3，因?yàn)椤?x1”“ x2”“ x3”是選擇的變量列表fastbw()這是我嘗試但無法正常工作的代碼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語言的新手，仍然沒有走過陡峭的學(xué)習(xí)曲線，因此對(duì)明顯的編程失誤表示歉意。
3 回答
侃侃無極
TA貢獻(xiàn)2051條經(jīng)驗(yàn) 獲得超10個(gè)贊
你快到了。你只需要paste在整個(gè)配方一起，像這樣：
paste("roll_pct ~ ",b,sep = "")
使用將其強(qiáng)制轉(zhuǎn)換為實(shí)際公式as.formula，然后將其傳遞給lm。從技術(shù)上講，我認(rèn)為lm可以強(qiáng)制字符串本身，但是通常自己更強(qiáng)制。（某些期望公式的函數(shù)不會(huì)為您提供強(qiáng)制性，其他函數(shù)則可以。）
ITMISS
TA貢獻(xiàn)1871條經(jīng)驗(yàn) 獲得超8個(gè)贊
我今天遇到了類似的問題，如果您想使其更加通用，甚至不需要固定的類名，則可以使用
frmla <- as.formula(paste(colnames(modelready)[1], paste(colnames(modelready)[2:ncol(modelready)], sep = "",
collapse = " + "), sep = " ~ "))
假定您在第一列中具有類變量或因變量，但索引可以很容易地切換到最后一列，如下所示：
frmla <- as.formula(paste(colnames(modelready)[ncol(modelready)], paste(colnames(modelready)[1:(ncol(modelready)-1)], sep = "",
collapse = " + "), sep = " ~ "))
然后繼續(xù)lm使用：
bestp.OAW0.r060 <- lm(frmla , data = modelready, ... )
- 3 回答
- 0 關(guān)注
- 817 瀏覽
添加回答
舉報(bào)
0/150
提交
取消