
Compute a BCa bootstrap confidence interval for the Mann-Whitney parameter. The code is based on the slides found here: http://users.stat.umn.edu/~helwig/notes/bootci-Notes.pdf
bca_mannwhitney(
treat,
covar,
out,
nboot,
treat_form,
out_levels,
out_form,
mannwhitney_est,
out_model,
alpha = 0.05
)
A numeric
vector containing treatment status. Missing
values are not allowed unless the corresponding entry in out
is also missing.
Only values of 0 or 1 are treated as actual treatment levels. Any other value is assumed
to encode a value for which the outcome is missing and the corresponding outcome value is
ignored.
A data.frame
containing the covariates to include in the working
proportional odds model.
A numeric
vector containing the outcomes. Missing outcomes are
allowed.
Number of bootstrap replicates used to compute bootstrap confidence intervals.
The right-hand side of a regression formula for the working model of treatment probability as a function of covariates
A numeric
vector containing all ordered levels of the
outcome.
The right-hand side of a regression formula for the working proportional odds model. NOTE: THIS FORMULA MUST NOT SUPPRESS THE INTERCEPT.
The point estimate of the Mann-Whitney parameter.
Which R function should be used to fit the proportional odds
model. Options are "polr"
(from the MASS
package),
"vglm" (from the VGAM
package), or "clm"
(from the ordinal
package).
Level of confidence interval.
Confidence interval for the Mann-Whitney parameter