fit_missingness_propensity function

Fit missingness propensity model P(delta=1 | X) from pooled data