[問題] 關於逐步迴歸

[ericchin0404 (凜魂)]

[問題類型]:

程式諮詢


[軟體熟悉度]:
請把以下不需要的部份刪除
新手

[問題敘述]:
就是我現在有一個反應變數向量Y
跟一個解釋變數矩陣X，每一行代表一個變數
然後我對他做迴歸
m=lm(Y~X)
然後我想用stepwise選變數，所以寫了
step(m)
不過做出來的結果就是原來的模型，他根本沒有挑選
但如果我把模型寫成
m=lm(Y~X[,1]+X[,2]+X[,3]+X[,4])
他就可以選

所以我的問題是，如果我要做逐步迴歸，可是迴歸模型想用矩陣表示，那程式應該怎麼寫
？
另外想問若想把模型改成有二次項跟交互作用項，要怎麼寫？
麻煩各位，謝謝



--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 42.77.230.44
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1505398422.A.CE2.html
→ andrew43: 你的第一個step(m)只是試試看把維一的X丟掉看看AIC變化
→ andrew43: 而發現AIC沒有變更小所以就維持不動。
→ andrew43: 不嫌棄的話請見拙作 https://youtu.be/en_GdbpZexE
→ andrew43: 等等，你的X是矩陣的話，那其實是把整個X都丟掉看看。
→ andrew43: 你可以as.data.frame(X)後再丟給lm()。
→ andrew43: 如果x有很多欄，就寫成lm(y ~ ., as.data.frame(x))
→ andrew43: 其中的 "." 就是指 as.data.frame(x) 的每個欄位
下面這樣寫可以了 感謝！
那我要讓模型加入二次項跟交互作用項呢？
一樣要能選變數的
比較簡單的寫法是？
※ 編輯: ericchin0404 (42.77.230.44), 09/14/2017 22:33:28
→ andrew43: lm(y ~ (.)^2, as.data.frame(x))
這樣好像只會有一次項跟交互作用項，但不會有平方項，請問還有解決辦法嗎？我試過這
樣寫
lm(y ~ (.)^2+I(X^2), as.data.frame(x))
但同樣沒辦法做逐步挑選
※ 編輯: ericchin0404 (42.77.230.44), 09/15/2017 02:07:32
→ andrew43: 那可以改變做法，先把所有自變項製作好並塞在新的matrix
→ andrew43: 或是塞在新的data frame
→ andrew43: 例如 model.matrix(~.^2, data = as.data.frame(x))[,1]
→ andrew43: 和x^2這二個matrix用cbind()接起來，看習慣。
→ andrew43: 另一種想法可以是利用字串處理生成很長的formula，直接
→ andrew43: 餵給lm()和step()，而不是預先製作變數。
→ ericchin0404: Ok 謝謝