示例代码:
require(datasets); data(InsectSprays)
model1 <- lm(count ~ spray, data = InsectSprays)
summary(model1)
但是喷雾的级别不是1、2、3、4、5,而是
我想知道如何将spray1/2/3/4/5显示为sprayB/C/D/E/F?我看到其他人可以使用相同的代码获得显示我想要的内容(显示真实水平)参见这里:
您可以尝试一下下面的
factor
吗
InsectSprays$spray <- factor(InsectSprays$spray, labels = c('1','2','3','4','5','6'))
model1 <- lm(count ~ spray, data = InsectSprays)
summary(model1)
Call:
lm(formula = count ~ spray, data = InsectSprays)
Residuals:
Min 1Q Median 3Q Max
-8.333 -1.958 -0.500 1.667 9.333
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 14.5000 1.1322 12.807 < 2e-16 ***
spray2 0.8333 1.6011 0.520 0.604
spray3 -12.4167 1.6011 -7.755 7.27e-11 ***
spray4 -9.5833 1.6011 -5.985 9.82e-08 ***
spray5 -11.0000 1.6011 -6.870 2.75e-09 ***
spray6 2.1667 1.6011 1.353 0.181
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 3.922 on 66 degrees of freedom
Multiple R-squared: 0.7244, Adjusted R-squared: 0.7036
F-statistic: 34.7 on 5 and 66 DF, p-value: < 2.2e-16