Skip to contents

plot_PPR_data_single_treatment() creates a categorical scatter plot with experimental state (i.e. baseline/before and after) on the x-axis and the paired-pulse ratio (PPR) on the y-axis. There are also lines connecting the "before" data point to the "after" data point for each letter.

Usage

plot_PPR_data_single_treatment(
  data,
  plot_treatment = "Control",
  plot_category = 2,
  included_sexes = "both",
  facet_by_sex = "no",
  male_label = "Male",
  female_label = "Female",
  baseline_label = "Baseline",
  post_hormone_label = "Post-hormone",
  y_axis_title = "PPR",
  test_type,
  plot_y_max = 3,
  map_signif_level_values = F,
  geom_signif_family = "",
  geom_signif_text_size = 5,
  large_axis_text = "no",
  mean_line_thickness = 1.2,
  mean_point_size = 2.5,
  geom_signif_size = 0.4,
  treatment_colour_theme,
  theme_options,
  save_plot_png = "no",
  ggplot_theme = patchclampplotteR_theme()
)

Arguments

data

Paired pulse ratio data generated from make_PPR_data().

plot_treatment

A character value specifying the treatment you would like to plot (e.g. "Control"). plot_treatment represents antagonists that were present on the brain slice, or the animals were fasted, etc.

plot_category

A numeric value specifying the category, which can be used to differentiate different protocol types. In the sample dataset for this package, plot_category == 2 represents experiments where insulin was applied continuously after a 5-minute baseline period.

included_sexes

A character value ("both", "male" or "female"). Useful if you want to have a plot with data from one sex only. Defaults to "both". If you choose a single sex, the resulting plot will have "-males-only" or "-females-only" in the file name. WARNING!! If you choose "male" or "female", you MUST ensure that the t_test_df contains data that has been filtered to only include one sex. Otherwise, the significance stars will represent both sexes and it will be inaccurate.

facet_by_sex

A character value ("yes" or "no") describing if the plots should be faceted by sex. This is only available if included_sexes is "both". The resulting plot will be split in two, with male data on the left and female data on the right.

male_label

A character value used to describe how males are encoded in the sex column of the dataframe used in data. This MUST match the value for male data in the sex column, and it must be consistent across data sheets. Defaults to "Male".

female_label

A character value used to describe how females are encoded in the sex column of the dataframe used in data. This MUST match the value for female data in the sex column, and it must be consistent across data sheets. This must be consistent in all data sheets. Defaults to "Female".

baseline_label

A character value for the x-axis label applied to the pre-hormone state. Defaults to "Baseline".

post_hormone_label

A character value for x-axis label applied to the post-hormone or post-protocol state. Defaults to "Post-hormone" but you will likely change this to the hormone or protocol name.

y_axis_title

A character value describing the y-axis title text. Defaults to "PPR" but could be expanded (e.g. "Paired pulse ratio").

test_type

A character (must be "wilcox.test", "t.test" or "none") describing the statistical model used to create a significance bracket comparing the pre- and post-hormone groups.

plot_y_max

A numeric value describing the maximum value of the y-axis. Defaults to 3.

map_signif_level_values

A TRUE/FALSE value or a list of character values for mapping p-values. If TRUE, p-values will be mapped with asterisks (e.g. \* for p < 0.05, for p < 0.01). If FALSE, raw p-values will display. You can also insert a list of custom mappings or a function. For example, use map_signif_level_values = function(p) if (p < 0.1) {round(p, 3)} else {"ns"} to only display the p-values when they are below 0.1.

geom_signif_family

A character value describing the font family used for the p-value annotations used by ggsignif::geom_signif(). Defaults to "" (empty value, will be replaced with default system font), but can be replaced with a named font. Use a package like extrafont to load system fonts into R.

geom_signif_text_size

A numeric value describing the size of the text annotations (significance stars or p-values) on the plot. Defaults to 8.

large_axis_text

A character ("yes" or "no"). If "yes", a ggplot theme layer will be applied which increases the size of the axis text.

mean_line_thickness

A numeric value describing the thickness of the line used to indicate the mean for a group. Defaults to 1.2.

mean_point_size

A numeric value describing the size of the points used to indicate the means. Defaults to 2.5.

geom_signif_size

A numeric value describing the size of the geom_signif bracket size. Defaults to 0.4, which is a good thickness for most applications.

treatment_colour_theme

A dataframe containing treatment names and their associated colours as hex values. See sample_treatment_names_and_colours for an example of what this dataframe should look like.

theme_options

A dataframe containing theme options. See sample_theme_options for an example of what this dataframe should look like.

save_plot_png

A character ("yes" or "no"). If "yes", the plot will be saved as a .png using ggsave. The filepath depends on the current type, but they will all go in subfolders below Figures/ in your project directory.

ggplot_theme

The name of a ggplot theme or your custom theme. This will be added as a layer to a ggplot object. The default is patchclampplotteR_theme(), but other valid entries include theme_bw(), theme_classic() or the name of a custom ggplot theme stored as an object.

Value

A ggplot object. If save_plot_png == "yes", it will also generate a .png file in the folder Figures/Evoked-currents/PPR relative to the project directory. The treatment will be included in the filename.

Details

If you specify a test_type, the function will perform a paired t-test or paired wilcox test and add brackets with significance stars through ggsignif::geom_signif().

See also

plot_PPR_data_multiple_treatments() to plot changes in PPR for multiple treatments. See make_PPR_data() for the function used to create the PPR data.

Examples

plot_PPR_data_single_treatment(
  data = sample_PPR_df,
  plot_treatment = "Control",
  plot_category = 2,
  baseline_label = "Baseline",
  post_hormone_label = "Insulin",
  included_sexes = "both",
  facet_by_sex = "no",
  test_type = "t.test",
  large_axis_text = "no",
  treatment_colour_theme = sample_treatment_names_and_colours,
  theme_options = sample_theme_options,
  save_plot_png = "no"
)


# Facet by sex

plot_PPR_data_single_treatment(
  data = sample_PPR_df,
  plot_treatment = "Control",
  plot_category = 2,
  baseline_label = "Baseline",
  post_hormone_label = "Insulin",
  included_sexes = "both",
  facet_by_sex = "yes",
  test_type = "t.test",
  large_axis_text = "no",
  treatment_colour_theme = sample_treatment_names_and_colours,
  theme_options = sample_theme_options,
  save_plot_png = "no"
)
#> Warning: Computation failed in `stat_signif()`.
#> Caused by error in `t.test.default()`:
#> ! not enough 'x' observations