ggstatsplot: ggplot2 Based Plots with Statistical Details

# <code>ggstatsplot</code>: <code>ggplot2</code> Based Plots with Statistical Details
## An Introductory Tutorial (0.9.0)
### Indrajeet Patil (Software Engineer, esqLABS)
### 2021-11-13

---

---

# Plan

---

- Why `ggstatsplot`?

- Primary functions

- Customizability

- Benefits

- Misconceptions

- Limitations

- Future plans

---

# Why *ggstatsplot*?

---

# Raison d'être

---

]

.content-box-yellow[
In short, `ggstatsplot` returns 
 
📊 .blue[information-rich] plots with .blue[statistical details] 
 
📝 suitable for .blue[faster] (exploratory) data analysis and scholarly publications
]

]
]

---

# Quickest overview: Publication (~ 10 min read)

---

Can be read [here](https://joss.theoj.org/papers/10.21105/joss.03167).

![](images/citation.png)

---

# Simpler/faster data analysis workflow

---

In a typical *exploratory* data analysis workflow, .blue[data visualization] and .blue[statistical modeling] are two
different phases: visualization informs modeling, and modeling can suggest a
different visualization, and so on and so forth.

💡 The central idea of `ggstatsplot` is simple: combine these two phases into one!

---

# Information-rich graphic is worth a thousand words

---

.footnote[[(Matejka & Fitzmaurice, *Autodesk Research*, 2017)](https://www.autodeskresearch.com/publications/samestats)]

.blue[Graphical] summaries can reveal problems not visible from .blue[numerical] statistics.
 
_I plotted my data and what I found surprised me!_ - BuzzFeed

---

# Ready-made plot = no customization

The .blue[grammar of graphics] is a powerful framework [(Wilkinson, 2011)](https://www.google.com/books/edition/_/iI1kcgAACAAJ?hl=en&sa=X&ved=2ahUKEwiGl8rJ2KztAhWyElkFHa8NAvkQre8FMBR6BAgMEAc) and can help you make *any* graphics fitting your specific data visualization needs! But...

`$\sum_{time}$` (Needed time ↑ + Likelihood to graphical explore data ↓) = Avoidance habit

---

# And a LOT more!
...but we will come back to that later 📌

Let's get started first!

---

# Installation

Install the stable version of `ggstatsplot` from 
[CRAN](https://cran.r-project.org/web/packages/ggstatsplot/index.html):

```r
install.packages("ggstatsplot")
```
--

You can get the development version of the package from
[Github](https://github.com/IndrajeetPatil/ggstatsplot):

```r
remotes::install_github("IndrajeetPatil/ggstatsplot")
```

Load the needed packages-

```r
library(ggstatsplot)
library(ggplot2)
```

---

# Primary functions

---

# Hypothesis about group differences

---

# ggbetweenstats - For between group comparisons

---

```r
ggbetweenstats(
  data = movies_long,
  x = mpaa,
  y = rating
)
```

- *t*-test if `2` groups
- ANOVA if `> 2` groups

📝 .blue[Defaults] return

✅ raw data + distributions 
✅ descriptive statistics 
✅ inferential statistics 
✅ effect size + CIs 
✅ pairwise comparisons 
✅ Bayesian hypothesis-testing 
✅ Bayesian estimation
]

]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggbetweenstats_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# ggbetweenstats - pairwise comparisons

---

```r
ggbetweenstats(
  data = movies_long,
  x = mpaa,
  y = rating,
* type = "np",
* pairwise.display = "ns"
)
```

✅ `"p"` &nbsp;&nbsp;→ **parametric** 
✅ `"np"` → **non-parametric** 
✅ `"r"` &nbsp;&nbsp;→ **robust** 
✅ `"bf"` → **Bayesian**

Changing pairwise comparisons displayed

ℹ️ `"ns"` nbsp;&→ only **non-significant** 
ℹ️ `"s"` nbsp;&nbsp;→→ only **significant** 
ℹ️ `"all"` → **all**
]

]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggbetweenstats_2-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# ggbetweenstats - outlier tagging

---

```r
ggbetweenstats(
  data = movies_long,
  x = mpaa,
  y = rating,
  type = "r",
  pairwise.comparisons = FALSE,
* outlier.tagging = TRUE,
* outlier.label = title
)
```

.font70[
[Tukey's fences](https://en.wikipedia.org/wiki/Outlier#Tukey's_fences) method using 
interquartile range flags outliers.

Centrality measures

✅ `"p"` &nbsp;&nbsp;→ `$\mu_{mean}$` 
✅ `"np"` → `$\mu_{median}$` 
✅ `"r"` &nbsp;&nbsp;→ `$\mu_{trimmed}$` 
✅ `"bf"` → `$\mu_{MAP}$`
]

]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggbetweenstats_3-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# ggwithinstats - repeated measures equivalent

---

```r
ggwithinstats(
  data = WRS2::WineTasting,
  x = Wine,
  y = Taste
)
```

Changing the `type` of test

✅ `"p"` &nbsp;&nbsp;→ **parametric** 
✅ `"np"` → **non-parametric** 
✅ `"r"` &nbsp;&nbsp;→ **robust** 
✅ `"bf"` → **Bayesian**
]
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggwithinstats_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# gghistostats - Distribution of a numeric variable

