作者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/cn.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