Statalist


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

st: [Fwd: Heteroskedastic probit]


From   mmolina@uniroma3.it
To   statalist@hsphsun2.harvard.edu
Subject   st: [Fwd: Heteroskedastic probit]
Date   Thu, 30 Apr 2009 17:27:33 +0200 (CEST)

Dear Statilist,
Richard, normscapab is not dichotomy. I constructed it as a score so it
varies from 0 and 1.You are right that then it has small value compared
with the other variables...
I tried Marteen's advice and obtained these results:

hetprob newproc normscapab e6 d12b if a4b==1, het(normscapab)

Fitting probit model:

Iteration 0:   log likelihood = -105.70439
Iteration 1:   log likelihood = -102.70201
Iteration 2:   log likelihood = -102.69757
Iteration 3:   log likelihood = -102.69757

Fitting full model:

Iteration 0:   log likelihood = -102.69757
Iteration 1:   log likelihood = -102.51681
Iteration 2:   log likelihood = -102.41164
Iteration 3:   log likelihood = -102.39453
Iteration 4:   log likelihood = -102.39175
Iteration 5:   log likelihood = -102.39174

Heteroskedastic probit model                    Number of obs     =       
165
                                                Zero outcomes     =       
 56
                                                Nonzero outcomes  =       
109

                                                Wald chi2(3)      =      
1.11
Log likelihood = -102.3917                      Prob > chi2       =    
0.7758

------------------------------------------------------------------------------
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
Interval]
-------------+----------------------------------------------------------------
newproc      |
  normscapab |   3.722963   3.615273     1.03   0.303    -3.362842   
10.80877
          e6 |  -.2859454   .6005927    -0.48   0.634    -1.463085   
.8911945
        d12b |  -.0033544   .0098363    -0.34   0.733    -.0226332   
.0159244
       _cons |  -.1114548    1.10271    -0.10   0.919    -2.272727   
2.049817
-------------+----------------------------------------------------------------
lnsigma2     |
  normscapab |   1.320455   1.518283     0.87   0.384    -1.655325   
4.296234
------------------------------------------------------------------------------
Likelihood-ratio test of lnsigma2=0: chi2(1) =     0.61   Prob > chi2 =
0.4342

The model now converged but the prob is very high. Should I suppose the
presence of Heteroskedasticity now and before?
Best,
Alejandra


-------------------------- Messaggio originale ---------------------------
Oggetto: Heteroskedastic probit
Da:      mmolina@uniroma3.it
Data:    Gio, 30 Aprile 2009 4:19 pm
A:       statalist@hsphsun2.harvard.edu
--------------------------------------------------------------------------

Dear Statalist,
I run this model:
hetprob newprod normscapab n5a e6 d12b Le8c if a4b==5, het(normscapab)
and obtained this result: convergence not achieved.
Is tis anyway to rescue the model?
Best,
Alejandra Molina

Heteroskedastic probit model                    Number of obs     =
165
                                                Zero outcomes     =
 56
                                                Nonzero outcomes  =
109

                                                Wald chi2(5)      =
10.61
Log likelihood = -97.15464                      Prob > chi2       =
0.0596

------------------------------------------------------------------------------
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
Interval]
-------------+----------------------------------------------------------------
newproc      |
  normscapab |  -.1146763   .0552553    -2.08   0.038    -.2229748
-.0063779
         n5a |   4.15e-10   4.03e-10     1.03   0.303    -3.75e-10
1.20e-09
          e6 |   .0098537   .0189325     0.52   0.603    -.0272533
.0469607
        d12b |  -.0004825   .0004277    -1.13   0.259    -.0013208
.0003557
        Le8c |   7.37e-07   4.60e-07     1.60   0.110    -1.66e-07
1.64e-06
       _cons |   .0818565   .0498956     1.64   0.101     -.015937
.1796501
-------------+----------------------------------------------------------------
lnsigma2     |
  normscapab |  -5.064362          .        .       .            .
  .
------------------------------------------------------------------------------
convergence not achieved










*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



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