Get one bootstrap computation of the CDF and PMF estimates

one_boot_marg_dist(
  treat,
  covar,
  out,
  treat_form,
  out_levels,
  out_form,
  out_model
)

Arguments

treat

A numeric vector containing treatment status. Should only assume a value 0 or 1.

covar

A data.frame containing the covariates to include in the working proportional odds model.

out

A numeric vector containing the outcomes. Missing outcomes are allowed.

treat_form

The right-hand side of a regression formula for the working model of treatment probability as a function of covariates

out_levels

A numeric vector containing all ordered levels of the outcome.

out_form

The right-hand side of a regression formula for the working proportional odds model. NOTE: THIS FORMULA MUST NOT SUPPRESS THE INTERCEPT.

out_model

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).

Value

Estimates of CDF and PMF for a particular bootstrap sample.