假设给我一个数据框df在运行时，如何使用多项式回归拟合多项式模型，每个预测变量都是 df 中的一列，并且常数 k >= 2 的程度
困难在于，“df”是在运行时读取的，因此在编写脚本时其列的数量和名称是未知的。（但我确实知道响应变量是第一列）所以当我调用lm我不知道怎么写formula.
如果 k = 1，那么我可以简单地写出一个通用的线性公式
names(df)[1] <- "y"
lm(y ~ ., data = df)
我可以对多项式公式做类似的事情吗？
一种相当复杂的方法是创建一个formula为了lm通过将术语粘贴在一起进行回归调用。
# some data
dat <- data.frame(replicate(10, rnorm(20)))
# Create formula - apply f function to all columns names excluding the first
form <- formula(paste(names(dat)[1], " ~ ",
paste0("poly(", names(dat)[-1], ", 2)", collapse="+")))
# run regression
lm(form , data=dat)
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
