Compute alarm threshold of Bernoulli cumulative sum control charts using simulation.
bcusum_crit_sim(L0, df, R0 = 1, RA = 2, m = 100, nc = 1,
jmax = 4, verbose = FALSE)Double. Prespecified in-control Average Run Length.
Data Frame. First column are Parsonnet Score values within a range of 0 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.
Integer. Number of digits for grid search.
Logical. 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 bcusum_crit_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 bcusum_arl_sim using mclapply.