Hi all,

I am trying to compute inverse probability weight for multiple treatments. But I am struggling to do so.

The common way to generate ipw for binary treatments is as follows.


* Estimate the propensity score

qui logit treatment x1 x2 x3, nolog

predict double pscore if e(sample)

* Create IP weight

gen ipw = 1/pscore if treatment == 1

replace ipw = 1/(1-pscore) if treatment ==0


But I’m trying to compute inverse probability weighting for multiple treatment groups. My treatment values are 0, 1, 2, 3.

I have obtained propensity score (pscore) for each group. Does anyone know how to compute “ipw” for each treatment group, please?

Thank you so much for your time.