---

```r
gghistostats(
  data = movies_long,
  x = budget,
* test.value = 30
)
```

✅ counts + proportion for bins 
✅ descriptive statistics 
✅ inferential statistics 
✅ effect size + CIs 
✅ Bayesian hypothesis-testing 
✅ Bayesian estimation

Centrality measures

✅ `"p"` &nbsp;&nbsp;→ `$\mu_{mean}$` 
✅ `"np"` → `$\mu_{median}$` 
✅ `"r"` &nbsp;&nbsp;→ `$\mu_{trimmed}$` 
✅ `"bf"` → `$\mu_{MAP}$`
]
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/gghistostats_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# ggdotplotstats - Labeled numeric variable

---

```r
ggdotplotstats(
  data = movies_long,
  x = budget,
  y = genre,
* test.value = 30
)
```

✅ descriptive statistics 
✅ inferential statistics 
✅ effect size + CIs 
✅ Bayesian hypothesis-testing 
✅ Bayesian estimation

Centrality measures

✅ `"p"` &nbsp;&nbsp;→ `$\mu_{mean}$` 
✅ `"np"` → `$\mu_{median}$` 
✅ `"r"` &nbsp;&nbsp;→ `$\mu_{trimmed}$` 
✅ `"bf"` → `$\mu_{MAP}$`
]
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggdotplotstats_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# Hypothesis about correlation

---

# ggscatterstats - Two numeric variables

---

```r
ggscatterstats(
  data = movies_long,
  x = budget,
  y = rating
)
```

✅ joint distribution 
✅ marginal distributions 
✅ inferential statistics 
✅ effect size + CIs 
✅ Bayesian hypothesis-testing 
✅ Bayesian estimation

Changing the `type` of test

✅ `"p"` &nbsp;&nbsp;→ **parametric** 
✅ `"np"` → **non-parametric** 
✅ `"r"` &nbsp;&nbsp;→ **robust** 
✅ `"bf"` → **Bayesian**
]
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggscatterstats_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# ggscatterstats - conditional point tagging

---

```r
ggscatterstats(
  data = movies_long,
  x = budget,
  y = rating,
  type = "r",
* label.var = title,
* label.expression = budget > 150
  & rating > 7.5
)
```
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggscatterstats_2-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# ggcorrmat - multiple numeric variables

---

```r
ggcorrmat(dplyr::starwars)
```

✅ effect size + significance 
✅ careful handling of `NA`s

Changing the `type` of test

✅ `"p"` &nbsp;&nbsp;→ **parametric** 
✅ `"np"` → **non-parametric** 
✅ `"r"` &nbsp;&nbsp;→ **robust** 
✅ `"bf"` → **Bayesian**

]
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggcorrmat_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# ggcorrmat - getting a dataframe

