Calculates mahalanobis distance and rejection of the Null operating model
Source:R/Ancillary_Indicators.R
Probs.RdCalculates mahalanobis distance and rejection of the Null operating model, used by wrapping function PRBcalc.
Arguments
- 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
Value
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