如何使用for循环查找R中几个变量之间的交互?

时间:2015-05-07 01:19:37

标签: r regression linear

我有一个包含17个变量的数据集 enter image description here

此链接提供了数据 http://www.uwyo.edu/crawford/stat3050/final%20project/maxwellchandler.txt

我想找到变量之间的重要交互。

例如

    fitcivilian<-lm(Civilian~Stock+Terrorism+log(Firepower)+Payload+Bombs*Temperature+FirstAid+Spies+Personnel+IG88, data=data)

enter image description here

其中炸弹*温度显着

我想做的是针对每个其他变量测试每个变量,

喜欢做

炸弹*温度

炸弹*凝固汽油弹

IG88 *武器

飞弹*火力

等。直到两个人的组合才被挖出

这样,我可以发现每个变量之间是否存在重要的相互作用。

我知道如何手动完成,创建线性模型,然后对该模型进行总结,但我希望能够创建一个测试每个变量的循环,因为测试所有内容的条目很多。 / p>

1 个答案:

答案 0 :(得分:5)

我做过类似的事情。您需要根据需要修改循环。如果您需要帮助,请告诉我。

    vars=colnames(mydata)[-1] 
    for (i in vars)  {
    for (j in vars) {
    if (i != j) {
    factor= paste(i,j,sep='*')}
    lm.fit <- lm(paste("Sales ~", factor), data=mydata)
    print(summary(lm.fit))
    }}