My gologit2 program, first introduced in 2006, works fine in Stata 11.2 through Stata 14. For example, in Stata 14, if I give the commands
Code:
use https://www3.nd.edu/~rwilliam/statafiles/ordwarm2, clear gologit2 warm i.male i.white, auto waldforce nolabel nofvlabel
Code:
Generalized Ordered Logit Estimates Number of obs = 2,293 Wald chi2(4) = 97.78 Prob > chi2 = 0.0000 Log likelihood = -2944.6035 Pseudo R2 = 0.0171 ( 1) [eq1]1.white - [eq2]1.white = 0 ( 2) [eq2]1.white - [eq3]1.white = 0 ------------------------------------------------------------------------------ warm | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- eq1 | 1.male | -.3042833 .1247377 -2.44 0.015 -.5487647 -.0598019 1.white | -.3331343 .1159147 -2.87 0.004 -.5603229 -.1059457 _cons | 2.351311 .1366675 17.20 0.000 2.083447 2.619174 -------------+---------------------------------------------------------------- eq2 | 1.male | -.5922755 .0851572 -6.96 0.000 -.7591805 -.4253705 1.white | -.3331343 .1159147 -2.87 0.004 -.5603229 -.1059457 _cons | .7943304 .117114 6.78 0.000 .5647912 1.02387 -------------+---------------------------------------------------------------- eq3 | 1.male | -1.01099 .1195953 -8.45 0.000 -1.245393 -.7765879 1.white | -.3331343 .1159147 -2.87 0.004 -.5603229 -.1059457 _cons | -.8269653 .1184391 -6.98 0.000 -1.059102 -.5948289 ------------------------------------------------------------------------------
Code:
Generalized Ordered Logit Estimates Number of obs = 2,293 Wald chi2(4) = 97.78 Prob > chi2 = 0.0000 Log likelihood = -2944.6035 Pseudo R2 = 0.0171 ( 1) [eq1]1.white - [eq2]1.white = 0 ( 2) [eq2]1.white - [eq3]1.white = 0 ------------------------------------------------------------------------------ warm | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- eq1 | 1.male | -.3042833 .1247377 -2.44 0.015 -.5487647 -.0598019 1.white | -.3331343 .1159147 -2.87 0.004 -.5603229 -.1059457 _cons | 2.351311 .1366675 17.20 0.000 2.083447 2.619174 -------------+---------------------------------------------------------------- eq2 | 0.male | -3.38e-16 3.76e-17 -8.99 0.000 -4.12e-16 -2.64e-16 1.male | -.5922755 .0851572 -6.96 0.000 -.7591805 -.4253705 0.white | 9.66e-17 1.72e-17 5.62 0.000 6.29e-17 1.30e-16 1.white | -.3331343 .1159147 -2.87 0.004 -.5603229 -.1059457 _cons | .7943304 .117114 6.78 0.000 .5647912 1.02387 -------------+---------------------------------------------------------------- eq3 | 1.male | -1.01099 .1195953 -8.45 0.000 -1.245393 -.7765879 1.white | -.3331343 .1159147 -2.87 0.004 -.5603229 -.1059457 _cons | -.8269653 .1184391 -6.98 0.000 -1.059102 -.5948289 ------------------------------------------------------------------------------
gologit2 calls the ml command. My guess is that something changed in ml in Stata 15 that causes it to add these junk parameters. That may need fixes on Stata's end, but I am happy to program around it on my side if I can figure out how.
Also, you can bypass gologit2 and call ml directly, and the problem is the same. I've tried tweaking the ml command (e.g. drop the collinear option) but so far no luck. Dropping the 0b. variables doesn't work if there are more than 2 categories to a variable. If there is some option I can add to ml that will make the problem go away I would love to know what it is! Thanks
Code:
use https://www3.nd.edu/~rwilliam/statafiles/ordwarm2, clear constraint 1998 [eq1]1.white = [eq2]1.white constraint 1999 [eq1]1.white = [eq3]1.white ml model lf gologit2_ll /// (eq1:warm= 0b.male 1.male 0b.white 1.white) (eq2: 0b.male 1.male 0b.white 1.white) /// (eq3: 0b.male 1.male 0b.white 1.white) , /// constraints(1999 1998) waldtest(-3) /// lf0(3 -2995.770427280555) collinear /// title(Generalized Ordered Logit Estimates) missing maximize nocnsnotes nolog collinear ml display, nofvlabel
0 Response to gologit2 suddenly having problems in Stata 15
Post a Comment