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: [email protected]
Data: Gio, 30 Aprile 2009 4:19 pm
A: [email protected]
--------------------------------------------------------------------------
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/
