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