I need some help in estimating the inverse Mills ration for use in a
Heckman 2-stage model using SAS.
I have a repeated neasures dataset and use the proc GenMod to create
the probabiilties for an event to occur (binary variable) using a
tobit link function. Once the probabilities for each observation are
obtained, I calculate the =93inverse Mills ratio=94 which will then used
in the 2nd stage as an additional explanatory variale.
The question is how to calculate the inverse mills ratio (IMR) which
is defined as the ratio between the Probability density function
(PDF) over the cumulative density function (CDF) associated with the
Does the SAS code below create the IMR properly?
/* First Stage PROBIT Model With Panel Data & Firm Fixed effects*/
PROC GENMOD DATA=3D myData DESCENDING;
CLASS FID Year;
Model BinVar=3D exVar1 exVar2
/ dist =3D Binomial link=3Dprobit MAXITER =3D 1000 ;
REPEATED subject=3DFID / WITHIN=3DYear;
OUTPUT OUT=3DmyData PROB=3DPbinVar;
/*calculate inverse mills ratio */
data myData ;
set myData ;
IMR =3D pdf('NORMAL',PbinVar ) / cdf('NORMAL',PbinVar );
Thanks in advance!
||7/19/2010 8:40:45 AM