---

```r
library(ggplot2) # for data
options(digits = 2)

ggcorrmat(
  data = dplyr::select(msleep, sleep_rem, awake, brainwt),
  type = "bayes",
* output = "dataframe"
)
```

<table>
 <thead>
 <tr>
 <th style="text-align:left;"> parameter1 </th>
 <th style="text-align:left;"> parameter2 </th>
 <th style="text-align:right;"> estimate </th>
 <th style="text-align:right;"> conf.level </th>
 <th style="text-align:right;"> conf.low </th>
 <th style="text-align:right;"> conf.high </th>
 <th style="text-align:right;"> pd </th>
 <th style="text-align:right;"> rope.percentage </th>
 <th style="text-align:left;"> prior.distribution </th>
 <th style="text-align:right;"> prior.location </th>
 <th style="text-align:right;"> prior.scale </th>
 <th style="text-align:right;"> bayes.factor </th>
 <th style="text-align:left;"> method </th>
 <th style="text-align:right;"> n.obs </th>
 </tr>
 </thead>
<tbody>
 <tr>
 <td style="text-align:left;"> sleep_rem </td>
 <td style="text-align:left;"> awake </td>
 <td style="text-align:right;"> -0.73 </td>
 <td style="text-align:right;"> 0.95 </td>
 <td style="text-align:right;"> -0.84 </td>
 <td style="text-align:right;"> -0.60 </td>
 <td style="text-align:right;"> 1.00 </td>
 <td style="text-align:right;"> 0.00 </td>
 <td style="text-align:left;"> beta </td>
 <td style="text-align:right;"> 1.4 </td>
 <td style="text-align:right;"> 1.4 </td>
 <td style="text-align:right;"> 3.0e+09 </td>
 <td style="text-align:left;"> Bayesian Pearson correlation </td>
 <td style="text-align:right;"> 61 </td>
 </tr>
 <tr>
 <td style="text-align:left;"> sleep_rem </td>
 <td style="text-align:left;"> brainwt </td>
 <td style="text-align:right;"> -0.21 </td>
 <td style="text-align:right;"> 0.95 </td>
 <td style="text-align:right;"> -0.45 </td>
 <td style="text-align:right;"> 0.06 </td>
 <td style="text-align:right;"> 0.93 </td>
 <td style="text-align:right;"> 0.21 </td>
 <td style="text-align:left;"> beta </td>
 <td style="text-align:right;"> 1.4 </td>
 <td style="text-align:right;"> 1.4 </td>
 <td style="text-align:right;"> 6.5e-01 </td>
 <td style="text-align:left;"> Bayesian Pearson correlation </td>
 <td style="text-align:right;"> 48 </td>
 </tr>
 <tr>
 <td style="text-align:left;"> awake </td>
 <td style="text-align:left;"> brainwt </td>
 <td style="text-align:right;"> 0.34 </td>
 <td style="text-align:right;"> 0.95 </td>
 <td style="text-align:right;"> 0.09 </td>
 <td style="text-align:right;"> 0.55 </td>
 <td style="text-align:right;"> 0.99 </td>
 <td style="text-align:right;"> 0.03 </td>
 <td style="text-align:left;"> beta </td>
 <td style="text-align:right;"> 1.4 </td>
 <td style="text-align:right;"> 1.4 </td>
 <td style="text-align:right;"> 7.3e+00 </td>
 <td style="text-align:left;"> Bayesian Pearson correlation </td>
 <td style="text-align:right;"> 56 </td>
 </tr>
</tbody>
</table>

</div>
]

---

# Hypothesis of composition of categorical variables

---

# ggpiestats - association between categorical variables

---

```r
# let's use subset of data
ggpiestats(
  data = dplyr::filter(
    movies_long,
    genre %in% c("Drama", "Comedy")
  ),
  x = mpaa,
  y = genre
)
```

