我有一个向量:
vector <- c(12, 17, 24, 35, 23, 34, 56)
如何计算 R 中该向量的置信区间(90%、99%、95%)?
这是我想要的结果示例:在此输入图像描述 https://i.stack.imgur.com/wau32.jpg
这是一个根据 t 分布计算置信区间的函数:
confidence_interval <- function(vector, interval) {
# Standard deviation of sample
vec_sd <- sd(vector)
# Sample size
n <- length(vector)
# Mean of sample
vec_mean <- mean(vector)
# Error according to t distribution
error <- qt((interval + 1)/2, df = n - 1) * vec_sd / sqrt(n)
# Confidence interval as a vector
result <- c("lower" = vec_mean - error, "upper" = vec_mean + error)
return(result)
}
以及您提供的向量和间隔的示例用法:
> vector <- c(12, 17, 24, 35, 23, 34, 56)
> confidence_interval(vector, 0.90)
lower upper
17.97255 39.45602
> confidence_interval(vector, 0.95)
lower upper
15.18797 42.24060
> confidence_interval(vector, 0.99)
lower upper
8.219946 49.208626
And 这是教程 http://www.cyclismo.org/tutorial/R/confidence.html我从中开发了这个方法。
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)