看板 R_Language 關於我們 聯絡資訊
[問題類型]: 我想用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
andrew43: a和b是因子或是數字? 11/28 18:06
andrew43: 數字的話,a*b改成I(a*b) 11/28 18:06
andrew43: 因子的話,a*b改成interaction(X$a, X$b) 11/28 18:07
andrew43: I(a*b)就是指數學上的乘法。interaction()是攤平因子 11/28 18:09
tony1331: 那請問一個因子一個數值呢? 11/28 18:17
andrew43: 好奇你這麼做的目的是什麼。通常交互作用不能少了主效果 11/28 19:44
clsmbstu: lm(y ~ a:b, data = X) 11/28 21:59
tony1331: 用a:b的話 還是出現a b係數 11/28 22:31
tony1331: 因為有交乘項和沒交乘項 a b係數有時候會差很大 甚至 11/28 22:33
tony1331: 正負變了 所以想說有沒有辦法只跑交乘項 11/28 22:33
clsmbstu: 我好奇為什麼你還會出現a b係數 @@ 11/28 22:35
clsmbstu: lm(Sepal.Width ~ Sepal.Length:Species, data = iris) 11/28 22:36
clsmbstu: 技術上是可行的才對 11/28 22:36
clsmbstu: 上面的結果裡也沒有main effects 11/28 22:37
clsmbstu: 統計上你則是要想清楚為什麼你不需要main effects 11/28 22:39
clsmbstu: 如果有無交乘項會影響到a和b的係數,那應該是從解釋交互 11/28 22:40
clsmbstu: 作用下手? 11/28 22:40
tony1331: 抱歉 我搞錯了 確實沒有a b係數 11/28 23:31
tony1331: 謝謝哦! 11/28 23:31
andrew43: 建議原po還是先了解「只有交乘項」的模型實際是什麼意義 11/29 01:04
wearytolove: 沒有主效果交成像根本沒有意義啊 11/29 02:38
wearytolove: 你乾脆先乘完創造一個新變數放進回歸不就好了 11/29 02:39
tony1331: 可是我有一個是因子 可以直接乘嗎 11/29 10:39
andrew43: 直接乘要先把因子轉成數字(例如dummy),但如果不明白 11/29 11:36
andrew43: 它實際的意義要小心解讀。好啦我直接說,你這樣做沒意義 11/29 11:37
andrew43: (我指一因子一共變數但只考慮相乘項的情況) 11/29 11:39
andrew43: 至於相乘的過程可以自己做,但在R中其實指定交互作用後 11/29 11:43
andrew43: R就會自動幫你乘好了。 11/29 11:44
jyung: 如果你是要用來解釋假說,那麼沒主效果比較難解釋,如果是 11/29 11:44
jyung: 要用來預測而已還好 11/29 11:44
andrew43: 另外,從說明中看來交互作用可能很重要。你可以先畫個 11/29 11:45
andrew43: 交互作用圖看看,說不定就可以解釋你看到的「矛盾」原因 11/29 11:45
andrew43: 單只有預測也要小心一些情況,例如是不平衡設計。 11/29 11:55
tony1331: 謝謝大家~ 11/29 21:28
tony1331: 再請問一下,我用a:b方式,但跑出來有NA是為什麼呢?http 11/29 21:59
tony1331: ://i.imgur.com/B8ArFFT.jpg 11/29 21:59
tony1331: http://i.imgur.com/F9NkOy0.jpg 11/29 22:01