[問題] 迴歸交互作用問題

看板R_Language作者 (BLUE)時間6年前 (2017/11/28 18:01), 編輯推噓1(1推 0噓 36→)
留言37則, 5人參與, 6年前最新討論串1/1
[問題類型]: 我想用R 做某件事情，但是我不知道要怎麼用R 寫出來 [軟體熟悉度]: 沒寫過程式，R 是我的第一次 [問題敘述]: 我想跑回歸的交互作用，跑出來會有2個變數的個別係數和交乘項的係數，請問有什麼方法可以只跑交乘項的係數，不要跑個別的， 因為這樣結果會被影響，謝謝! [程式範例]: lm(y~a*b,data=X) 這個會跑出a,b個別係數 ----- Sent from JPTT on my Asus ASUS_Z012DA. -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 42.76.200.42 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1511863306.A.78D.html
→
11/28 18:06, 6年前 , 1F
a和b是因子或是數字？
11/28 18:06, 1F
→
11/28 18:06, 6年前 , 2F
數字的話，a*b改成I(a*b)
11/28 18:06, 2F
→
11/28 18:07, 6年前 , 3F
因子的話，a*b改成interaction(X$a, X$b)
11/28 18:07, 3F
→
11/28 18:09, 6年前 , 4F
I(a*b)就是指數學上的乘法。interaction()是攤平因子
11/28 18:09, 4F
→
11/28 18:17, 6年前 , 5F
那請問一個因子一個數值呢?
11/28 18:17, 5F
→
11/28 19:44, 6年前 , 6F
好奇你這麼做的目的是什麼。通常交互作用不能少了主效果
11/28 19:44, 6F
推
11/28 21:59, 6年前 , 7F
lm(y ~ a:b, data = X)
11/28 21:59, 7F
→
11/28 22:31, 6年前 , 8F
用a:b的話 還是出現a b係數
11/28 22:31, 8F
→
11/28 22:33, 6年前 , 9F
因為有交乘項和沒交乘項 a b係數有時候會差很大 甚至
11/28 22:33, 9F
→
11/28 22:33, 6年前 , 10F
正負變了 所以想說有沒有辦法只跑交乘項
11/28 22:33, 10F
→
11/28 22:35, 6年前 , 11F
我好奇為什麼你還會出現a b係數 @@
11/28 22:35, 11F
→
11/28 22:36, 6年前 , 12F
lm(Sepal.Width ~ Sepal.Length:Species, data = iris)
11/28 22:36, 12F
→
11/28 22:36, 6年前 , 13F
技術上是可行的才對
11/28 22:36, 13F
→
11/28 22:37, 6年前 , 14F
上面的結果裡也沒有main effects
11/28 22:37, 14F
→
11/28 22:39, 6年前 , 15F
統計上你則是要想清楚為什麼你不需要main effects
11/28 22:39, 15F
→
11/28 22:40, 6年前 , 16F
如果有無交乘項會影響到a和b的係數，那應該是從解釋交互
11/28 22:40, 16F
→
11/28 22:40, 6年前 , 17F
作用下手？
11/28 22:40, 17F
→
11/28 23:31, 6年前 , 18F
抱歉 我搞錯了 確實沒有a b係數
11/28 23:31, 18F
→
11/28 23:31, 6年前 , 19F
謝謝哦！
11/28 23:31, 19F
→
11/29 01:04, 6年前 , 20F
建議原po還是先了解「只有交乘項」的模型實際是什麼意義
11/29 01:04, 20F
→
11/29 02:38, 6年前 , 21F
沒有主效果交成像根本沒有意義啊
11/29 02:38, 21F
→
11/29 02:39, 6年前 , 22F
你乾脆先乘完創造一個新變數放進回歸不就好了
11/29 02:39, 22F
→
11/29 10:39, 6年前 , 23F
可是我有一個是因子 可以直接乘嗎
11/29 10:39, 23F
→
11/29 11:36, 6年前 , 24F
直接乘要先把因子轉成數字（例如dummy），但如果不明白
11/29 11:36, 24F
→
11/29 11:37, 6年前 , 25F
它實際的意義要小心解讀。好啦我直接說，你這樣做沒意義
11/29 11:37, 25F
→
11/29 11:39, 6年前 , 26F
（我指一因子一共變數但只考慮相乘項的情況）
11/29 11:39, 26F
→
11/29 11:43, 6年前 , 27F
至於相乘的過程可以自己做，但在R中其實指定交互作用後
11/29 11:43, 27F
→
11/29 11:44, 6年前 , 28F
R就會自動幫你乘好了。
11/29 11:44, 28F
→
11/29 11:44, 6年前 , 29F
如果你是要用來解釋假說，那麼沒主效果比較難解釋，如果是
11/29 11:44, 29F
→
11/29 11:44, 6年前 , 30F
要用來預測而已還好
11/29 11:44, 30F
→
11/29 11:45, 6年前 , 31F
另外，從說明中看來交互作用可能很重要。你可以先畫個
11/29 11:45, 31F
→
11/29 11:45, 6年前 , 32F
交互作用圖看看，說不定就可以解釋你看到的「矛盾」原因
11/29 11:45, 32F
→
11/29 11:55, 6年前 , 33F
單只有預測也要小心一些情況，例如是不平衡設計。
11/29 11:55, 33F
→
11/29 21:28, 6年前 , 34F
謝謝大家~
11/29 21:28, 34F
→
11/29 21:59, 6年前 , 35F
再請問一下，我用a:b方式，但跑出來有NA是為什麼呢?http
11/29 21:59, 35F
→
11/29 21:59, 6年前 , 36F
://i.imgur.com/B8ArFFT.jpg
11/29 21:59, 36F
→
11/29 22:01, 6年前 , 37F
11/29 22:01, 37F
‣ 更多 tony1331 的文章
文章代碼(AID): #1Q7JGAUD (R_Language)