Calculates AU p-values for testing independence in 2x2 case-control tables, while adjusting for categorical covariates.
Inputs are given as a vector of counts in each strata defined by the covariate(s). Note that computational time can be extremely high.