我有一个包含 16 个变量的数据框。当我进行多元线性回归时，我会执行以下操作：
fit <- lm(y ~ .,data=data)
现在，我知道如何添加其中一个变量的二阶项：
fit2 <- lm(y ~ poly(x1,2) + .,data=data)
但现在我不想为所有 16 个变量都写出这个。如何以简单的方式对所有变量执行此操作？
当假设数据中的第一个变量是我们的“y”时，我们得到：
as.formula(
paste('y ~',paste('poly(',colnames(data[-1]),',2)',collapse = ' + ')
)
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
