在下面的 R 会话中,summary(model)将公式显示为model_str。如何让它显示为mpg ~ cyl + hp同时仍然能够通过字符串设置模型公式?
> data(mtcars)
> names(mtcars)
[1] "mpg" "cyl" "disp" "hp" "drat" "wt" "qsec" "vs" "am" "gear" "carb"
> model_str <- 'mpg ~ cyl + hp'
> model <- lm(model_str, data=mtcars)
> summary(model)
Call:
lm(formula = model_str, data = mtcars)
Residuals:
Min 1Q Median 3Q Max
-4.4948 -2.4901 -0.1828 1.9777 7.2934
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 36.90833 2.19080 16.847 < 2e-16 ***
cyl -2.26469 0.57589 -3.933 0.00048 ***
hp -0.01912 0.01500 -1.275 0.21253
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 3.173 on 29 degrees of freedom
Multiple R-squared: 0.7407, Adjusted R-squared: 0.7228
F-statistic: 41.42 on 2 and 29 DF, p-value: 3.162e-09
Run Code Online (Sandbox Code Playgroud)
使用do.call以便model_str在发送到之前进行评估lm但引用mtcars以便它不是(否则会有大量输出显示 中的实际值mtcars)。
do.call("lm", list(as.formula(model_str), data = quote(mtcars)))
Run Code Online (Sandbox Code Playgroud)
给予:
Call:
lm(formula = mpg ~ cyl + hp, data = mtcars)
Coefficients:
(Intercept) cyl hp
36.90833 -2.26469 -0.01912
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
108 次 |
| 最近记录: |