作者horby (horby)
看板Statistics
標題Re: [問題][R]關於參數輸入的格式
時間Mon Apr 7 15:37:45 2008
※ 引述《ablazekandyh (ablazekandyh)》之銘言：
: ※ 引述《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 熟棯的朋友 給予幫助。
: e.g.
: List = c("x1", "x1" "x3", "x1" "x3" "x7")
: for(i in List){
: b = get(i) [這裡好像是 b=get('i')才對]
: update(a, ~ . + b)
: }
你好！ 我一開始也是這樣想的想法，不過不行耶
會出現
Error in model.frame.default(formula = y ~ x1 + b, data = D,
drop.unused.levels = TRUE) : variable lengths differ (found for 'b')
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.109.74.65