SSmase: SSmase() computes MASE for one-step ahead hindcasting...

View source: R/SSmase.R

SSmaseR Documentation

SSmase() computes MASE for one-step ahead hindcasting cross-validations of indices

Description

MASE for one-step ahead hindcasting cross-validations and computes MASE from prediction redisuals. MASE is calculated the average ratio of mean absolute error (MAE) of prediction residuals (MAE.PR) and Naive Predictions (MAE.base) MASE.adj sets the MAE.base to a minimum MAE.base.adj (default=0.1) MASE.adj allow passing (MASE<1) if MAE.PE < 0.1 and thus accurate, when obs show extremely little variation

Usage

SSmase(
  retroSummary,
  quants = c("cpue", "len", "age"),
  Season = "default",
  models = "all",
  endyrvec = "default",
  indexselect = NULL,
  MAE.base.adj = 0.1,
  residuals = FALSE,
  verbose = FALSE
)

Arguments

retroSummary

List created by r4ss::SSsummarize()

quants

data type c("cpue","len","age)

Season

option to specify Season - Default uses first available, i.e. usual Seas = 1

models

Optional subset of the models described in r4ss function summaryoutput(). Either "all" or a vector of numbers indicating columns in summary tables.

endyrvec

Optional single year or vector of years representing the final year of values to show for each model. By default it is set to the ending year specified in each model.

indexselect

= Vector of fleet numbers for each model for which to compare

MAE.base.adj

minimum MASE demoninator (naive predictions) for MASE.adj (default = 0.1)

residuals

if TRUE, outputs individual prediction and naive residuals

verbose

Report progress to R GUI?

indexfleets

CHECK IF NEEDED or how to adjust indexfleets

Value

MASE and hcxval statistic

Author(s)

Henning Winker (JRC-EC) and Laurence Kell (Sea++)


jabbamodel/ss3diags documentation built on Oct. 8, 2024, 11:49 p.m.