パネル設定でHeckman選択モデルにgllammを使用

Question

パネルデータにHeckman選択モデルを使用します。私はグーグルグーグルとStataのgllammはそれを行うことができるようだ。

しかし、私はそれを使用する方法の適切なチュートリアルを見つけることができませんでした。私はSophia Rabe-Heskethでそれに従うことを試みるが、それらのステップを理解できなかった。

私は現在、断面データに限定しています。 Stataの組み込みのheckmanコマンドと同等でなければなりません。つまり、

use http://www.stata-press.com/data/r13/womenwk 
gen gotowork=1 
replace gotowork=0 if wage==. 
heckman wage educ age, select(gotowork=married children educ age) 

しかし、これらの変数をgllammチュートリアルにマッピングすることは難しいです。具体的には、スライド10では、heckmanコマンドにy1とy2があります。しかし、gllammには1つしかありません。 gotoworkがy1で、wageがy2の場合、このyをどのように定義するのですか？ wageである必要がありますか？

そして私は、

reshape long y, i(id) j(var) 

reshape long wage, i(id) j(var) 

として次のステップを実装しようとするとエラーが

variable var contains all missing values 

なぜ

を言ってしまいましたか？

現在、私は次のステップ

tab gotowork, gen(i) 

をすることによってこの問題を回避し、私は推定ステップ

gen married_i1 = married*i1 
gen children_i1 = children*i1 
gen educ_i1 = educ*i1 
gen age_i1 = age*i1 

gen wage_i2 = wage*i2 
gen educ_i2 = educa*i2 
gen age_i2 = age*i2 

eq load: i1 i2 
constraint define 1 [id1_1]i1 = 1 
gllamm wage married_i1 children_i1 educ_i1 age_i1 i1 wage_i2 educ_i2 age_i2 i2, i(id) eqs(load) nocons constr(1) 

エラーメッセージのために別のエラーを得た：

initial values not feasible 
(error occurred in ML computation) 
(use trace option and check correctness of initial model) 

ことができます誰もこれらのエラーについて説明し、Heckman選択モデルにgllammを正しく使用する方法を教えてください。

私の究極の目標は、パネルのヘッ​​クマン選択モデルを実装することです。これを行うことができる他のSTATA（またはR）パッケージがありますか？

ありがとうございました。

Answer 1

これは、賃金y1とgotowork y2の名前を変更すると機能します。ワイドからロングまでの再フォーマット関数は、新しい変数varを作成します。