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

3306

收藏 2017-05-27

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

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

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

全部回复

2017-5-28 03:10:09

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

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

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

2017-6-2 17:43:07

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

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

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

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

2017-6-3 02:25:05

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

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

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

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

2017-6-3 02:26:26

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

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

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

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

分享