R语言如何选择显著的二阶交互作用？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我现在有一个400*2000的矩阵，即有400个样本，2000个变量，我现在想要求这2000个变量两两之间的交互作用，然后选择出有交互作用显著的项，请问用R语言该怎么操作？？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：交互作用 R语言 怎么操作 R语言 交互作用

自己写for loop把二阶项算出来，然后用glmnet这个包去fit一个LASSO/Ridge/elasticNet，这样就能知道经过regularize之后的model～

foozhencheng 发表于 2017-5-28 03:10
自己写for loop把二阶项算出来，然后用glmnet这个包去fit一个LASSO/Ridge/elasticNet，这样就能知道经过reg ...

要用到for循环吗?我之前试过直接用x=model.matrix(y~.^2,data)来产生有主效应和二阶交互的矩阵，然后考虑用lasso来选择变量，但是由于x矩阵的维数很大，电脑根本跑不出来

Heavenlr 发表于 2017-6-2 17:43
要用到for循环吗?我之前试过直接用x=model.matrix(y~.^2,data)来产生有主效应和二阶交互的矩阵，然后考虑 ...

这个直接平方design matrix啊，那可能是有些困难。不过你变量数有多少？

Heavenlr 发表于 2017-6-2 17:43
要用到for循环吗?我之前试过直接用x=model.matrix(y~.^2,data)来产生有主效应和二阶交互的矩阵，然后考虑 ...

哦，看到了。2000个变量、400个样本，直接跑LASSO啊，别再弄interaction了，那样的话，真的跑不动～