plot_profile_verif: Plot verification scores for vertical profiles

In andrew-MET/harpVis: Visualisation functions for harp.

Plot verification scores for vertical profiles

Description

This is a wrapper for plot_point_verif that ensures the coordinate system is properly set up for plots of verification scores for vertical profiles.

Usage

plot_profile_verif(
  verif_data,
  score,
  y_axis = "p",
  lead_time = NA,
  plot_num_cases = TRUE,
  num_cases_position = c("right", "left"),
  ...
)

Arguments

verif_data	The verification data.
score	The score to be plotted. Must be unquoted.
y_axis	The y axis for the plot. Must be unquoted. Be default this is p, for pressure on the y axis, otherwise it should be set to the column name in verif_data that you wish to use for the y axis.
lead_time	The leadtime (or leadtimes) to plot. By default all lead times in the data are plotted. If more than one lead time is to be plotted, they must be separated by setting colour_by = leadtime and / or facet_by = vars(leadtime). See plot_point_verif for details of these arguments.
plot_num_cases	Logical of whether to inlclude a panel for the number of cases. Note that the number of cases will not be plotted if facet_by is set.
num_cases_position	The position of the number of cases panel relative to the score panel. For profile scores, this can only be "left", or "right".
...	Other arguments for plot_point_verif.

Value

A ggplot object of vertical profiles of verification scores.

Examples

# Not filtering or faceting results in overplotting
plot_profile_verif(verif_data_prf, bias)

# Facet by lead_time
plot_profile_verif(verif_data_prf, bias, facet_by = vars(lead_time))

# Remove rows with few cases using filter_by
plot_profile_verif(
  verif_data_prf,
  bias,
  facet_by = vars(lead_time),
  filter_by = vars(num_cases > 1)
)

# Filter to only a single panel
plot_profile_verif(verif_data_prf, bias, filter_by = vars(lead_time == 24))

andrew-MET/harpVis documentation built on March 11, 2024, 9:34 a.m.