R：所有变量之间的相互作用项测试

Question

我正在尝试对包含23个变量的数据集执行交互测试。以下是所有变量之间的交互项测试循环，但代码之前未保留变量测试。如何进行校正，允许汇总分别显示所有变量（如x y z）之间的交互，而不是x y和y x？

vars=colnames(happiness_logper1)[-1] 
varHap <- rep(0,100)
for (i in vars)  {
for (j in vars) {
if (i != j) {
varHap= paste(i,j,sep='*')}
lm.fit <- lm(paste("Happiness.in.life ~", varHap), data=happiness_logper1)
     print(summary(lm.fit))
}
  }