Skip to contents

Brier score and its decomposition for an ensemble.

Usage

ens_brier(
  .fcst,
  parameter,
  thresholds,
  groupings = "lead_time",
  climatology = "sample",
  rel_probs = NA,
  num_ref_members = NA,
  keep_score = c("both", "brier", "reliability"),
  show_progress = TRUE,
  ...
)

# S3 method for harp_ens_point_df
ens_brier(
  .fcst,
  parameter,
  thresholds,
  groupings = "lead_time",
  climatology = "sample",
  rel_probs = NA,
  num_ref_members = NA,
  keep_score = "both",
  show_progress = TRUE,
  fcst_model = NULL,
  ...
)

Arguments

.fcst

A harp_df or harp_list object with tables that have a column for observations, or a single forecast table.

parameter

The name of the column for the observations data. Can be the column name, quoted, or unquoted. If a variable it should be embraced - i.e. wrapped in {{}}.

thresholds

A numeric vector of thresholds for which to compute the threshold based scores. Set to NULL (the default) to only compute summary scores.

groupings

The groups for which to compute the scores. See group_by for more information of how grouping works.

climatology

The climatology to use for the Brier Skill Score. Can be "sample" for the sample climatology (the default), a named list with elements eps_model and member to use a member of an eps model in the harp_fcst object for the climatology, or a data frame with columns for threshold and climatology and also optionally lead_time.

rel_probs

Probabilities to use for reliability diagrams. Set to NA (the default) to select automatically.

num_ref_members

For "fair" scores, the score is scaled to be valid for this number of ensemble members. Set to NA (the default) to not modify the score.

show_progress

Logical - whether to show progress bars. Defaults to TRUE.

...

Reserved for methods.

fcst_model

The name of the forecast model to use in the fcst_model column of the output. If the function is dispatched on a harp_list object, the names of the harp_list are automatically used.

Value

A data frame with data grouped for the groupings column(s) and columns for brier_score, brier_skill_score and the deomposition of the brier score - brier_score_reliability, brier_score_resolution and brier_score_uncertainty.