R语言如何计算一个数据集中每行数据可信区间的下限和上限
时间: 2024-05-28 13:08:51 浏览: 209
可以使用R语言中的t分布来计算每行数据的置信区间。下限和上限可以分别通过均值减去和加上t乘以标准误差得到。具体代码如下:
confidence_interval <- function(x, confidence_level) {
n <- length(x)
mean_x <- mean(x)
std_err <- sd(x) / sqrt(n)
t_value <- qt(1 - (1 - confidence_level) / 2, n - 1)
lower <- mean_x - t_value * std_err
upper <- mean_x + t_value * std_err
return(list(lower = lower, upper = upper))
}
# Example usage
set.seed(123)
data <- data.frame(matrix(rnorm(100, 5, 2), ncol = 10))
results <- apply(data, 1, confidence_interval, confidence_level = 0.95)
# Display results
head(results)
#> [[1]]
#> [1] 3.502921 6.837278
#>
#> [[2]]
#> [1] 3.846519 6.044373
#>
#> [[3]]
#> [1] 3.652838 6.520422
#>
#> [[4]]
#> [1] 2.361098 6.076619
#>
#> [[5]]
#> [1] 3.085592 6.035831
#>
#> [[6]]
#> [1] 3.289702 6.793412
阅读全文