[问题] 回归交互作用问题

楼主: 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