[問題類型]:
程式諮詢
[軟體熟悉度]:
請把以下不需要的部份刪除
新手
[問題敘述]:
就是我現在有一個反應變數向量Y
跟一個解釋變數矩陣X，每一行代表一個變數
然後我對他做迴歸
m=lm(Y~X)
然後我想用stepwise選變數，所以寫了
step(m)
不過做出來的結果就是原來的模型，他根本沒有挑選
但如果我把模型寫成
m=lm(Y~X[,1]+X[,2]+X[,3]+X[,4])
他就可以選
所以我的問題是，如果我要做逐步迴歸，可是迴歸模型想用矩陣表示，那程式應該怎麼寫
？
另外想問若想把模型改成有二次項跟交互作用項，要怎麼寫？
麻煩各位，謝謝

你的第一個step(m)只是試試看把維一的X丟掉看看AIC變化而發現AIC沒有變更小所以就維持不動。不嫌棄的話請見拙作 https://youtu.be/en_GdbpZexE等等，你的X是矩陣的話，那其實是把整個X都丟掉看看。你可以as.data.frame(X)後再丟給lm()。如果x有很多欄，就寫成lm(y ~ ., as.data.frame(x))其中的 "." 就是指 as.data.frame(x) 的每個欄位lm(y ~ (.)^2, as.data.frame(x))那可以改變做法，先把所有自變項製作好並塞在新的matrix或是塞在新的data frame例如 model.matrix(~.^2, data = as.data.frame(x))[,1]和x^2這二個matrix用cbind()接起來，看習慣。另一種想法可以是利用字串處理生成很長的formula，直接餵給lm()和step()，而不是預先製作變數。

Ok 謝謝