-
Notifications
You must be signed in to change notification settings - Fork 7
/
F-rsimsumtidyverse.Rmd
53 lines (47 loc) · 1.22 KB
/
F-rsimsumtidyverse.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
---
title: "rsimsum and the tidyverse"
author: "Alessandro Gasparini"
date: "`r Sys.Date()`"
output: rmarkdown::html_vignette
vignette: >
%\VignetteIndexEntry{rsimsum and the tidyverse}
%\VignetteEngine{knitr::rmarkdown}
%\VignetteEncoding{UTF-8}
---
```{r setup, include = FALSE}
options(width = 150)
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.align = "center", fig.height = 6, fig.width = 6,
out.width = "75%"
)
```
`rsimsum` plays nice with the tidyverse.
```{r packages}
library(rsimsum)
library(dplyr)
library(ggplot2)
library(knitr)
```
For instance, it is possible to chain functions using the piping operator `%>%` to obtain tables and plots with a single call:
```{r table}
data("MIsim", package = "rsimsum")
MIsim %>%
simsum(estvarname = "b", se = "se", methodvar = "method", true = 0.5) %>%
summary() %>%
tidy() %>%
kable()
```
```{r plot}
MIsim %>%
simsum(estvarname = "b", se = "se", methodvar = "method", true = 0.5) %>%
summary() %>%
tidy(stats = "bias") %>%
ggplot(aes(x = method, y = est, ymin = lower, ymax = upper)) +
geom_hline(yintercept = 0, color = "red", lty = "dashed") +
geom_point() +
geom_errorbar(width = 1 / 3) +
theme_bw() +
labs(x = "Method", y = "Bias")
```