✅ descriptive statistics 
✅ inferential statistics 
✅ effect size + CIs 
✅ Goodness-of-fit tests 
✅ Bayesian hypothesis-testing 
✅ Bayesian estimation

]
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggpiestats_2-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# ggbarstats - association between categorical variables

---

```r
ggbarstats(
  data = dplyr::filter(
    movies_long,
    genre %in% c("Drama", "Comedy")
  ),
  x = mpaa,
  y = genre,
* label = "both"
)
```

✅ descriptive statistics 
✅ inferential statistics 
✅ effect size + CIs 
✅ Goodness-of-fit tests 
✅ Bayesian hypothesis-testing 
✅ Bayesian estimation 
]
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggbarstats_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# Hypothesis about regression coefficients

---

# ggcoefstats

---

```r
# model
mod <- lm(
 formula = rating ~ mpaa,
 data = movies_long
)

# plot
ggcoefstats(mod)
```

✅ estimate + CIs 
✅ inferential statistics ($t$, `$z$`, `$F$`, `$\chi^2$`) 
✅ model fit indices (AIC + BIC)
]

]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggcoefstats_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# ggcoefstats: Supported models

---

.font60[
`aareg`, `afex_aov`, `anova`, `anova.mlm`, `anova`, `aov`, `aovlist`, `Arima`,
`bam`, `bayesx`, `bayesGARCH`, `bayesQR`, `BBmm`, `BBreg`, `bcplm`, `betamfx`,
`betaor`, `BFBayesFactor`, `BGGM`, `bglmerMod`, `bife`, `bigglm`, `biglm`,
`blavaan`, `bmlm`, `blmerMod`, `blrm`, `bracl`, `brglm`, `brglm2`, `brmsfit`,
`brmultinom`, `btergm`, `cch`, `censReg`, `cgam`, `cgamm`, `cglm`, `clm`,
`clm2`, `clmm`, `clmm2`, `coeftest`, `complmrob`, `confusionMatrix`, `coxme`,
`coxph`, `coxr`, `coxph.penal`, `cpglm`, `cpglmm`, `crch`, `crq`, `crr`,
`DirichReg`, `drc`, `eglm`, `elm`, `emmGrid`, `epi.2by2`, `ergm`, `feis`,
`felm`, `fitdistr`, `fixest`, `flexsurvreg`, `gam`, `Gam`, `gamlss`, `garch`,
`geeglm`, `gjrm`, `glmc`, `glmerMod`, `glmmTMB`, `gls`, `glht`, `glm`, `glmm`,
`glmmadmb`, `glmmPQL`, `glmRob`, `glmrob`, `glmx`, `gmm`, `HLfit`, `hurdle`,
`ivFixed`, `ivprobit`, `ivreg`, `iv_robust`, `lavaan`, `lm`, `lm.beta`,
`lmerMod`, `lmerModLmerTest`, `lmodel2`, `lmRob`, `lmrob`, `lm_robust`,
`logitmfx`, `logitor`, `logitsf`, `LORgee`, `lqm`, `lqmm`, `lrm`, `manova`,
`maov`, `margins`, `mcmc`, `mcmc.list`, `MCMCglmm`, `mclogit`, `mice`,
`mmclogit`, `mediate`, `metafor`, `merMod`, `merModList`, `metaplus`, `mhurdle`,
`mixor`, `mjoint`, `mle2`, `mlm`, `multinom`, `mvord`, `negbin`, `negbinmfx`,
`negbinirr`, `nlmerMod`, `nlrq`, `nlreg`, `nls`, `orcutt`, `orm`, `plm`,
`poissonmfx`, `poissonirr`, `polr`, `probitmfx`, `ridgelm`, `riskRegression`,
`rjags`, `rlm`, `rlmerMod`, `robmixglm`, `rq`, `rqs`, `rqss`, `rrvglm`, `scam`,
`selection`, `semLm`, `semLme`, `slm`, `speedglm`, `speedlm`, `stanfit`,
`stanreg`, `summary.lm`, `survreg`, `svyglm`, `svy_vglm`, `svyolr`, `tobit`,
`truncreg`, `varest`, `vgam`, `vglm`, `wbgee`, `wblm`, `zeroinfl`, etc.
]

