estimate_gps function

Estimate generalized propensity score (GPS) values