vlad (version 0.2.2)

racusum_betabinomial_crit_sim: Compute alarm threshold of RA-CUSUM control charts using simulation

Description

Compute alarm threshold of risk-adjusted cumulative sum control charts using simulation.

Usage

racusum_betabinomial_crit_sim(L0, shape1, shape2, coeff, RA = 2, rs = 71,
  RQ = 1, m = 10000, nc = 1, hmax = 30, jmax = 4, verbose = FALSE)

Arguments

L0

Double. Prespecified in-control Average Run Length.

shape1

Double. Shape parameter \(\alpha\) > 0 of the beta-binomial distribution.

shape2

Double. Shape parameter \(\beta\) > 0 of the beta-binomial distribution.

coeff

Numeric Vector. Estimated intercept and slope coefficients from a binary logistic regression model.

RA

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.

rs

Integer. Maximum risk score.

RQ

Double. Defines the performance of a surgeon with the odds ratio ratio of death Q.

m

Integer. Number of simulation runs.

nc

Integer. Number of cores used for parallel processing. Value is passed to parSapply.

hmax

Integer. Maximum value of h for the grid search.

jmax

Integer. Number of digits for grid search.

verbose

Logical. If TRUE verbose output is included, if FALSE a quiet calculation of h is done.

Value

Returns a single value which is the control limit h for a given in-control ARL.

Details

Determines the control limit ("h") for given in-control ARL ("L0") applying a grid search using racusum_betabinomial_arl_sim and parSapply.

Examples

# NOT RUN {
library(vlad)
racusum_betabinomial_crit_sim(L0=100, shape1=1, shape2=3, coeff=c(-3.6798, 0.0768), RA = 2,
rs = 71, RQ = 1, verbose=TRUE)
# }
# NOT RUN {
# }