我正在寻找B1(newx)线性模型系数的值,而不是名称。我只想要0.5值。我不想要名字“newx”。
newx <- c(0.5,1.5.2.5)
newy <- c(2,3,4)
out <- lm(newy ~ newx)
out
看起来像:
Call:
lm(formula = newy ~ newx)
Coefficients:
(Intercept) newx
1.5 1.0
我来到这里。但现在我被卡住了。
out$coefficients["newx"]
newx
1.0
答案 0 :(得分:91)
对于像这样的单个元素,请使用[[
而不是[
。比较:
coefficients(out)["newx"]
# newx
# 1
coefficients(out)[["newx"]]
# [1] 1
更一般地说,使用unname()
:
unname(coefficients(out)[c("newx", "(Intercept)")])
# [1] 1.0 1.5
head(unname(mtcars))
# NA NA NA NA NA NA NA NA NA NA NA
# Mazda RX4 21.0 6 160 110 3.90 2.620 16.46 0 1 4 4
# Mazda RX4 Wag 21.0 6 160 110 3.90 2.875 17.02 0 1 4 4
# Datsun 710 22.8 4 108 93 3.85 2.320 18.61 1 1 4 1
# Hornet 4 Drive 21.4 6 258 110 3.08 3.215 19.44 1 0 3 1
# Hornet Sportabout 18.7 8 360 175 3.15 3.440 17.02 0 0 3 2
# Valiant 18.1 6 225 105 2.76 3.460 20.22 1 0 3 1
## etc.
答案 1 :(得分:4)
如果有关删除姓名的问题,这是另一种方法
my_vec <- # some quantile function(returns named vector)
names(my_vec) <- NULL
my_vec
## [1] 1 2 3
答案 2 :(得分:0)
另一种方法是使用broom
软件包:
broom::tidy(out)$estimate[1]
#1.5
答案 3 :(得分:0)
一个简单而直接的方法是
as.numeric(out$coefficients["newx"])