从 lm 对象中提取因变量
Extracting dependent variable from lm object
是否有从 lm 对象中提取 Y 的函数?
我使用 residual(m) 和 predict(m) 但我使用对象内部结构来提取 Y...
m = lm(Y ~ X1, d)
head(m$model$Y)
[1] -0.791214 -1.291986 -0.472839 1.940940 -0.977910 -1.705539
您可以使用 model.frame()
,如下所示:
# From the stats::lm documentation
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
lm1 <- lm(weight ~ group)
model.frame(lm1)$weight
## [1] 4.17 5.58 5.18 6.11 4.50 4.61 5.17 4.53 5.33 5.14 4.81 4.17 4.41
## 3.59 5.87 3.83 6.03 4.89 4.32 4.69
如果您对公式中的一个或多个变量调用函数,例如
lm2 <- lm(log(weight) ~ group)
您可以使用 get_all_vars(lm2)$weight
获得未转换的值(model.frame()
returns 转换后的值)。
如果您想查看特定 class 可用的函数(特别是提取函数),您可以使用 methods(class = "lm")
(或您感兴趣的任何对象 class)进行检查中).
是否有从 lm 对象中提取 Y 的函数?
我使用 residual(m) 和 predict(m) 但我使用对象内部结构来提取 Y...
m = lm(Y ~ X1, d)
head(m$model$Y)
[1] -0.791214 -1.291986 -0.472839 1.940940 -0.977910 -1.705539
您可以使用 model.frame()
,如下所示:
# From the stats::lm documentation
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
lm1 <- lm(weight ~ group)
model.frame(lm1)$weight
## [1] 4.17 5.58 5.18 6.11 4.50 4.61 5.17 4.53 5.33 5.14 4.81 4.17 4.41
## 3.59 5.87 3.83 6.03 4.89 4.32 4.69
如果您对公式中的一个或多个变量调用函数,例如
lm2 <- lm(log(weight) ~ group)
您可以使用 get_all_vars(lm2)$weight
获得未转换的值(model.frame()
returns 转换后的值)。
如果您想查看特定 class 可用的函数(特别是提取函数),您可以使用 methods(class = "lm")
(或您感兴趣的任何对象 class)进行检查中).