控制与线性模型函数相关的置信区间的打印输出

Posted 2023-02-14

【中文标题】控制与线性模型函数相关的置信区间的打印输出

【英文标题】：Control the printout of confidence intervals related to a linear model function

【发布时间】：2022-01-17 03:14:49

【问题描述】：

我正在为线性模型进行引导，但如何编辑截距和 x 变量的打印输出名称？

这是模拟数据

set.seed(42) 
n <- 100
x <- rnorm(n)
e <- rnorm(n)
y <- as.numeric(50 + 25*x + e)
dd <- data.frame(id=1:n, x=x, y=y)

这是模型：

mo <- lm(y ~ x, data=dd)

找到拟合和残差：

fit <- fitted(mo)
resi <- residuals(mo)

基于残差自举检索置信区间的函数：

FUN <- function() 
  X <- model.matrix(mo)
  ressampy <- fit + sample(resi, length(resi), replace = TRUE)
  bootmod <- lm(ressampy ~ X-1)
  confint(bootmod, level = 0.95)

1 次运行的输出（注意打印输出是 X(Intercept) 和 Xx，但我只希望它们是 (Intercept) 和 x）

FUN()
                2.5 %   97.5 %
X(Intercept) 49.74439 50.07817
Xx           24.92904 25.25103

这可能是一个简单的解决方法，但我无法让它工作。任何帮助将不胜感激！

【参考方案1】：

只需使用rownames() 更改包含置信区间的矩阵的行名，如下所示：

FUN <- function() 
  X <- model.matrix(mo)
  ressampy <- fit + sample(resi, length(resi), replace = TRUE)
  bootmod <- lm(ressampy ~ X-1)
  ci <- confint(bootmod, level = 0.95)
  rownames(ci) <- c("(Intercept)", "x")
  return(ci)

【讨论】：

哇，就这么简单。谢谢！