`R/Ancillary_Indicators.R`

`Probs.Rd`

Calculates mahalanobis distance and rejection of the Null operating model, used by wrapping function PRBcalc.

`Probs(indPPD, indData, alpha = 0.05, removedat = FALSE, removethresh = 0.05)`

- indPPD
A 3D array of results arising from running getind on an MSE of the Null operating model (type of data/stat (e.g. mean catches),time period (chunk), simulation)

- indData
A 3D array of results arising from running getind on an MSE of the Alternative operating model (type of data/stat (e.g. mean catches),time period (chunk), simulation)

- alpha
Positive fraction: rate of type I error, alpha

- removedat
Logical, should data not contributing to the mahalanobis distance be removed?

- removethresh
Positive fraction: the cumulative percentage of removed data (removedat=TRUE) that contribute to the mahalanobis distance

A list object. Position 1 is an array of the mahalanobis distances. Dimension 1 is length 2 for the Null OM (indPPD) and the alternative OM (indData). Dimension 2 is the time block (same length as indPPD dim 2). Dimension 3 is the simulation number (same length at indPPD dim 3.), Position 2 is a matrix (2 rows, ntimeblock columns) which is (row 1) alpha: the rate of false positives, and row 2 the power (1-beta) the rate of true positives

Carruthers and Hordyk 2018