plot_MADE.power {POMADE}

Plot function for a 'power' object

Description

Creates a faceted plot or plots for power analyses conducted with power_MADE.

Usage

## S3 method for class 'power'
plot_MADE(
  data,
  v_lines = NULL,
  legend_position = "bottom",
  color = TRUE,
  numbers = TRUE,
  number_size = 2.5,
  numbers_ynudge = 0,
  caption = TRUE,
  x_lab = NULL,
  x_breaks = NULL,
  x_limits = NULL,
  y_breaks = seq(0, 1, 0.2),
  y_limits = c(0, 1),
  y_expand = NULL,
  warning = TRUE,
  traffic_light_assumptions = NULL,
  power_min = NULL,
  expected_studies = NULL,
  model_comparison = FALSE,
  ...
)

Arguments

data

Data/object for which the plot should be made.

v_lines

Integer or vector to specify vertical line(s) in within each plot. Default is NULL.

legend_position

Character string to specify position of legend. Default is "bottom".

color

Logical indicating whether to use color in the plot(s). Default is TRUE.

numbers

Logical indicating whether to number the plots. Default is TRUE.

number_size

Integer value specifying the size of the (optional) plot numbers. Default is 2.5.

numbers_ynudge

Integer value for vertical nudge of the (optional) plot numbers.

caption

Logical indicating whether to include a caption with detailed information regarding the analysis. Default is TRUE.

x_lab

Title for the x-axis. If NULL (the default), the x_lab is specified automatically.

x_breaks

Optional vector to specify breaks on the x-axis. Default is NULL.

x_limits

Optional vector of length 2 to specify the limits of the x-axis. Default is NULL, which allows limits to be determined automatically from the data.

y_breaks

Optional vector to specify breaks on the y-axis.

y_limits

Optional vector of length 2 to specify the limits of the y-axis.

y_expand

Optional vector to expand the limits of the y-axis. Default is NULL.

warning

Logical indicating whether warnings should be returned when multiple models appear in the data. Default is TRUE.

traffic_light_assumptions

Optional logical to specify coloring of strips of the facet grids to emphasize assumptions about the likelihood the given analytical scenario. See Vembye, Pustejovsky, & Pigott (In preparation) for further details.

power_min

Either an integer specify a horizontal line or a length-2 vector to specify an interval, indicating a benchmark level of power (default is NULL).

expected_studies

Optional vector of length 2 specifying a range for the number of studies one expects to include in the meta-analysis. If specified, this interval will be shaded across facet_grip plots (default is NULL).

model_comparison

Logical indicating whether power estimates should be plotted across different working models for dependent effect size estimates (default is FALSE) instead of across values for the sampling correlation.

...

Additional arguments available for some classes of objects.

Details

In general, it can be rather difficult to guess/approximate the true model parameters and sample characteristics a priori. Calculating power under only a single set of assumptions can easily be misleading even if the true model and data structure only slightly diverge from the yielded data and model assumptions. To maximize the informativeness of the power approximations, Vembye, Pustejovsky, & Pigott (In preparation) suggest accommodating the uncertainty of the power approximations by reporting or plotting power estimates across a range of possible scenarios, which can be done using plot_MADE.power.

Value

A ggplot plot showing power across the expected number of studies, faceted by the between-study and within-study SDs, with different colors, lines, and shapes corresponding to different values of the assumed sample correlation. If model_comparison = TRUE, it returns a ggplot plot showing power across the expected number of studies, faceted by the between-study and within-study SDs, with different colors, lines, and shapes corresponding to different working models for dependent effect size estimates

References

Vembye, M. H., Pustejovsky, J. E., & Pigott, T. D. (In preparation). Conducting power analysis for meta-analysis of dependent effect sizes: Common guidelines and an introduction to the POMADE R package.

See Also

plot_MADE

Examples

power_dat <-
  power_MADE(
    J = c(50, 56),
    mu = 0.15,
    tau = 0.1,
    omega = 0.05,
    rho = 0,
    sigma2_dist = 4 / 200,
    n_ES_dist = 6
  )

power_example <-
  plot_MADE(
   data = power_dat,
   power_min = 0.8,
   expected_studies = c(52, 54),
   warning = FALSE,
   caption = TRUE,
   color = TRUE,
   model_comparison = FALSE,
   numbers = FALSE
   )

power_example