Skip to contents

This function creates a generic clinical significance plot by plotting the patients' pre intervention value on the x-axis and the post intervention score on the y-axis.

Usage

# S3 method for cs_statistical
plot(
  x,
  x_lab = "Pre",
  y_lab = "Post",
  color_lab = "Group",
  include_cutoff = TRUE,
  lower_limit,
  upper_limit,
  show,
  point_alpha = 1,
  overplotting = 0.02,
  ...
)

Arguments

x

An object of class cs_statistical

x_lab

String, x axis label. Default is "Pre".

y_lab

String, x axis label. Default is "Post".

color_lab

String, color label (if colors are displayed). Default is "Group"

include_cutoff

Logical, whether to include the population cutoff. Default is TRUE.

lower_limit

Numeric, lower plotting limit. Defaults to 2% smaller than minimum instrument score

upper_limit

Numeric, upper plotting limit. Defaults to 2% larger than maximum instrument score

show

Unquoted category name. You have several options to color different features. Available are

  • category (shows all categories at once)

  • clinical_pre (shows participants with clinical scores pre intervention)

  • functional_post (shows participants with functional scores post intervention)

  • unchanged (shows unchanged participants)

point_alpha

Numeric, transparency adjustment for points. A value between 0 and 1 where 1 corresponds to not transparent at all and 0 to fully transparent.

overplotting

Numeric, control amount of overplotting. Defaults to 0.02 (i.e., 2% of range between lower and upper limit).

...

Additional arguments

Value

A ggplot2 plot

Examples

cs_results <- antidepressants |>
  cs_statistical(
    patient,
    measurement,
    pre = "Before",
    mom_di,
    m_functional = 15,
    sd_functional = 8,
    cutoff_type = "c"
  )


# Plot the results "as is"
plot(cs_results)



# Change the axis labels
plot(cs_results, x_lab = "Before Intervention", y_lab = "After Intervention")



# Show the individual categories
plot(cs_results, show = category)



# Show groups as specified in the data
cs_results_grouped <- antidepressants |>
  cs_statistical(
    patient,
    measurement,
    pre = "Before",
    mom_di,
    m_functional = 15,
    sd_functional = 8,
    cutoff_type = "c",
    group = condition
  )

plot(cs_results_grouped)



# To avoid overplotting, generic ggplot2 code can be used to facet the plot
library(ggplot2)
plot(cs_results_grouped) +
  facet_wrap(~ group)



# Adjust the transparency of individual data points
plot(cs_results, point_alpha = 0.3)



# Control the overplotting
plot(cs_results, overplotting = 0.1)



# Or adjust the axis limits by hand
plot(cs_results, lower_limit = 0, upper_limit = 80)