前一种，目前不是采用配对方法解决的，而是基于 treatment effect 模型的思路解决。
findit mtreatreg

参见：
Deb, P., P. Trivedi, 2006, Maximum simulated likelihood estimation of a negative binomial regression model with multinomial endogenous treatment, Stata Journal, 6 (2), pp. 246-255.

后一种，就转化成普通的内生变量了，采用 ivreg 或 ivreg2 即可。

got it . thanks