如何从 t 测试函数中提取 R 置信度的下限和上限？

Question

我使用以下代码检索数据的置信度：

out <- t.test(my_data$my_col, conf.level = 0.95)
out

这 returns 类似于：

data:  my_data$my_column
t = 30, df = 20, p-value < 2.1e-14
alternative hypothesis: true mean is not equal to 0
95 percent confidence interval:
 62.23191 80.11201
sample estimates:
mean of x 
 75.10457

我试过：

out[4][1]

但是这个returns:

$conf.int
[1]  62.23191 80.11201
attr(,"conf.level")
[1] 0.95

如何分别从中获取下限和上限？（即如何提取 62.23191 和 80.11201 作为变量？）

Answer 1

t.test() 的输出是一个列表。置信区间存储为 $conf.int 列表元素中的向量。

要访问各个置信区间，请使用 out$conf.int[1] & out$conf.int[2]

示例：

out <- t.test(1:10, y=c(7:20)) 
out$conf.int
#[1] -11.052802  -4.947198
#attr(,"conf.level")
#[1] 0.95
out$conf.int[1]
#[1] -11.0528
out$conf.int[2]
#[1] -4.947198

Answer 2

您可以使用 broom 包，不确定它是否比公认的答案更有效，但它是另一种选择。

library(broom)

tidy(t.test(1:10, y = c(7:20),conf.int = TRUE)

 estimate estimate1 estimate2 statistic   p.value parameter conf.low conf.high method                  alternative
     <dbl>     <dbl>     <dbl>     <dbl>     <dbl>     <dbl>    <dbl>     <dbl> <chr>                   <chr>      
1       -8       5.5      13.5     -5.43 0.0000186      22.0    -11.1     -4.95 Welch Two Sample t-test two.sided

如何从 t 测试函数中提取 R 置信度的下限和上限？

How to extract lower and upper bounds from confidence level in R from t test function?

r

statistics

confidence-interval

t-test