对于作业,需要使用数据集的交互项来确定线性模型,我可以通过创建不同的线性模型然后执行ANOVA测试来手动完成。
我想知道的是,如果R中有任何函数可以检查给定主效应模型的线性模型交互的所有可能组合吗? (与处理逐步变量选择的方式类似的函数)
我看过了,但是我对R的了解并不强烈,我可能已经浏览了这个功能。
答案 0 :(得分:0)
第一列是因变量(y)(mainClass),我们制作所有可能的线性回归模型
有三个独立变量(X,Y和Z)
mainClass~X,mainClass~X + Y,mainClass~Y + Z以及所有其他组合被探索。
- 代码说明
创建一个数据框,其中第一列是mainClass,其他列是变量
查找不包括第一列(mainClass)的列名称
具有列号的组合
具有用于线性回归的公式组合
DF<- data.frame(mainClass=1:10,X=rnorm(10,3,1),Y=seq(4,40,4),Z=seq(100,10,-10))
Cols <- names(DF)
Cols <- Cols[! Cols %in% "mainClass"]
comb_id <- unlist(lapply(1:n, function(i)combn(1:n,i,simplify=F)),recursive=F)
公式&lt; - sapply(comb_id,function(i)paste(&#34; mainClass~&#34;,paste(Cols [i],collapse =&#34; +&#34;)))<登记/> lapply(公式,函数(i)lm(as.formula(i),data = DF))