estimate_pdiff_two function

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