使用stata估计Heckman-style selection models很方便
即使考虑时间效应（年份效应）操作起来也很方便，不需要做很大的改动。


但是选择R软件来做
认真看了官网下载的package说明（Package ‘sampleSelection’  ）和一篇关于具体操作的文章“Sample Selection Models in R:Package~sampleSelection”，如附件所示。
还是没找到直接简单的方法。


难道非要手动生成 时间哑变量，作为具体的控制变量，考虑进每个阶段的回归模型中去吗？


请教一下各位大神，非常感谢^_^