Description Usage Arguments Value Author(s) Examples
View source: R/racusum_betabinomial_sim.R
Compute ARLs of RA-CUSUM control charts using simulation.
1 2 | racusum_betabinomial_arl_sim(r, shape1, shape2, coeff, h, RA = 2, rs = 71,
RQ = 1)
|
r |
Integer Vector. Number of runs. |
shape1 |
Double. Shape parameter alpha |
shape2 |
Double. Shape parameter beta |
coeff |
Numeric Vector. Estimated intercept and slope coefficients from a binary logistic regression model. |
h |
Double. Control Chart limit for detecting deterioration/improvement. |
RA |
Double. Odds ratio of death under the alternative hypotheses. Detecting deterioration
in performance with increased mortality risk by doubling the odds Ratio |
rs |
Integer. Maximum risk score. |
RQ |
Double. Defines the performance of a surgeon with the odds ratio ratio of death
|
Returns a single value which is the Run Length.
Philipp Wittenberg
1 2 3 4 5 6 7 8 | ## Not run:
library(vlad)
m <- 1e3
RLS <- sapply(1:m, racusum_betabinomial_arl_sim, shape1=1, shape2=3, coeff=c(-3.6798, 0.0768),
h=4.5, RA=2, rs=71, RQ=1)
data.frame(cbind(ARL=mean(RLS), ARLSE=sd(RLS)/sqrt(m)))
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.