[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

RE: st: Censored ... Heteroskedastic regression

From   Fabrice <>
Subject   RE: st: Censored ... Heteroskedastic regression
Date   Fri, 21 Dec 2007 15:39:39 +0100


Thanks a lot for very practical contribution. I have to think about it and
test it (I'm still not fully at ease with those likelihood estimation

Great, really! :-)

Fabrice Cavarretta
Mobile: +33 6 09 59 46 74

-----Original Message-----
[] On Behalf Of Maarten buis
Sent: Friday, December 21, 2007 12:31 PM
Subject: Re: st: Censored ... Heteroskedastic regression

--- Fabrice <> wrote:
> I am trying to estimate an effect on variability of a dependent
> variable that happens to be censored.
> Normally, without censoring, I use a simple Maximum Likelihood
> Estimation, modeling in Stata the residual as I see fit. For
> instance, I can parameterize a simple linear effect both on mean and
> variance by:
> ml model lf mynormal1_lf (mu: Y=X) (sigma: X)
> For censoring, if I were to just seek the main effect, I would use a
> normal tobit.
> QUESTION: would anyone have a suggestion how to estimate (in Stata!)
> effect on dependent variable variance in the case of censored DV?
> Obviously, I could write my own likelihood function, but do not feel
> confident enough to do so yet (and maybe that is already done
> somewhere).

Below is a quick stab at a heteroskedastic tobit: the global tau is the
value at which the dependent variable is censored.

*-------------- begin example -------------------
set more off
sysuse nlsw88, clear
gen lnw = ln(wage)

capture program drop hettobit_lf
program define hettobit_lf
	version 9.2
	args lnf mu lnsigma
	quietly {
		replace `lnf' =                         ///
                lnnormalden($ML_y1,`mu',exp(`lnsigma')) ///
                if $ML_y1 > $tau

		replace `lnf' =                         /// 
                lnnormal(($tau -`mu')/exp(`lnsigma'))   ///
		if $ML_y <= $tau

global tau = 1.4

ml model lf hettobit_lf                                 ///
         (mu: lnw=grade ttl_exp union)                  ///
         (lnsigma: union)
ml check
ml search
ml maximize
*--------------- end example -------------------------
(For more on how to use examples I sent to the Statalist, see )

Hope this helps,

Maarten L. Buis
Department of Social Research Methodology
Vrije Universiteit Amsterdam
Boelelaan 1081
1081 HV Amsterdam
The Netherlands

visiting address:
Buitenveldertselaan 3 (Metropolitan), room Z434

+31 20 5986715

Sent from Yahoo! Mail - a smarter inbox

*   For searches and help try:

*   For searches and help try:

© Copyright 1996–2021 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index