Compute alarm threshold of the Bernoulli CUSUM control charts using simulation.
cusum_arl_h_sim(L0, df, R0 = 1, RA = 2, m = 100, nc = 1,
verbose = FALSE)double. Prespecified in-control Average Run Length.
DataFrame. First column are Parsonnet Score values within a range of zero to 100 representing the preoperative patient risk. The second column are binary (0/1) outcome values of each operation.
double. Odds ratio of death under the null hypotheses.
double. Odds ratio of death under the alternative hypotheses. Detecting deterioration in performance with increased mortality risk by doubling the odds Ratio RA=2. Detecting improvement in performance with decreased mortality risk by halving the odds ratio of death RA=1/2.
integer. Number of simulation runs.
integer. Number of cores.
boolean. If TRUE verbose output is included, if FALSE a quiet calculation of h is done.
Returns a single value which is the control limit h for a given in-control ARL.
The function cusum_arl_h_sim determines the control limit for given in-control ARL (L0) by applying a
multi-stage search procedure which includes secant rule and the parallel version of cusum_arl_sim
using mclapply.