Hmisc的'mean_cl_normal函数'是如何工作的？

Question

mean_cl_normal function的内部结构是什么，可以安全地用于置信区间计算？它是如何工作的？

Answer 1

mean_cl_normal是ggplot的一部分。它在Hmisc中使用smean.cl.normal(...)，它使用t分布的分位数来计算置信限。

set.seed(1)   # for reproducibility
x <- rnorm(100, mean=5)
library(Hmisc)
smean.cl.normal(x)
#     Mean    Lower    Upper 
# 5.108887 4.930665 5.287110 
mean(x) + c(0,1,-1) * qt(0.05/2,(length(x)-1))*sd(x)/sqrt(length(x))
# [1] 5.108887 4.930665 5.287110