如何使用`r`仅将R命令的某些部分嵌入到R markdown文档中?

时间:2018-09-22 19:33:56

标签: r markdown r-markdown

我正在使用以下命令创建假设检验和置信区间:

```{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 

1 个答案:

答案 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)值。