Thanks to [`easystats`](https://easystats.github.io/easystats/)! 🙏

![](images/easystats.jpg)

---

# *grouped_* variants of all functions
Running the same function for 
all levels of a single grouping variable

---

# *grouped_* functions

---

```r
grouped_ggpiestats(
  data = mtcars,
  x = cyl,
* grouping.var = am
)
```

.font70[
Available `grouped_` variants
- `grouped_ggbetweenstats`
- `grouped_ggwithinstats`
- `grouped_gghistostats`
- `grouped_ggdotplotstats`
- `grouped_ggscatterstats`
- `grouped_ggcorrmat`
- `grouped_ggpiestats`
- `grouped_ggbarstats`
]
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/grouped_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# Customizability of *ggstatsplot*
"What if I don't like the default plots?" 🤔

---

# Changing aesthetics (themes + palettes) 🖼🎨

---

Aesthetic preferences are not an excuse to not use `ggstatsplot`! 😻

```r
ggbetweenstats(
  data = movies_long,
  x = mpaa,
  y = rating,
* ggtheme = hrbrthemes::theme_ipsum_tw(),
* palette = "Darjeeling2",
* package = "wesanderson"
)
```

.font70[
The default palette 
is .blue[[colorblind-friendly](https://cran.r-project.org/web/packages/colorBlindness/vignettes/colorBlindness.html).
]
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggbetweenstats_4-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# Further modification with *ggplot2* 🛠

---

You can modify `ggstatsplot` plots further using `ggplot2` functions. 🎉

```r
ggbetweenstats(
  data = mtcars,
  x = am,
  y = wt,
  type = "bayes"
) +
* scale_y_continuous(sec.axis = dup_axis())
```

]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/ggbetweenstats_5-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# Too much information 🙈

---

`ggstatsplot` can be used to get .blue[only plots].

```r
# using `ggstatsplot` only for plot
ggbetweenstats(
  data = iris,
  x = Species,
  y = Sepal.Length,
  # turn off centrality measure
* centrality.plotting = FALSE,
  # turn off statistical analysis
* results.subtitle = FALSE,
  # turn off Bayesian message
* bf.message = FALSE,
  # turn off pairwise comparisons
* pairwise.comparisons = FALSE
)
```
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/only_plot-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# Expressions for custom plots 🀄️

---

`ggstatsplot` can be used to get .blue[only expressions].

```r
# using `ggstatsplot` for stats
results <-
 ggpiestats(
 data = Titanic_full,
 x = Survived,
 y = Sex,
* output = "subtitle"
 )

# using `ggiraphExtra` for plot
*ggiraphExtra::ggSpine(
  data = Titanic_full,
  aes(x = Sex, fill = Survived),
  addlabel = TRUE,
  interactive = FALSE
) +
* labs(subtitle = results)
```
]

.right-plot[
<img src="ggstatsplot_presentation_files/figure-html/subtitle_1-1.png" width="100%" style="display: block; margin: auto;" />
]

---

# Only dataframes

---

[`statsExpressions`](https://indrajeetpatil.github.io/statsExpressions/), statistical processing backend for `ggstatsplot`, can provide .blue[only dataframes].

]

```r
library(statsExpressions)

# for example
one_sample_test(
  data = mtcars,
  x = wt,
  test.value = 3
)
```
]

---

# Why use *ggstatsplot*? 👍️

---

# Supports different statistical approaches

