.-
help for ^hetprob^                                    (work in progress 7mar199
> 7)
.-

Estimate heteroskedastic probit model
--------------------------------------

        ^hetprob^ depvar [indepvars] [^if^ exp] [^in^ range], variance(varlist)
                                        [ ^nochi l^evel^(^#^)^ maximize_options
>  ]


^hetprob^ shares the features of all estimation commands; see help @weights@.

^hetprob^, typed without arguments, redisplays previous results.

To reset problem-size limits; see help @matsize@.


Description
------------

^hetprob^ estimates heteroskedastic corrected probits.  If you type 

        . ^hetprob y x1 x2 x3, variance(z1 z2)^

^hetprob^ estimates b0, b1, b2, b3, g1, and g2 of

        Pr(y~=0) = F( b0 + b1*x1 + b2*x2 + b3*x3  where 
                                    s^^2 = exp(g1*x1 + g2*z2) )

                 = F( [b0 + b1*x1 + b2*x2 + b3*x3] / 
                                                exp[(g1*x1 + g2*x2)] )


Examples
---------

        . ^hetprob outcome bp age, v(group1 group2)^

        . ^hetprob works educ sex, v(age)^

        . ^hetprob works age, v(age)^


Calculating predicted values
-----------------------------

After ^hetprob^ estimation, you can calculate predicted probabilities of 
a positive outcome by typing 

        . ^predict i1^
        . ^predict s, eq(ln_var)^
        . ^replace s = exp(s/2)^
        . ^gen p = normprob(i1/s)^


Authors
--------

William Gould <wgould@@stata.com>
James Hardin  <jhardin@@stata.com>

This is work in progress.  Results should be viewed with suspicion.  Direct
comments either to the authors or Statalist.  When work on this project is
complete, the authors will seek publication in the Stata Technical Bulletin;
see help @stb@.

Code improvements are invited.  The two ado-files associated with this 
command are

        hetprob.ado 
        hetpr_lf.ado

Optimization is via @ml@ method ^lf^.
        


Also see
---------

On-line:  help for @est@; @probit@, @xtprobit@