Hi,

I am a first-time user of the generalized propensity score. I have 3 treatment levels that are qualitatively different. I used mlogit and predict commands to estimate three sets of gpscore. But I am not sure how to match the pairs of treatments. I have read in several papers that the researchers used nearest neighbour matching to match the treated and control units. Can I use the same if I have a vector of propensity scores, and can I do it in Stata? If not, is there any other way to match the units based on the gpscores in Stata?

Will appreciate any help on this matter.

Thanks,
Nadia