estimate_pdiff_ind_contrast function

Estimates for a multi-group study with a categorical outcome variable