Re: [問題][R]關於參數輸入的格式
※ 引述《horby (horby)》之銘言：
: 我最主要要問的問題是，當我在原有的回歸模型中要加入一個新的解釋變數
: 比如：D是一個data.frame， 其中 names(D)= "y" "x1" "x2" "x3" ... "x20"
: a=lm(y~x1,data=D)
: 則我們要加入一個新的解釋變數 x5 的話，以 update 來說可以寫成
: update(a, ~ . + x5)
: 但是如果我想改成下面這樣的輸入格式，請問一下大家要怎麼處理
: update(a, ~ . + b)
: 其中 b="x5"
: 註：
: 這是一個關於變數選擇的程式，有點像是 forward selection。將對 y 解釋能力最高的
: x 依序加入回歸模型之中。 在一些運算之後，找出某一個解釋變數（在上例中，就是
: column name 中那個叫做 "x5" 的變數）並將他的名字暫存在 b 中。 而 b 會隨著
: 新加入的解釋變數，而變化。 比如 "x1", "x1" "x3", "x1" "x3" "x7" ....
: 請對 s-plus 或是 R 熟棯的朋友 給予幫助。
可以直接用
D[,which(names(D)=="x5")]
假設你每次找出來最好的變數暫存在 temp2 裡
temp1 <- numeric(20)
temp2 <- [x1,x3,x7....]
for (i in 1:20){
temp1[i] <- which(names(D)==temp2[i])
update(a, ~ . + D[,temp1[1:i]])
}
我想你還是得自己改一下以符合你的需要
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 68.42.72.37
推
04/08 03:01, , 1F
04/08 03:01, 1F
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 4 之 4 篇)：