我有一个包含17个变量的数据集
此链接提供了数据 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)
其中炸弹*温度显着
我想做的是针对每个其他变量测试每个变量,
喜欢做
炸弹*温度
炸弹*凝固汽油弹
IG88 *武器
飞弹*火力
等。直到两个人的组合才被挖出
这样,我可以发现每个变量之间是否存在重要的相互作用。
我知道如何手动完成,创建线性模型,然后对该模型进行总结,但我希望能够创建一个测试每个变量的循环,因为测试所有内容的条目很多。 / p>
答案 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))
}}