3 回答
TA贡献2051条经验 获得超10个赞
你快到了。你只需要paste在整个配方一起,像这样:
paste("roll_pct ~ ",b,sep = "")
使用将其强制转换为实际公式as.formula,然后将其传递给lm。从技术上讲,我认为lm可以强制字符串本身,但是通常自己更强制。(某些期望公式的函数不会为您提供强制性,其他函数则可以。)
TA贡献1871条经验 获得超8个赞
我今天遇到了类似的问题,如果您想使其更加通用,甚至不需要固定的类名,则可以使用
frmla <- as.formula(paste(colnames(modelready)[1], paste(colnames(modelready)[2:ncol(modelready)], sep = "",
collapse = " + "), sep = " ~ "))
假定您在第一列中具有类变量或因变量,但索引可以很容易地切换到最后一列,如下所示:
frmla <- as.formula(paste(colnames(modelready)[ncol(modelready)], paste(colnames(modelready)[1:(ncol(modelready)-1)], sep = "",
collapse = " + "), sep = " ~ "))
然后继续lm使用:
bestp.OAW0.r060 <- lm(frmla , data = modelready, ... )
- 3 回答
- 0 关注
- 769 浏览
添加回答
举报