首先,让我介绍一个示例数据。
set.seed(1)
x1=rnorm(10)
y=as.factor(sample(c(1,0),10,replace=TRUE))
x2=sample(c('Young','Middle','Old'),10,replace=TRUE)
model1 <- glm(y~as.factor(x1>=0)+as.factor(x2),binomial)
当我输入
summary(model1)
时,会得到以下结果。 Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.1835 1.0926 -0.168 0.867
as.factor(x1 >= 0)TRUE 0.7470 1.7287 0.432 0.666
as.factor(x2)Old 0.7470 1.7287 0.432 0.666
as.factor(x2)Young 18.0026 4612.2023 0.004 0.997
请忽略模型估计结果,因为数据是假的。
在R中,有没有一种方法可以更清楚地更改出现在最左列的估计量名称?例如,删除“as.factor”,并在因子水平前加上
_
。输出应如下所示: Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.1835 1.0926 -0.168 0.867
(x1 >= 0)_TRUE 0.7470 1.7287 0.432 0.666
(x2)_Old 0.7470 1.7287 0.432 0.666
(x2)_Young 18.0026 4612.2023 0.004 0.997
summary.(g)lm
对象的打印方法,还是仅仅想要重命名这个特定summary.glm
对象的coefficients
组件的rownames
? - BenBarnescontr.treatment
)或明确设置的对比函数。可以查看car::contr.Treatment
以了解如何更改名称,可以使用该函数或者按照该函数的代码编写自己的函数来更改名称。请注意不要改变原文的意思。 - Brian Diggs