Function to produce plots of the distribution of partial eta-squared.

```
plot_pes(
Fstat,
df1,
df2,
type = c("c", "cd"),
levels = c(0.68, 0.9, 0.95, 0.999)
)
```

## Arguments

- Fstat
The F-statistic from the F-test.

- df1
Degrees of freedom for the numerator.

- df2
Degrees of freedom for the denominator.

- type
Choose whether to plot a "consonance" function ("c"), consonance density ("cd"), or both (c("c","cd"); default option).

- levels
Numeric vector of confidence levels to display

## Value

Returns plot of the distribution of partial eta-squared

## Details

This function was created so that users could create consonance plots of partial eta-squared from ANOVA-level effects.
These types of plots are discussed by Schweder T, Hjort NL. (2016, ISBN:9781316445051) and Rafi Z, Greenland S. (2020) doi:10.1186/s12874-020-01105-9.