Functions | Description | Parametric | Non-parametric | Robust | Bayesian
------- | ------------------ | ---- | ----- | ----| ----- 
`ggbetweenstats` | Between group comparisons | ✅ | ✅ | ✅ | ✅
`ggwithinstats` | Within group comparisons | ✅ | ✅ | ✅ | ✅
`gghistostats`, `ggdotplotstats` | Distribution of a numeric variable | ✅ | ✅ | ✅ | ✅
`ggcorrmat` | Correlation matrix | ✅ | ✅ | ✅ | ✅
`ggscatterstats` | Correlation between two variables | ✅ | ✅ | ✅ | ✅
`ggpiestats`, `ggbarstats` | Association between categorical variables | ✅ | `NA` | `NA` | ✅
`ggpiestats`, `ggbarstats` | Equal proportions for categorical variable levels | ✅ | `NA` | `NA` | ✅
`ggcoefstats` | Regression modeling | ✅ | ✅ | ✅ | ✅
`ggcoefstats` | Random-effects meta-analysis | ✅ | `NA` | ✅ | ✅

---

# Alternative workflow to the following

📦 to carry out the test (e.g. `stats`) 
📦 computing effect size + CIs (e.g. `effectsize`) 
📦 for descriptives (e.g. `skimr`) 
📦 pairwise comparisons (e.g. `multcomp`) 
📦 Bayesian hypothesis testing (e.g. `BayesFactor`) 
📦 Bayesian estimation (e.g. `bayestestR`) 
📦 . 
📦 . 
]

]

🤔 accepts dataframe, vectors, matrix? 
🤔 long/wide format data? 
🤔 works with `NA`s? 
🤔 returns list, dataframe, arrays? 
🤔 works with tibbles? 
🤔 . 
🤔 .

]

---

# Toggling between statistical approaches 🔀

**.blue[Parametric]**

```r
# anova
ggbetweenstats(
  data = mtcars,
  x = cyl,
  y = wt,
* type = "p"
)

# correlation analysis
ggscatterstats(
  data = mtcars,
  x = wt,
  y = mpg,
* type = "p"
)

# t-test
gghistostats(
  data = mtcars,
  x = wt,
  test.value = 2,
* type = "p"
)
```
]

**.orange[Non-parametric]**

```r
# anova
ggbetweenstats(
  data = mtcars,
  x = cyl,
  y = wt,
* type = "np"
)

# correlation analysis
ggscatterstats(
  data = mtcars,
  x = wt,
  y = mpg,
* type = "np"
)

# t-test
gghistostats(
  data = mtcars,
  x = wt,
  test.value = 2,
* type = "np"
)
```
]

---

# Results *in context* of the underlying data 🕵️

**Standard approach**

Pearson's correlation test revealed that, across 142 participants, variable `x`
was negatively correlated with variable `y`: `$t(140)=-0.76, p=.446$`. The
effect size `$(r=-0.06, 95\% CI [-.23,.10])$` was small, as per Cohen’s (1988)
conventions. The Bayes Factor for the same analysis revealed that the data were
5.81 times more probable under the null hypothesis as
compared to the alternative hypothesis. This can be considered moderate evidence
(Jeffreys, 1961) in favor of the null hypothesis (absence of any correlation
between `x` and `y`).

]

**`ggstatsplot` approach**

![](images/after_ggstats.PNG)

]

---

# Best practices in statistical reporting 🏆

The expression template tries to follow the gold standard for statistical
reporting.

![](images/stats_reporting_format.png)

---

# Avoiding reporting errors

.content-box-green[
 "half of all published psychology papers that use NHST contained at least one
*p*-value that was inconsistent with its test statistic and degrees of freedom.
One in eight papers contained a grossly inconsistent *p*-value that may have
affected the statistical conclusion"

[(Nuijten et al., *Behavior Research Methods*, 2016)](https://link.springer.com/article/10.3758/s13428-015-0664-2)
]

Since the plot and the statistical analysis are yoked together, the chances of
making an error in reporting the results are minimized.

No need to worry about updating figures and statistical details **separately**. 🔗

---

# Making sense of null results

`$p > 0.05$`: The null hypothesis (`H0`) can't be rejected

But can it be **accepted**?! Null Hypothesis Significance Testing 🤫

.content-box-green[
"In 72% of cases, nonsignificant results were misinterpreted, in that the
authors inferred that the effect was absent. A Bayesian reanalysis revealed
that fewer than 5% of the nonsignificant findings provided strong evidence
(i.e., `$BF_{01} > 10$`) in favor of the null hypothesis over the alternative
hypothesis."

[(Aczel et al., *AMPPS*, 2018)](https://journals.sagepub.com/doi/pdf/10.1177/2515245918773742)
]

Juxtaposing frequentist and Bayesian statistics for the same analysis helps to
properly interpret the null results.

---

# A few other benefits

---

.content-box-green[
Minimal code needed (`data`, `x`, `y`): minimizes chances of error + tidy scripts. 💅
]

---

# No more excuses not to explore data! 😉

---

In summary, the `ggstatsplot` approach- 
 
(*a*) avoids errors in statistical reporting, 
 
(*b*) highlights the importance of the effect by providing effect size measures by default, 
 
(*c*) provides an easy way to evaluate *absence* of an effect using Bayesian framework, 
 
(*d*) demands to evaluate statistical analysis in the context of the underlying
data, 
 
and is (*e*) easy and (*f*) simple enough that somebody with little coding
experience can use it without making an error.

]

---

# Misconceptions, limitations, future plans

---

# Misconceptions: This package is...

---

❌ an alternative to learning `ggplot2` 
--
✅ (the more you know `ggplot2`, the better you can modify the
defaults to your liking)

❌ meant to be used in talks/presentations 
--
✅ (defaults too complicated for effectively communicating
results in time-constrained presentation settings, e.g. conference talks)

❌ only relevant when used in publications 
--
✅ not necessary; can also be useful *only* during exploratory phase

❌ the only game in town 
--
✅ (excellent GUI open-source softwares: [JASP](https://jasp-stats.org/) and [jamovi](https://www.jamovi.org/))

---

# Limitations of *ggstatsplot* 👎️

---

.content-box-red[
Limited no. of **plots** and **statistical tests** available
(no support for multivariate statistics, time series, etc.). This will
**always** be the case. 🤷
]

.content-box-red[
Expects a non-trivial level of statistical proficiency (but
plots without statistics can still be useful).
]

.content-box-red[
**Faceting** does not work (since there are no corresponding
`geom_` s). For the same reason, plots are not `gganimate`-friendly.
]

---

# Overcoming these limitations 👥

---

Contributions (big or small) welcome!

![](images/needs_you.jpg)

]

- Raise issues about bugs/features 🐞

- Cite if used in a publication 📝

- Star on GitHub (increases visibility) ⭐

- Review code 🕵

- Add new functionality 👩‍💻
]

]

---

# Future plans 📅

---

.pull-left[
.content-box-purple[
**New plots**
- density plot
- connected scatterplot
- sankey diagram
- bar plot

**New statistical tests**
- complex (multiway) ANOVA
- survival analysis
- cluster analysis
]
]

![](images/taytay.gif)

]

---

# Acknowledgments

Developer friends 🙌

[Daniel Lüdecke](https://github.com/strengejacke), [Dominique Makowski](https://github.com/DominiqueMakowski), [Mattan S. Ben-Shachar](https://github.com/mattansb), [Brenton M. Wiernik](https://github.com/bwiernik)

Support 💰

[Mina Cikara](http://www.intergroupneurosciencelaboratory.com/), [Fiery Cushman](http://cushmanlab.fas.harvard.edu/index.php), [Iyad Rahwan](https://rahwan.me/)

Community 🙏

Contributors to *ggstatsplot* &
*rstats* users and developers

---

# Find me at...

[🐦 @patilindrajeets](http://twitter.com/patilindrajeets)

[💻 @IndrajeetPatil](http://github.com/IndrajeetPatil)

[🔗 https://sites.google.com/site/indrajeetspatilmorality/](https://sites.google.com/site/indrajeetspatilmorality/)

[📭 patilindrajeet.science@gmail.com](mailto:patilindrajeet.science@gmail.com)

]

---

# The End 👋

For more information, see-

<https://indrajeetpatil.github.io/ggstatsplot/>

To access code for these slides, see-

<https://github.com/IndrajeetPatil/ggstatsplot_slides/>