我正在使用以下命令创建假设检验和置信区间:
```{r}
t.test(urtak1$fermetraverd ~ urtak1$matssvaedi)
```
我需要将t值和p值以及置信区间(分别)嵌入到我的R降价文档的文本中。如何使用“ r”执行此操作?
这是我的输出:
Welch Two Sample t-test
data: urtak1$fermetraverd by urtak1$matssvaedi
t = 1.0812, df = 96.784, p-value = 0.2823
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-12.02648 40.80856
sample estimates:
mean in group Hagar mean in group Kringlan
348.5697 334.1787
答案 0 :(得分:1)
您可以使用broom
软件包进行此操作。我没有您正在使用的问题数据,因此我使用iris
进行了快速设置。
在第一块中,您将进行分析并创建一个对象,该对象保存t检验输出的“整齐”版本。例如
```{r}
library(dplyr)
library(broom)
example <- iris %>%
filter(Species != "setosa") %>%
droplevels()
result <- t.test(example$Sepal.Length ~ example$Species)
tidy_result <- tidy(result)
```
看看tidy_result
:扫帚软件包的神奇之处在于,它可以在整洁的数据框中提取t检验(以及许多其他检验)的所有输出,您可以在其他地方使用。
现在您可以在文本中使用它,对r tidy_result$statistic
进行评估以显示t统计量,对r tidy_result$p.value
进行评估以显示p值(反引号未正确显示,但应放置在放在tidy_result对象的'r'表达式前面。
在您的Rmarkdown中,它可能看起来像这样:
the t-test resulted in a t-statistic of `r tidy_result$statistic`
with a p-value of `r tidy_result$p.value`
要查看可以包含的所有参数,请在控制台中评估tidy_result
。您将看到:
> tidy_result
estimate estimate1 estimate2 statistic p.value parameter conf.low conf.high method alternative
1 -0.652 5.936 6.588 -5.629165 1.866144e-07 94.02549 -0.8819731 -0.4220269 Welch Two Sample t-test two.sided
因此您可以选择以下任意一项:
"estimate" "estimate1" "estimate2"
"statistic" "p.value" "parameter"
"conf.low" "conf.high" "method" "alternative"
请注意,为了获得置信区间,您可以使用低(conf.low
)和高(conf.high
)值。