Question

您好，我想知道执行以下问题的命令是什么

为获胜游戏的平均数量确定95％的置信区间。游泳，在Rstudio中限制2 000

viex<-c(2205,2096, 1847, 1903, 1457, 1848, 1564, 1821, 2577, 2476, 1984, 1917, 1761, 1709, 1901, 2288, 2072, 2861, 2411, 2289, 2203 ,2592, 2053 ,1979, 2048, 1786, 2876, 2560)

clasy<-c(10, 11, 11, 13, 10, 11, 10 ,11,  4 , 2,  7, 10,  9,  9,  6,  5,  5,  5,  6,  4 , 3 , 3 , 4, 10,  6 , 8,  2 , 0)

modelo = lm(clasy~ viex)
modelo

这就是我所做的

Answer 1

这直接来自置信区间的定义。

ci_regression <- function(x, conf = 0.95){
  est <- coef(summary(x))[, 1]
  se <- coef(summary(x))[, 2]
  qq <- qt(1 - (1 - conf)/2, df = x$df.residual)
  cbind(lower = est - qq*se, upper = est + qq*se)
}

ci_regression(modelo)
#                   lower        upper
#(Intercept) 16.246064040 27.330437725
#viex        -0.009614347 -0.004435854

对于线性回归，有一个confint.lm方法用于类"lm"的对象。
只需运行

confint(modelo)
#                   2.5 %       97.5 %
#(Intercept) 16.246064040 27.330437725
#viex        -0.009614347 -0.004435854

平均数量的置信区间

1 个答案: