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

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

Links booklink

Contact Us: admin [ a t ] ucptt.com