Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down on April 23, and its replacement, is already up and running.

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

st: Consistent estimation of the Double Hurdle model

From   verdes <>
Subject   st: Consistent estimation of the Double Hurdle model
Date   Wed, 8 Feb 2012 02:19:34 -0800 (PST)

Dear Statalisters, 

I've been trying to estimate a Double Hurdle model on alcohol expenditure

-craggit- by William Burke allows for hetereschedasticity but not for
non-normality. So, I guess estimates would still be subject to bias.

I then came across -dh- written by Moffatt (2005) in the Journal of the
Operational Research Society, Vol.56, No.9 (*see below the command). -dh-
allows for non-normality by applying a Box-Cox transformation to the
dependant variable.  

I'd like to ask you the following question:
I was told that, to deal with non-normality, the Inverse Hyperbolic Sine
transformation is a better solution than the Box-Cox transformation. Is that
true? If so, why?
Also, Moffatt's Stata code implies homoscedasticity (sigma is an invariant
parameter). Will the estimates be consistent still?

I'd greatly appreciate your help.

Stefano Verde

program define dh
version 6
args Inf thetal theta2 theta3 theta4
tempvar d p z p0 pll yt
quietly gen double 'd'= $ML_yl >0
quietly gen double 'p'= normprob('theta3')
quietly gen double '1'= 'theta4'
quietly gen double 'yt'= ($MLyl^'I'-1)/'l'
quietly gen double 'z'= ('yt'-'thetal' )/('theta2')
quietly gen double 'pO'= 1-('p'*normprob(-'z'))
quietly gen double 'pl' = (($ML_yl + (1-'d'))^
quietly replace 'lnf = ln((1-'d')*'p0+' 'd'*'pl')
ml model If dh (y = 'listy') () (d - 'listd') ()
ml init b, copy
ml maximize

Notes: 'listy' is a previously defined list of variables
appearing in the second hurdle; 'listd' contains the variables
of the first hurdle. 'thetal' corresponds to xi'f# in (14),
'theta2' to a, 'theta3' to zi'a, and 'theta4' to .. b is a vector of
suitable starting values.

View this message in context:
Sent from the Statalist mailing list archive at
*   For searches and help try:

© Copyright 1996–2015 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index