许多文献都在Heckman selection model中加入了固定效应，请问大神们在stata里这是如何操作的呢？是在已有的官方代码（MLE方法或twostep法）中直接加入虚拟变量吗，如果加的话 selection model里需不需要也加入呢？
还是大家都是手动计算的Heckman两步法，那么第一阶段Probit方程是如何加入固定效应的呢？（据我所知Probit方程并不适用固定效应）。如果第一阶段不加固定效应，第二阶段才加，又与文献中模型的设定（第一阶段Probit模型中是有固定效应项的）不符。
不管上述二者中的哪种，总而言之，我是对Probit方程如何加入固定效应感到不解。求助厉害的计量大神们帮助解答，感激不尽！！