. * prchange: version 1.6.5 . * verified: 3/25/01 . which prchange c:\stata6\ado\stbplus\p\prchange.ado *! version 1.6.5 1/10/01 uncen option changed to uncentered . which _pemarg c:\stata6\ado\stbplus\_\_pemarg.ado *! version 1.0.2 2/12/01 . . use gsstestdata.dta, clear . * logit and logistic . qui logit feminist age educ lesbian liberal . prchange educ, x(age 40 educ 12) fromto h logit: Changes in Predicted Probabilities for feminist from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x-1/2 educ 0.1368 0.3659 0.2291 0.1368 0.1446 0.0078 0.2499 to: dif: from: to: dif: x+1/2 -+1/2 x-1/2sd x+1/2sd -+sd/2 MargEfct educ 0.2622 0.0123 0.2394 0.2734 0.0340 0.0123 No Yes Pr(y|x) 0.7440 0.2560 age educ lesbian liberal x= 40 12 .024523 2.90191 sd(x)= 16.8692 2.75988 .154772 1.35255 Pr(y|x): probability of observing each y for specified x values Avg|Chg|: average of absolute value of the change across categories Min->Max: change in predicted probability as x changes from its minimum to its maximum 0->1: change in predicted probability as x changes from 0 to 1 -+1/2: change in predicted probability as x changes from 1/2 unit below base value to 1/2 unit above -+sd/2: change in predicted probability as x changes from 1/2 standard dev below base to 1/2 standard dev above MargEfct: the partial derivative of the predicted probability/rate with respect to a given independent variable . prchange educ, x(age 40 educ 12) fromto uncentered h logit: Changes in Predicted Probabilities for feminist from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x educ 0.1368 0.3659 0.2291 0.1368 0.1446 0.0078 0.2560 to: dif: from: to: dif: x+1 +1 x x+sd +sd MargEfct educ 0.2685 0.0125 0.2560 0.2914 0.0354 0.0123 No Yes Pr(y|x) 0.7440 0.2560 age educ lesbian liberal x= 40 12 .024523 2.90191 sd(x)= 16.8692 2.75988 .154772 1.35255 Pr(y|x): probability of observing each y for specified x values Avg|Chg|: average of absolute value of the change across categories Min->Max: change in predicted probability as x changes from its minimum to its maximum 0->1: change in predicted probability as x changes from 0 to 1 +1: change in predicted probability as x increases by 1 unit +sd: change in predicted probability as x increases by a standard dev MargEfct: the partial derivative of the predicted probability/rate with respect to a given independent variable . prchange educ, x(age 40 educ 12) fromto d(4) h logit: Changes in Predicted Probabilities for feminist (Note: d = 4) from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x-d/2 educ 0.1368 0.3659 0.2291 0.1368 0.1446 0.0078 0.2322 to: dif: from: to: dif: x+d/2 -+d/2 x-1/2sd x+1/2sd -+sd/2 MargEfct educ 0.2814 0.0492 0.2394 0.2734 0.0340 0.0123 No Yes Pr(y|x) 0.7440 0.2560 age educ lesbian liberal x= 40 12 .024523 2.90191 sd(x)= 16.8692 2.75988 .154772 1.35255 Pr(y|x): probability of observing each y for specified x values Avg|Chg|: average of absolute value of the change across categories Min->Max: change in predicted probability as x changes from its minimum to its maximum 0->1: change in predicted probability as x changes from 0 to 1 -+d/2: change in predicted probability as x changes from 2 units below base value to 2 units above -+sd/2: change in predicted probability as x changes from 1/2 standard dev below base to 1/2 standard dev above MargEfct: the partial derivative of the predicted probability/rate with respect to a given independent variable . prchange educ, x(age 40 educ 12) fromto uncentered d(4) h logit: Changes in Predicted Probabilities for feminist (Note: delta = 4) from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x educ 0.1368 0.3659 0.2291 0.1368 0.1446 0.0078 0.2560 to: dif: from: to: dif: x+delta +delta x-sd x+sd +sd MargEfct educ 0.3082 0.0522 0.2560 0.2914 0.0354 0.0123 No Yes Pr(y|x) 0.7440 0.2560 age educ lesbian liberal x= 40 12 .024523 2.90191 sd(x)= 16.8692 2.75988 .154772 1.35255 Pr(y|x): probability of observing each y for specified x values Avg|Chg|: average of absolute value of the change across categories Min->Max: change in predicted probability as x changes from its minimum to its maximum 0->1: change in predicted probability as x changes from 0 to 1 +delta: change in predicted probability as x increases by 4 units +sd: change in predicted probability as x increases by a standard dev MargEfct: the partial derivative of the predicted probability/rate with respect to a given independent variable . prchange educ lesbian liberal, x(age 40 educ 12) logit: Changes in Predicted Probabilities for feminist min->max 0->1 -+1/2 -+sd/2 MargEfct educ 0.2291 0.0078 0.0123 0.0340 0.0123 lesbian 0.0074 0.0074 0.0073 0.0011 0.0073 liberal 0.3234 0.0359 0.0543 0.0734 0.0543 No Yes Pr(y|x) 0.7440 0.2560 age educ lesbian liberal x= 40 12 .024523 2.90191 sd(x)= 16.8692 2.75988 .154772 1.35255 . qui logistic feminist age educ lesbian liberal . prchange educ lesbian liberal, x(age 40 educ 12) logistic: Changes in Predicted Probabilities for feminist min->max 0->1 -+1/2 -+sd/2 MargEfct educ 0.2291 0.0078 0.0123 0.0340 0.0123 lesbian 0.0074 0.0074 0.0073 0.0011 0.0073 liberal 0.3234 0.0359 0.0543 0.0734 0.0543 No Yes Pr(y|x) 0.7440 0.2560 age educ lesbian liberal x= 40 12 .024523 2.90191 sd(x)= 16.8692 2.75988 .154772 1.35255 . . * probit . qui probit feminist age educ lesbian liberal . prchange, x(age 40 educ max) rest(median) probit: Changes in Predicted Probabilities for feminist min->max 0->1 -+1/2 -+sd/2 MargEfct age 0.1512 0.0019 0.0021 0.0353 0.0021 educ 0.2248 0.0084 0.0137 0.0379 0.0137 lesbian 0.0079 0.0079 0.0078 0.0012 0.0078 liberal 0.3627 0.0499 0.0626 0.0846 0.0627 No Yes Pr(y|x) 0.6320 0.3680 age educ lesbian liberal x= 40 20 0 3 sd(x)= 16.8692 2.75988 .154772 1.35255 . . * cloglog . qui cloglog feminist age educ lesbian liberal . prchange, rest(upper) cloglog: Changes in Predicted Probabilities for feminist min->max 0->1 -+1/2 -+sd/2 age 0.2026 0.0028 0.0026 0.0447 educ 0.4060 0.0174 0.0196 0.0539 lesbian 0.0090 0.0090 0.0089 0.0014 liberal 0.4865 0.0650 0.0817 0.1102 No Yes Pr(y|x) 0.2209 0.7791 age educ lesbian liberal x= 89 20 1 6 sd(x)= 16.8692 2.75988 .154772 1.35255 . . * ologit . qui ologit mchilds age educ lesbian liberal . prchange, x(age max educ min) ologit: Changes in Predicted Probabilities for mchilds age Avg|Chg| no_child one_chil two_chil three_ch Min->Max .2496515 -.10344535 -.12676268 -.26909496 .49930301 -+1/2 .00123718 -.00026703 -.00042816 -.00177916 .00247437 -+sd/2 .02159167 -.00469205 -.00750726 -.03098403 .04318333 MargEfct .00494847 -.00026701 -.00042814 -.00177909 .00247424 educ Avg|Chg| no_child one_chil two_chil three_ch Min->Max .22322204 .08428647 .1084691 .25368847 -.44644409 -+1/2 .00407586 .00087992 .00141079 .00586101 -.00815171 -+sd/2 .01155446 .00249877 .00400414 .01660602 -.0231089 MargEfct .01629535 .00087926 .00140985 .00585856 -.00814767 lesbian Avg|Chg| no_child one_chil two_chil three_ch 0->1 .02885085 .00657841 .01037129 .04075199 -.05770171 liberal Avg|Chg| no_child one_chil two_chil three_ch Min->Max .01984755 .00431491 .00690305 .02847714 -.03969508 -+1/2 .00323465 .00069825 .00111955 .00465151 -.00646931 -+sd/2 .00435272 .00093972 .00150666 .00625907 -.00870544 MargEfct .01293457 .00069792 .00111908 .00465028 -.00646729 no_child one_chil two_chil three_ch Pr(y|x) .00649185 .01059144 .04675842 .9361583 age educ lesbian liberal x= 89 0 .024503 2.87682 sd(x)= 17.2202 2.82507 .154657 1.34531 . qui ologit mchilds age educ lesbian liberal faminc . prchange, x(age 40 educ min) fromto outcome(1) d(10) ologit: Changes in Predicted Probabilities for mchilds Outcome: 1 (one_child) (Note: d = 10) from: to: dif: from: to: x=min x=max min->max x=0 x=1 age .06609333 .0036208 -.06247252 .13061032 .12593468 educ .02742963 .48687944 .45944981 .02742963 .03252982 lesbian .02694573 .05401887 .02707315 .02694573 .05401887 liberal .02059999 .03728127 .01668127 .02059999 .02275573 faminc .03878095 .01574046 -.02304049 .03897412 .0385887 dif: from: to: dif: from: 0->1 x-d/2 x+d/2 -+d/2 x-1/2sd age -.00467564 .03359435 .02236997 -.01122439 .03859771 educ .00510019 .01157489 .06360397 .05202908 .02155268 lesbian .02707315 .00075582 .51257533 .51181951 .02595374 liberal .00215574 .016676 .04480184 .02812584 .02566357 faminc -.00038542 .02884298 .02608367 -.00275932 .03113415 to: dif: x+1/2sd -+sd/2 MargEfct age .01942768 -.01917003 -.00168521 educ .03485201 .01329933 .00708404 lesbian .02898694 .00303319 .02916475 liberal .02931356 .00364999 .00409951 faminc .0241549 -.00697925 -.00041666 no_child one_chil two_chil three_ch Pr(y|x) .02742963 .0447565 .17205495 .75575894 age educ lesbian liberal faminc x= 40 0 .025281 2.88343 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . prchange, x(age 40 educ min) fromto uncentered outcome(1) d(10) ologit: Changes in Predicted Probabilities for mchilds Outcome: 1 (one_child) (Note: delta = 10) from: to: dif: from: to: x=min x=max min->max x=0 x=1 age .06609333 .0036208 -.06247252 .13061032 .12593468 educ .02742963 .48687944 .45944981 .02742963 .03252982 lesbian .02694573 .05401887 .02707315 .02694573 .05401887 liberal .02059999 .03728127 .01668127 .02059999 .02275573 faminc .03878095 .01574046 -.02304049 .03897412 .0385887 dif: from: to: dif: from: 0->1 x x+delta +delta x age -.00467564 .02742963 .01822612 -.0092035 .02742963 educ .00510019 .02742963 .140589 .11315937 .02742963 lesbian .02707315 .02742963 .9751308 .94770117 .02742963 liberal .00215574 .02742963 .07235795 .04492833 .02742963 faminc -.00038542 .02742963 .02480207 -.00262756 .02742963 to: dif: x+sd +sd MargEfct age .01372716 -.01370246 -.00168521 educ .04419161 .01676198 .00708404 lesbian .03062988 .00320025 .02916475 liberal .03132273 .0038931 .00409951 faminc .02126258 -.00616704 -.00041666 no_child one_chil two_chil three_ch Pr(y|x) .02742963 .0447565 .17205495 .75575894 age educ lesbian liberal faminc x= 40 0 .025281 2.88343 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . *check unusual values . qui ologit dmchilds age educ lesbian liberal faminc . prchange ologit: Changes in Predicted Probabilities for dmchilds age Avg|Chg| 0 1_5 2_200000 3_5 Min->Max .2945401 -.39238448 -.18630595 -.01038978 .58908018 -+1/2 .00504413 -.00666231 -.00342596 .00203991 .00804836 -+sd/2 .08440455 -.11260611 -.05620301 .03332883 .13548027 MargEfct .02017718 -.00666229 -.0034263 .00204014 .00804845 educ Avg|Chg| 0 1_5 2_200000 3_5 Min->Max .34451257 .41687278 .20667724 .06547514 -.68902514 -+1/2 .02119228 .02800758 .01437698 -.00855836 -.03382619 -+sd/2 .05935335 .07878189 .0399248 -.02372485 -.09498185 MargEfct .08481797 .028006 .01440299 -.00857606 -.03383293 lesbian Avg|Chg| 0 1_5 2_200000 3_5 0->1 .08946461 .13851535 .04041387 -.06321964 -.11570957 liberal Avg|Chg| 0 1_5 2_200000 3_5 Min->Max .07327475 .09798834 .04856116 -.0300166 -.11653289 -+1/2 .0122686 .01620729 .00832991 -.00495949 -.01957771 -+sd/2 .01648718 .02178392 .01119046 -.00666216 -.0263122 MargEfct .04908389 .01620698 .00833496 -.00496294 -.01957901 faminc Avg|Chg| 0 1_5 2_200000 3_5 Min->Max .10757682 -.13865639 -.07649724 .02970704 .18544659 -+1/2 .00124717 -.00164722 -.00084713 .0005044 .00198993 -+sd/2 .03144219 -.04158539 -.02129899 .01267505 .05020933 MargEfct .00498872 -.00164722 -.00084714 .00050442 .00198994 0 1_5 2_200000 3_5 Pr(y|x) .19885878 .20758004 .3334344 .26012677 age educ lesbian liberal faminc x= 44.5442 13.4529 .025281 2.88343 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . qui ologit ngchilds age educ lesbian liberal faminc . prchange ologit: Changes in Predicted Probabilities for ngchilds age Avg|Chg| -17 -4 12 78 Min->Max .2945401 -.39238448 -.18630595 -.01038978 .58908018 -+1/2 .00504413 -.00666231 -.00342596 .00203991 .00804836 -+sd/2 .08440455 -.11260611 -.05620301 .03332883 .13548027 MargEfct .02017718 -.00666229 -.0034263 .00204014 .00804845 educ Avg|Chg| -17 -4 12 78 Min->Max .34451257 .41687278 .20667724 .06547514 -.68902514 -+1/2 .02119228 .02800758 .01437698 -.00855836 -.03382619 -+sd/2 .05935335 .07878189 .0399248 -.02372485 -.09498185 MargEfct .08481797 .028006 .01440299 -.00857606 -.03383293 lesbian Avg|Chg| -17 -4 12 78 0->1 .08946461 .13851535 .04041387 -.06321964 -.11570957 liberal Avg|Chg| -17 -4 12 78 Min->Max .07327475 .09798834 .04856116 -.0300166 -.11653289 -+1/2 .0122686 .01620729 .00832991 -.00495949 -.01957771 -+sd/2 .01648718 .02178392 .01119046 -.00666216 -.0263122 MargEfct .04908389 .01620698 .00833496 -.00496294 -.01957901 faminc Avg|Chg| -17 -4 12 78 Min->Max .10757682 -.13865639 -.07649724 .02970704 .18544659 -+1/2 .00124717 -.00164722 -.00084713 .0005044 .00198993 -+sd/2 .03144219 -.04158539 -.02129899 .01267505 .05020933 MargEfct .00498872 -.00164722 -.00084714 .00050442 .00198994 -17 -4 12 78 Pr(y|x) .19885878 .20758004 .3334344 .26012677 age educ lesbian liberal faminc x= 44.5442 13.4529 .025281 2.88343 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . qui ologit irchilds age educ lesbian liberal faminc . prchange ologit: Changes in Predicted Probabilities for irchilds age Avg|Chg| 0 3 9 37 Min->Max .2945401 -.39238448 -.18630595 -.01038978 .58908018 -+1/2 .00504413 -.00666231 -.00342596 .00203991 .00804836 -+sd/2 .08440455 -.11260611 -.05620301 .03332883 .13548027 MargEfct .02017718 -.00666229 -.0034263 .00204014 .00804845 educ Avg|Chg| 0 3 9 37 Min->Max .34451257 .41687278 .20667724 .06547514 -.68902514 -+1/2 .02119228 .02800758 .01437698 -.00855836 -.03382619 -+sd/2 .05935335 .07878189 .0399248 -.02372485 -.09498185 MargEfct .08481797 .028006 .01440299 -.00857606 -.03383293 lesbian Avg|Chg| 0 3 9 37 0->1 .08946461 .13851535 .04041387 -.06321964 -.11570957 liberal Avg|Chg| 0 3 9 37 Min->Max .07327475 .09798834 .04856116 -.0300166 -.11653289 -+1/2 .0122686 .01620729 .00832991 -.00495949 -.01957771 -+sd/2 .01648718 .02178392 .01119046 -.00666216 -.0263122 MargEfct .04908389 .01620698 .00833496 -.00496294 -.01957901 faminc Avg|Chg| 0 3 9 37 Min->Max .10757682 -.13865639 -.07649724 .02970704 .18544659 -+1/2 .00124717 -.00164722 -.00084713 .0005044 .00198993 -+sd/2 .03144219 -.04158539 -.02129899 .01267505 .05020933 MargEfct .00498872 -.00164722 -.00084714 .00050442 .00198994 0 3 9 37 Pr(y|x) .19885878 .20758004 .3334344 .26012677 age educ lesbian liberal faminc x= 44.5442 13.4529 .025281 2.88343 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . . * oprobit . qui oprobit mchilds age educ lesbian liberal . prchange, x(age mean educ 12) rest(upper) oprobit: Changes in Predicted Probabilities for mchilds age Avg|Chg| no_child one_chil two_chil three_ch Min->Max .29206494 -.29268296 -.17956876 -.11187817 .58412988 -+1/2 .0044135 -.00490481 -.00302945 -.00089276 .008827 -+sd/2 .07554691 -.08459712 -.05140275 -.01509395 .15109381 MargEfct .01765436 -.00490477 -.00302959 -.00089282 .00882718 educ Avg|Chg| no_child one_chil two_chil three_ch Min->Max .2902947 .30147816 .17787656 .1012347 -.58058938 -+1/2 .01559585 .01733709 .01070139 .0031532 -.03119171 -+sd/2 .04398145 .04900215 .03010079 .00885996 -.0879629 MargEfct .06239925 .0173359 .01070807 .00315566 -.03119962 lesbian Avg|Chg| no_child one_chil two_chil three_ch 0->1 .0736718 .10258898 .04475462 -.00767118 -.13967243 liberal Avg|Chg| no_child one_chil two_chil three_ch Min->Max .07300512 .10155417 .04445608 -.00741476 -.13859549 -+1/2 .01248695 .01387948 .0085693 .00252512 -.0249739 -+sd/2 .01679664 .01867293 .01152465 .00339571 -.0335933 MargEfct .04995594 .01387887 .00857273 .00252638 -.02497797 no_child one_chil two_chil three_ch Pr(y|x) .13154861 .1634343 .31440631 .39061078 age educ lesbian liberal x= 45.0079 12 0 0 sd(x)= 17.2202 2.82507 .154657 1.34531 . qui oprobit mchilds age educ lesbian liberal faminc . prchange, x(age mean educ 12) rest(upper) outcome(0) fromto oprobit: Changes in Predicted Probabilities for mchilds Outcome: 0 (no_children) from: to: dif: from: to: x=min x=max min->max x=0 x=1 age 0 0 .28980349 0 0 educ 0 0 .34210843 0 0 lesbian 0 0 .07553129 0 0 liberal 0 0 .07190567 0 0 faminc 0 0 .10665264 0 0 dif: from: to: dif: from: 0->1 x-1/2 x+1/2 -+1/2 x-1/2sd age .00462893 0 0 .00466971 0 educ .00964441 0 0 .02061324 0 lesbian .07553129 0 0 .0762897 0 liberal .01216429 0 0 .0121537 0 faminc .0011716 0 0 .00123858 0 to: dif: x+1/2sd -+sd/2 MargEfct age 0 .0783537 -.00306768 educ 0 .05778337 .01354724 lesbian 0 .01205207 .05042547 liberal 0 .01633329 .00798522 faminc 0 .03123318 -.00081365 no_child one_chil two_chil three_ch Pr(y|x) .06773811 .11829851 .29650477 .51745862 age educ lesbian liberal faminc x= 44.5442 12 0 0 90 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . prchange age educ, x(age mean educ 12) rest(upper) d(10) oprobit: Changes in Predicted Probabilities for mchilds (Note: d = 10) age Avg|Chg| no_child one_chil two_chil three_ch Min->Max .28980349 -.1862952 -.18052895 -.21278282 .57960698 -+d/2 .04659174 -.03076282 -.0319759 -.03044477 .09318346 -+sd/2 .0783537 -.05221919 -.05364198 -.05084625 .15670738 MargEfct .01867926 -.00306768 -.00320912 -.00306283 .00933963 educ Avg|Chg| no_child one_chil two_chil three_ch Min->Max .34210843 .2499115 .20788272 .22642265 -.68421687 -+d/2 .19739672 .142473 .13202588 .12029456 -.39479345 -+sd/2 .05778337 .03825696 .03962828 .03768149 -.11556673 MargEfct .08248973 .01354724 .01417182 .0135258 -.04124487 no_child one_chil two_chil three_ch Pr(y|x) .06773811 .11829851 .29650477 .51745862 age educ lesbian liberal faminc x= 44.5442 12 0 0 90 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . prchange age educ, x(age mean educ 12) rest(upper) uncentered d(10) oprobit: Changes in Predicted Probabilities for mchilds (Note: delta = 10) age Avg|Chg| no_child one_chil two_chil three_ch Min->Max .28980349 -.1862952 -.18052895 -.21278282 .57960698 +delta .04603881 -.0256709 -.03047986 -.03592685 .09207761 +sd .07619995 -.03826519 -.0489581 -.06517659 .1523999 MargEfct .01867926 -.00306768 -.00320912 -.00306283 .00933963 educ Avg|Chg| no_child one_chil two_chil three_ch Min->Max .34210843 .2499115 .20788272 .22642265 -.68421687 +delta .18526104 .25573654 .11478554 -.0138891 -.35663299 +sd .05754206 .04696906 .04070968 .02740538 -.11508411 MargEfct .08248973 .01354724 .01417182 .0135258 -.04124487 no_child one_chil two_chil three_ch Pr(y|x) .06773811 .11829851 .29650477 .51745862 age educ lesbian liberal faminc x= 44.5442 12 0 0 90 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . . * mlogit . qui mlogit mchilds age educ lesbian liberal . prchange, x(age 40 educ 12) mlogit: Changes in Predicted Probabilities for mchilds age Avg|Chg| no_child one_chil two_chil three_ch Min->Max .28475921 -.36528246 -.20423594 .05252677 .51699166 -+1/2 .00528009 -.00723974 -.00332044 .0025537 .00800648 -+sd/2 .08999199 -.12401797 -.055966 .04339501 .13658898 MargEfct .02112111 -.00723987 -.00332069 .00255381 .00800675 educ Avg|Chg| no_child one_chil two_chil three_ch Min->Max .27622521 .47115378 -.02081214 .08129664 -.53163827 -+1/2 .01601754 .02867505 -.00315443 .00336003 -.02888066 -+sd/2 .04515528 .08095829 -.00891654 .00935227 -.08139402 MargEfct .06408953 .02867758 -.00315415 .00336719 -.02889062 lesbian Avg|Chg| no_child one_chil two_chil three_ch 0->1 .0735017 .1470034 -.00428884 -.07187283 -.07084173 liberal Avg|Chg| no_child one_chil two_chil three_ch Min->Max .07553434 .09324361 .05321932 .00460574 -.1510687 -+1/2 .01274423 .01557401 .00902894 .00088549 -.02548847 -+sd/2 .01714217 .02094899 .01214351 .00119185 -.03428432 MargEfct .0509872 .01557671 .00903192 .00088496 -.0254936 no_child one_chil two_chil three_ch Pr(y|x) .19485176 .21740308 .29509264 .29265252 age educ lesbian liberal x= 40 12 .024503 2.87682 sd(x)= 17.2202 2.82507 .154657 1.34531 . qui mlogit mchilds age educ lesbian liberal faminc . prchange, x(age 40 educ 12 liberal 5) outcome(2) mlogit: Changes in Predicted Probabilities for mchilds Outcome: 2 (two_children) Min->Max 0->1 -+1/2 -+sd/2 MargEfct age -.2324545 .0008592 -.00375699 -.06236547 .00332766 educ -.00704125 .01010171 -.00223902 -.00640428 -.00754496 lesbian -.00747065 -.00747065 .00548419 .00093047 -.08094006 liberal .05161102 .00949159 .00776719 .01043734 -.00176649 faminc -.09413111 -.00078398 -.00105049 -.0264225 .00188621 no_child one_chil two_chil three_ch Pr(y|x) .20901966 .23686318 .31303602 .24108115 age educ lesbian liberal faminc x= 40 12 .025281 5 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . prchange, x(age 40 educ 12 liberal 5) uncentered d(2) outcome(2) mlogit: Changes in Predicted Probabilities for mchilds Outcome: 2 (two_children) (Note: delta = 2) Min->Max 0->1 +delta +sd MargEfct age -.2324545 .0008592 -.00764541 -.06878462 .00332766 educ -.00704125 .01010171 -.00953369 -.01629095 -.00754496 lesbian -.00747065 -.00747065 -.04444364 .00060984 -.08094006 liberal .05161102 .00949159 .01444076 .00995506 -.00176649 faminc -.09413111 -.00078398 -.00211158 -.02785555 .00188621 no_child one_chil two_chil three_ch Pr(y|x) .20901966 .23686318 .31303602 .24108115 age educ lesbian liberal faminc x= 40 12 .025281 5 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . *check unusual values . qui mlogit dmchilds age educ lesbian liberal faminc . prchange mlogit: Changes in Predicted Probabilities for dmchilds age Avg|Chg| 0 1_5 2_200000 3_5 Min->Max .29928457 -.38933419 -.20923492 .08626884 .51230031 -+1/2 .00530963 -.00700068 -.00361858 .00272629 .00789297 -+sd/2 .08875254 -.11748452 -.06002054 .04538551 .13211958 MargEfct .02123931 -.00700085 -.00361881 .00272645 .00789321 educ Avg|Chg| 0 1_5 2_200000 3_5 Min->Max .28121454 .5199483 .03984125 .00263953 -.56242909 -+1/2 .0188775 .03775498 -.00069709 -.00578019 -.03127772 -+sd/2 .05305353 .10610706 -.00209478 -.01634637 -.08766592 MargEfct .07551504 .03775752 -.00069003 -.0057752 -.03129229 lesbian Avg|Chg| 0 1_5 2_200000 3_5 0->1 .07462971 .14925943 -.00313455 -.08305661 -.06306826 liberal Avg|Chg| 0 1_5 2_200000 3_5 Min->Max .06956485 .09390302 .04522666 -.00136077 -.13776894 -+1/2 .01166802 .01566859 .00766745 -.00012717 -.02320886 -+sd/2 .01567999 .02105704 .01030293 -.00017011 -.03118986 MargEfct .04668151 .01567089 .00766986 -.00012794 -.02321282 faminc Avg|Chg| 0 1_5 2_200000 3_5 Min->Max .12940267 -.17128362 -.08752172 .15182135 .106984 -+1/2 .00152638 -.00206244 -.00099032 .0017707 .00128207 -+sd/2 .03846151 -.05201344 -.02490957 .04461974 .0323033 MargEfct .00610554 -.00206244 -.00099033 .00177071 .00128206 0 1_5 2_200000 3_5 Pr(y|x) .19987729 .20161118 .31913504 .27937651 age educ lesbian liberal faminc x= 44.5442 13.4529 .025281 2.88343 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . qui mlogit ngchilds age educ lesbian liberal faminc . prchange mlogit: Changes in Predicted Probabilities for ngchilds age Avg|Chg| -17 -4 12 78 Min->Max .29928457 -.38933419 -.20923492 .08626884 .51230031 -+1/2 .00530963 -.00700068 -.00361858 .00272629 .00789297 -+sd/2 .08875254 -.11748452 -.06002054 .04538551 .13211958 MargEfct .02123931 -.00700085 -.00361881 .00272645 .00789321 educ Avg|Chg| -17 -4 12 78 Min->Max .28121454 .5199483 .03984125 .00263953 -.56242909 -+1/2 .0188775 .03775498 -.00069709 -.00578019 -.03127772 -+sd/2 .05305353 .10610706 -.00209478 -.01634637 -.08766592 MargEfct .07551504 .03775752 -.00069003 -.0057752 -.03129229 lesbian Avg|Chg| -17 -4 12 78 0->1 .07462971 .14925943 -.00313455 -.08305661 -.06306826 liberal Avg|Chg| -17 -4 12 78 Min->Max .06956485 .09390302 .04522666 -.00136077 -.13776894 -+1/2 .01166802 .01566859 .00766745 -.00012717 -.02320886 -+sd/2 .01567999 .02105704 .01030293 -.00017011 -.03118986 MargEfct .04668151 .01567089 .00766986 -.00012794 -.02321282 faminc Avg|Chg| -17 -4 12 78 Min->Max .12940267 -.17128362 -.08752172 .15182135 .106984 -+1/2 .00152638 -.00206244 -.00099032 .0017707 .00128207 -+sd/2 .03846151 -.05201344 -.02490957 .04461974 .0323033 MargEfct .00610554 -.00206244 -.00099033 .00177071 .00128206 -17 -4 12 78 Pr(y|x) .19987729 .20161118 .31913504 .27937651 age educ lesbian liberal faminc x= 44.5442 13.4529 .025281 2.88343 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . qui mlogit irchilds age educ lesbian liberal faminc . prchange mlogit: Changes in Predicted Probabilities for irchilds age Avg|Chg| 0 3 9 37 Min->Max .29928457 -.38933419 -.20923492 .08626884 .51230031 -+1/2 .00530963 -.00700068 -.00361858 .00272629 .00789297 -+sd/2 .08875254 -.11748452 -.06002054 .04538551 .13211958 MargEfct .02123931 -.00700085 -.00361881 .00272645 .00789321 educ Avg|Chg| 0 3 9 37 Min->Max .28121454 .5199483 .03984125 .00263953 -.56242909 -+1/2 .0188775 .03775498 -.00069709 -.00578019 -.03127772 -+sd/2 .05305353 .10610706 -.00209478 -.01634637 -.08766592 MargEfct .07551504 .03775752 -.00069003 -.0057752 -.03129229 lesbian Avg|Chg| 0 3 9 37 0->1 .07462971 .14925943 -.00313455 -.08305661 -.06306826 liberal Avg|Chg| 0 3 9 37 Min->Max .06956485 .09390302 .04522666 -.00136077 -.13776894 -+1/2 .01166802 .01566859 .00766745 -.00012717 -.02320886 -+sd/2 .01567999 .02105704 .01030293 -.00017011 -.03118986 MargEfct .04668151 .01567089 .00766986 -.00012794 -.02321282 faminc Avg|Chg| 0 3 9 37 Min->Max .12940267 -.17128362 -.08752172 .15182135 .106984 -+1/2 .00152638 -.00206244 -.00099032 .0017707 .00128207 -+sd/2 .03846151 -.05201344 -.02490957 .04461974 .0323033 MargEfct .00610554 -.00206244 -.00099033 .00177071 .00128206 0 3 9 37 Pr(y|x) .19987729 .20161118 .31913504 .27937651 age educ lesbian liberal faminc x= 44.5442 13.4529 .025281 2.88343 35.129 sd(x)= 16.8881 2.81184 .157032 1.34406 25.2427 . . * poisson . qui poisson ctchilds age educ lesbian liberal . prchange, x(age 40 educ 12 liberal 4) poisson: Changes in Predicted Rate for ctchilds min->max 0->1 -+1/2 -+sd/2 MargEfct age 2.1415 0.0137 0.0239 0.4132 0.0239 educ -2.2245 -0.1793 -0.0946 -0.2674 -0.0946 lesbian -0.3949 -0.3949 -0.4474 -0.0690 -0.4461 liberal -0.5131 -0.0961 -0.0813 -0.1093 -0.0813 exp(xb): 1.7006 age educ lesbian liberal x= 40 12 .024503 4 sd(x)= 17.2202 2.82507 .154657 1.34531 . qui poisson ctchilds age educ lesbian liberal . prchange, x(age 40 educ 12 liberal 4) outcome(3) poisson: Changes in Predicted Probabilities for ctchilds Outcome: 3 min->max 0->1 -+1/2 -+sd/2 age 0.1259 0.0017 0.0027 0.0466 educ -0.1481 0.0028 -0.0108 -0.0304 lesbian -0.0489 -0.0489 -0.0503 -0.0079 liberal -0.0544 -0.0086 -0.0093 -0.0125 0 1 2 3 4 5 6 7 Pr(y|x) 0.1826 0.3105 0.2640 0.1497 0.0636 0.0216 0.0061 0.0015 8 9 Pr(y|x) 0.0003 0.0001 age educ lesbian liberal x= 40 12 .024503 4 sd(x)= 17.2202 2.82507 .154657 1.34531 . qui poisson ctchilds age educ lesbian liberal feminist . prchange liberal educ age, x(age 40 educ 12 liberal 4) fromto uncentered poisson: Changes in Predicted Rate for ctchilds from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x liberal 2.0317 1.5439 -0.4878 2.0317 1.9409 -0.0909 1.6919 educ 3.2857 1.0869 -2.1988 3.2857 3.1089 -0.1768 1.6919 age 1.2550 3.2909 2.0359 0.9829 0.9963 0.0134 1.6919 to: dif: from: to: dif: x+1 -+1 x x+sd +sd MargEfct liberal 1.6162 -0.0757 1.6919 1.5904 -0.1015 -0.0774 educ 1.6009 -0.0910 1.6919 1.4522 -0.2397 -0.0936 age 1.7150 0.0231 1.6919 2.1246 0.4327 0.0230 exp(xb): 1.6919 age educ lesbian liberal feminist x= 40 12 .024658 4 .294521 sd(x)= 16.7729 2.76263 .155185 1.35218 .456139 . prchange liberal age educ, x(age 40 educ 12 liberal 4) fromto outcome(2) d(2) poisson: Changes in Predicted Probabilities for ctchilds Outcome: 2 (Note: d = 2) from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x-d/2 liberal 0.2706 0.2545 -0.0161 0.2706 0.2704 -0.0002 0.2669 age 0.2245 0.2015 -0.0230 0.1808 0.1833 0.0025 0.2625 educ 0.2020 0.1992 -0.0027 0.2020 0.2158 0.0138 0.2674 to: dif: from: to: dif: x+d/2 -+d/2 x-1/2sd x+1/2sd -+sd/2 liberal 0.2595 -0.0074 0.2659 0.2609 -0.0050 age 0.2647 0.0022 0.2518 0.2699 0.0181 educ 0.2585 -0.0089 0.2685 0.2562 -0.0123 0 1 2 3 4 5 6 7 Pr(y|x) 0.1842 0.3116 0.2636 0.1487 0.0629 0.0213 0.0060 0.0015 8 9 Pr(y|x) 0.0003 0.0001 age educ lesbian liberal feminist x= 40 12 .024658 4 .294521 sd(x)= 16.7729 2.76263 .155185 1.35218 .456139 . prchange liberal educ age, x(age 40 educ 12 liberal 4) fromto uncentered outc > ome(2) d(2) poisson: Changes in Predicted Probabilities for ctchilds Outcome: 2 (Note: delta = 2) from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x liberal 0.2706 0.2545 -0.0161 0.2706 0.2704 -0.0002 0.2636 educ 0.2020 0.1992 -0.0027 0.2020 0.2158 0.0138 0.2636 age 0.2245 0.2015 -0.0230 0.1808 0.1833 0.0025 0.2636 to: dif: from: to: dif: x+delta +delta x x+sd +sd liberal 0.2545 -0.0091 0.2636 0.2578 -0.0058 educ 0.2522 -0.0114 0.2636 0.2468 -0.0168 age 0.2656 0.0020 0.2636 0.2697 0.0061 0 1 2 3 4 5 6 7 Pr(y|x) 0.1842 0.3116 0.2636 0.1487 0.0629 0.0213 0.0060 0.0015 8 9 Pr(y|x) 0.0003 0.0001 age educ lesbian liberal feminist x= 40 12 .024658 4 .294521 sd(x)= 16.7729 2.76263 .155185 1.35218 .456139 . . * nbreg . qui nbreg ctchilds age educ lesbian liberal . prchange, x(age 40 educ 12 liberal 3) nbreg: Changes in Predicted Rate for ctchilds min->max 0->1 -+1/2 -+sd/2 MargEfct age 2.3033 0.0145 0.0256 0.4419 0.0256 educ -2.3522 -0.1902 -0.0998 -0.2823 -0.0998 lesbian -0.4116 -0.4116 -0.4660 -0.0719 -0.4647 liberal -0.5155 -0.0966 -0.0856 -0.1152 -0.0856 exp(xb): 1.7806 age educ lesbian liberal x= 40 12 .024503 3 sd(x)= 17.2202 2.82507 .154657 1.34531 . qui nbreg ctchilds age educ lesbian liberal feminist . prchange, x(age 40 educ 12 liberal 7) outcome(8) nbreg: Changes in Predicted Probabilities for ctchilds Outcome: 8 min->max 0->1 -+1/2 -+sd/2 age 0.0093 0.0000 0.0000 0.0004 educ -0.0089 -0.0020 -0.0001 -0.0002 lesbian -0.0001 -0.0001 -0.0002 0.0000 liberal -0.0012 -0.0004 -0.0001 -0.0001 feminist -0.0001 -0.0001 -0.0001 -0.0001 0 1 2 3 4 5 6 7 Pr(y|x) 0.2403 0.3328 0.2399 0.1198 0.0466 0.0150 0.0042 0.0010 8 9 Pr(y|x) 0.0002 0.0000 age educ lesbian liberal feminist x= 40 12 .024658 7 .294521 sd(x)= 16.7729 2.76263 .155185 1.35218 .456139 . qui nbreg ctchilds age educ lesbian liberal feminist . prchange, x(age 40 educ 12 liberal 7) fromto nbreg: Changes in Predicted Rate for ctchilds from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x-1/2 age 1.0811 2.9027 1.8215 0.8416 0.8534 0.0118 1.4580 educ 2.8731 0.9384 -1.9347 2.8731 2.7168 -0.1563 1.5098 lesbian 1.4722 1.3165 -0.1558 1.4722 1.3165 -0.1558 1.5526 liberal 2.0325 1.5380 -0.4945 2.0325 1.9402 -0.0923 1.5027 feminist 1.5100 1.3726 -0.1375 1.5100 1.3726 -0.1375 1.5399 to: dif: from: to: dif: x+1/2 -+1/2 x-1/2sd x+1/2sd -+sd/2 MargEfct age 1.4784 0.0204 1.3065 1.6498 0.3433 0.0204 educ 1.4277 -0.0822 1.5861 1.3590 -0.2272 -0.0821 lesbian 1.3883 -0.1643 1.4810 1.4555 -0.0255 -0.1642 liberal 1.4345 -0.0682 1.5150 1.4228 -0.0923 -0.0682 feminist 1.3997 -0.1402 1.5005 1.4366 -0.0639 -0.1401 exp(xb): 1.4682 age educ lesbian liberal feminist x= 40 12 .024658 7 .294521 sd(x)= 16.7729 2.76263 .155185 1.35218 .456139 . prchange educ age liberal, x(age 40 educ 12 liberal 3) d(2) nbreg: Changes in Predicted Rate for ctchilds (Note: d = 2) min->max 0->1 -+d/2 -+sd/2 MargEfct educ -2.3298 -0.1883 -0.1979 -0.2735 -0.0989 age 2.1936 0.0142 0.0492 0.4135 0.0246 liberal -0.4945 -0.0923 -0.1644 -0.1111 -0.0821 exp(xb): 1.7680 age educ lesbian liberal feminist x= 40 12 .024658 3 .294521 sd(x)= 16.7729 2.76263 .155185 1.35218 .456139 . prchange liberal age educ, x(age 40 educ 12 liberal 3) uncentered d(2) nbreg: Changes in Predicted Rate for ctchilds (Note: delta = 2) min->max 0->1 +delta +sd MargEfct liberal -0.4945 -0.0923 -0.1569 -0.1077 -0.0821 age 2.1936 0.0142 0.0499 0.4646 0.0246 educ -2.3298 -0.1883 -0.1872 -0.2532 -0.0989 exp(xb): 1.7680 age educ lesbian liberal feminist x= 40 12 .024658 3 .294521 sd(x)= 16.7729 2.76263 .155185 1.35218 .456139 . . * zip . qui zip zchilds age educ lesbian liberal, inf(age faminc liberal) . prchange, x(age 40 educ median faminc 45) zip: Changes in Predicted Rate for zchilds min->max 0->1 -+1/2 -+sd/2 age 1.1815 0.0038 0.0158 0.2655 educ -1.4232 -0.1585 -0.0444 -0.1251 lesbian 0.1017 0.1017 0.0920 0.0144 liberal -0.3233 -0.0682 -0.0540 -0.0726 exp(xb): 0.4355 base x values for count equation: age educ lesbian liberal x= 40 13 .025281 2.88343 sd(x)= 16.8881 2.81184 .157032 1.34406 base x values for binary equation: age faminc liberal x= 40 45 2.88343 sd(x)= 16.8881 25.2427 1.34406 . prchange, x(age max educ 12 faminc mean) outcome(0) zip: Changes in Predicted Probabilities for zchilds Outcome: 0 min->max 0->1 -+1/2 -+sd/2 age -0.6240 -0.0005 -0.0087 -0.1471 educ 0.4310 0.0035 0.0288 0.0807 lesbian -0.0564 -0.0564 -0.0596 -0.0094 liberal 0.1888 0.0282 0.0316 0.0424 0 1 2 3 4 5 6 7 Pr(y|x) 0.3523 0.2831 0.2090 0.1029 0.0380 0.0112 0.0028 0.0006 8 9 Pr(y|x) 0.0001 0.0000 base x values for count equation: age educ lesbian liberal x= 89 12 .025281 2.88343 sd(x)= 16.8881 2.81184 .157032 1.34406 base x values for binary equation: age faminc liberal x= 89 13.4529 2.88343 sd(x)= 16.8881 25.2427 1.34406 . prchange educ faminc, x(age max educ 12 faminc mean) uncentered outcome(0) zip: Changes in Predicted Probabilities for zchilds Outcome: 0 min->max 0->1 +1 +sd educ 0.4310 0.0035 0.0295 0.0852 faminc -0.0437 -0.0005 -0.0005 -0.0129 0 1 2 3 4 5 6 7 Pr(y|x) 0.3523 0.2831 0.2090 0.1029 0.0380 0.0112 0.0028 0.0006 8 9 Pr(y|x) 0.0001 0.0000 base x values for count equation: age educ lesbian liberal x= 89 12 .025281 2.88343 sd(x)= 16.8881 2.81184 .157032 1.34406 base x values for binary equation: age faminc liberal x= 89 13.4529 2.88343 sd(x)= 16.8881 25.2427 1.34406 . prchange educ faminc, x(age max educ 12 faminc mean) uncentered outcome(0) d( > 4) zip: Changes in Predicted Probabilities for zchilds Outcome: 0 (Note: delta = 4) min->max 0->1 +delta +sd educ 0.4310 0.0035 0.1225 0.0852 faminc -0.0437 -0.0005 -0.0021 -0.0129 0 1 2 3 4 5 6 7 Pr(y|x) 0.3523 0.2831 0.2090 0.1029 0.0380 0.0112 0.0028 0.0006 8 9 Pr(y|x) 0.0001 0.0000 base x values for count equation: age educ lesbian liberal x= 89 12 .025281 2.88343 sd(x)= 16.8881 2.81184 .157032 1.34406 base x values for binary equation: age faminc liberal x= 89 13.4529 2.88343 sd(x)= 16.8881 25.2427 1.34406 . . * zinb . qui zinb zchilds age educ lesbian liberal, inf(age educ faminc) . prchange, x(age 40 educ 12 faminc 45) zinb: Changes in Predicted Rate for zchilds min->max 0->1 -+1/2 -+sd/2 age 1.1371 0.0005 0.0266 0.4457 educ -2.4534 -0.2332 -0.1051 -0.2951 lesbian 0.0211 0.0211 0.0208 0.0033 liberal -0.4069 -0.0860 -0.0676 -0.0909 exp(xb): 0.6680 base x values for count equation: age educ lesbian liberal x= 40 12 .025281 2.88343 sd(x)= 16.8881 2.81184 .157032 1.34406 base x values for binary equation: age educ faminc x= 40 12 45 sd(x)= 16.8881 2.81184 25.2427 . qui zinb zchilds age educ lesbian liberal, inf(age educ faminc liberal) . prchange, x(age 40 educ 12 faminc 60) outcome(4) zinb: Changes in Predicted Probabilities for zchilds Outcome: 4 min->max 0->1 -+1/2 -+sd/2 age 0.0408 0.0000 0.0014 0.0202 educ -0.0730 -0.0034 -0.0048 -0.0134 lesbian 0.0011 0.0011 0.0011 0.0002 liberal -0.0196 -0.0041 -0.0034 -0.0046 0 1 2 3 4 5 6 7 Pr(y|x) 0.6649 0.1996 0.0807 0.0326 0.0132 0.0054 0.0022 0.0009 8 9 Pr(y|x) 0.0004 0.0001 base x values for count equation: age educ lesbian liberal x= 40 12 .025281 2.88343 sd(x)= 16.8881 2.81184 .157032 1.34406 base x values for binary equation: age educ faminc liberal x= 40 12 60 2.88343 sd(x)= 16.8881 2.81184 25.2427 1.34406 . prchange educ faminc, x(age 40 educ 12 faminc 45) fromto d(4) zinb: Changes in Predicted Rate for zchilds (Note: d = 4) from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x-d/2 educ 2.4249 0.0778 -2.3471 2.4249 2.2115 -0.2134 0.8895 faminc 0.6649 0.6834 0.0185 0.6648 0.6650 0.0002 0.6739 to: dif: from: to: dif: x+d/2 -+d/2 x-1/2sd x+1/2sd -+sd/2 educ 0.4723 -0.4172 0.8243 0.5306 -0.2937 faminc 0.6747 0.0008 0.6717 0.6769 0.0052 exp(xb): 0.6743 base x values for count equation: age educ lesbian liberal x= 40 12 .025281 2.88343 sd(x)= 16.8881 2.81184 .157032 1.34406 base x values for binary equation: age educ faminc liberal x= 40 12 45 2.88343 sd(x)= 16.8881 2.81184 25.2427 1.34406 . prchange educ faminc, x(age 40 educ 12 faminc 45) fromto d(4) uncentered zinb: Changes in Predicted Rate for zchilds (Note: delta = 4) from: to: dif: from: to: dif: from: x=min x=max min->max x=0 x=1 0->1 x educ 2.4249 0.0778 -2.3471 2.4249 2.2115 -0.2134 0.6743 faminc 0.6649 0.6834 0.0185 0.6648 0.6650 0.0002 0.6743 to: dif: from: to: dif: x+delta +delta x x+sd +sd educ 0.2948 -0.3795 0.6743 0.3962 -0.2781 faminc 0.6752 0.0008 0.6743 0.6795 0.0051 exp(xb): 0.6743 base x values for count equation: age educ lesbian liberal x= 40 12 .025281 2.88343 sd(x)= 16.8881 2.81184 .157032 1.34406 base x values for binary equation: age educ faminc liberal x= 40 12 45 2.88343 sd(x)= 16.8881 2.81184 25.2427 1.34406 . . * Testing marginal effects with output from LIMDEP 7 . * RM4CLDV Ch 3: Binary logit - verified 12/6/00 . use binlfp2,clear (PSID 1976 / T Mroz) . logit lfp k5 k618 age wc hc lwg inc, nolog Logit estimates Number of obs = 753 LR chi2(7) = 124.48 Prob > chi2 = 0.0000 Log likelihood = -452.63296 Pseudo R2 = 0.1209 ------------------------------------------------------------------------------ lfp | Coef. Std. Err. z P>|z| [95% Conf. Interval] ---------+-------------------------------------------------------------------- k5 | -1.462913 .1970006 -7.426 0.000 -1.849027 -1.076799 k618 | -.0645707 .0680008 -0.950 0.342 -.1978499 .0687085 age | -.0628706 .0127831 -4.918 0.000 -.0879249 -.0378162 wc | .8072738 .2299799 3.510 0.000 .3565215 1.258026 hc | .1117336 .2060397 0.542 0.588 -.2920969 .515564 lwg | .6046931 .1508176 4.009 0.000 .3090961 .9002901 inc | -.0344464 .0082084 -4.196 0.000 -.0505346 -.0183583 _cons | 3.18214 .6443751 4.938 0.000 1.919188 4.445092 ------------------------------------------------------------------------------ . prchange logit: Changes in Predicted Probabilities for lfp min->max 0->1 -+1/2 -+sd/2 MargEfct k5 -0.6361 -0.3499 -0.3428 -0.1849 -0.3569 k618 -0.1278 -0.0156 -0.0158 -0.0208 -0.0158 age -0.4372 -0.0030 -0.0153 -0.1232 -0.0153 wc 0.1881 0.1881 0.1945 0.0884 0.1969 hc 0.0272 0.0272 0.0273 0.0133 0.0273 lwg 0.6624 0.1499 0.1465 0.0865 0.1475 inc -0.6415 -0.0068 -0.0084 -0.0975 -0.0084 NotInLF inLF Pr(y|x) 0.4222 0.5778 k5 k618 age wc hc lwg inc x= .237716 1.35325 42.5378 .281541 .391766 1.09711 20.129 sd(x)= .523959 1.31987 8.07257 .450049 .488469 .587556 11.6348 . /* > +---------------------------------------------+ > | Multinomial Logit Model | > +---------------------------------------------+ > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Characteristics in numerator of Prob[Y = 1] > Constant 3.182140457 .64437510 4.938 .0000 > K5 -1.462913039 .19700061 -7.426 .0000 .23771580 > K618 -.6457068349E-01 .68000829E-01 -.950 .3423 1.3532537 > AGE -.6287055111E-01 .12783091E-01 -4.918 .0000 42.537849 > WC .8072737782 .22997989 3.510 .0004 .28154050 > HC .1117335717 .20603972 .542 .5876 .39176627 > LWG .6046931218 .15081757 4.009 .0001 1.0971148 > INC -.3444643073E-01 .82083765E-02 -4.196 .0000 20.128965 > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Marginal effects on Prob[Y = 1] > Constant .7762769957 .15584838 4.981 .0000 > K5 -.3568747999 .48210796E-01 -7.402 .0000 .23771580 > K618 -.1575189306E-01 .16589340E-01 -.950 .3424 1.3532537 > AGE -.1533714906E-01 .31108547E-02 -4.930 .0000 42.537849 > WC .1969328732 .55947366E-01 3.520 .0004 .28154050 > HC .2725718822E-01 .50271630E-01 .542 .5877 .39176627 > LWG .1475137148 .36744550E-01 4.015 .0001 1.0971148 > INC -.8403140000E-02 .20046214E-02 -4.192 .0000 20.128965 > */ . * RM4CLDV Ch 3: Binary probit - verified 12/6/00 . use binlfp2,clear (PSID 1976 / T Mroz) . probit lfp k5 k618 age wc hc lwg inc, nolog Probit estimates Number of obs = 753 LR chi2(7) = 124.36 Prob > chi2 = 0.0000 Log likelihood = -452.69496 Pseudo R2 = 0.1208 ------------------------------------------------------------------------------ lfp | Coef. Std. Err. z P>|z| [95% Conf. Interval] ---------+-------------------------------------------------------------------- k5 | -.8747112 .1135583 -7.703 0.000 -1.097281 -.6521411 k618 | -.0385945 .0404893 -0.953 0.340 -.117952 .0407631 age | -.0378235 .0076093 -4.971 0.000 -.0527375 -.0229095 wc | .4883144 .1354873 3.604 0.000 .2227642 .7538645 hc | .0571704 .1240052 0.461 0.645 -.1858754 .3002161 lwg | .3656287 .0877792 4.165 0.000 .1935847 .5376727 inc | -.020525 .0047769 -4.297 0.000 -.0298875 -.0111626 _cons | 1.918422 .3806536 5.040 0.000 1.172355 2.66449 ------------------------------------------------------------------------------ . prchange probit: Changes in Predicted Probabilities for lfp min->max 0->1 -+1/2 -+sd/2 MargEfct k5 -0.6441 -0.3380 -0.3320 -0.1778 -0.3422 k618 -0.1221 -0.0150 -0.0151 -0.0199 -0.0151 age -0.4274 -0.0031 -0.0148 -0.1190 -0.0148 wc 0.1844 0.1844 0.1892 0.0858 0.1911 hc 0.0223 0.0223 0.0224 0.0109 0.0224 lwg 0.6649 0.1450 0.1423 0.0839 0.1431 inc -0.6425 -0.0068 -0.0080 -0.0932 -0.0080 NotInLF inLF Pr(y|x) 0.4218 0.5782 k5 k618 age wc hc lwg inc x= .237716 1.35325 42.5378 .281541 .391766 1.09711 20.129 sd(x)= .523959 1.31987 8.07257 .450049 .488469 .587556 11.6348 . /* > +---------------------------------------------+ > | Binomial Probit Model | > +---------------------------------------------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Index function for probability > Constant 1.918422308 .38065387 5.040 .0000 > K5 -.8747111768 .11355839 -7.703 .0000 .23771580 > K618 -.3859449198E-01 .40489311E-01 -.953 .3405 1.3532537 > AGE -.3782350247E-01 .76093431E-02 -4.971 .0000 42.537849 > WC .4883143923 .13548734 3.604 .0003 .28154050 > HC .5717035056E-01 .12400528 .461 .6448 .39176627 > LWG .3656287134 .87779191E-01 4.165 .0000 1.0971148 > INC -.2052502534E-01 .47768634E-02 -4.297 .0000 20.128965 > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Index function for probability > Constant .7506026112 .14798760 5.072 .0000 > K5 -.3422398139 .44525077E-01 -7.686 .0000 .23771580 > K618 -.1510049500E-01 .15840334E-01 -.953 .3404 1.3532537 > AGE -.1479883736E-01 .29733362E-02 -4.977 .0000 42.537849 > WC .1910580671 .52960432E-01 3.608 .0003 .28154050 > HC .2236849220E-01 .48523420E-01 .461 .6448 .39176627 > LWG .1430560237 .34259780E-01 4.176 .0000 1.0971148 > INC -.8030628895E-02 .18697205E-02 -4.295 .0000 20.128965 > */ . * RM4CLDV Ch 5: Ordinal Logit - verified 12/6/00 . use ordwarm2,clear (77 & 89 General Social Survey) . ologit warm yr89 male white age ed prst, nolog Ordered logit estimates Number of obs = 2293 LR chi2(6) = 301.72 Prob > chi2 = 0.0000 Log likelihood = -2844.9123 Pseudo R2 = 0.0504 ------------------------------------------------------------------------------ warm | Coef. Std. Err. z P>|z| [95% Conf. Interval] ---------+-------------------------------------------------------------------- yr89 | .5239025 .0798988 6.557 0.000 .3673037 .6805013 male | -.7332997 .0784827 -9.343 0.000 -.8871229 -.5794766 white | -.3911595 .1183808 -3.304 0.001 -.6231815 -.1591374 age | -.0216655 .0024683 -8.778 0.000 -.0265032 -.0168278 ed | .0671728 .015975 4.205 0.000 .0358624 .0984831 prst | .0060727 .0032929 1.844 0.065 -.0003813 .0125267 ---------+-------------------------------------------------------------------- _cut1 | -2.465362 .2389126 (Ancillary parameters) _cut2 | -.630904 .2333155 _cut3 | 1.261854 .2340179 ------------------------------------------------------------------------------ . prchange ologit: Changes in Predicted Probabilities for warm yr89 Avg|Chg| SD D A SA 0->1 .06372773 -.04993673 -.07751873 .05391335 .0735421 male Avg|Chg| SD D A SA 0->1 .08961766 .07461427 .10462105 -.08137083 -.09786447 white Avg|Chg| SD D A SA 0->1 .04694162 .03452424 .05935898 -.03557652 -.05830672 age Avg|Chg| SD D A SA Min->Max .18319855 .18012119 .1862759 -.17905769 -.18733941 -+1/2 .00266841 .00214228 .00319454 -.00240716 -.00292964 -+sd/2 .0446563 .03602537 .05328724 -.0401054 -.0492072 MargEfct .01067375 .00214226 .00319461 -.00240723 -.00292964 ed Avg|Chg| SD D A SA Min->Max .16156666 -.15235534 -.17077798 .15631917 .16681415 -+1/2 .00827261 -.00664248 -.00990272 .0074617 .00908352 -+sd/2 .02612717 -.02101015 -.03124419 .02353376 .02872056 MargEfct .03309338 -.00664198 -.00990472 .0074635 .00908319 prst Avg|Chg| SD D A SA Min->Max .05186236 -.04070839 -.06301633 .04440692 .05931778 -+1/2 .00074795 -.00060046 -.00089544 .00067475 .00082116 -+sd/2 .01083777 -.00870322 -.01297233 .00977433 .0119012 MargEfct .00299179 -.00060046 -.00089543 .00067473 .00082116 SD D A SA Pr(y|x) .11125716 .32816544 .39936733 .16121005 yr89 male white age ed prst x= .398604 .464893 .876581 44.9355 12.2181 39.5853 sd(x)= .489718 .498875 .328989 16.779 3.16083 14.4923 . /* > +---------------------------------------------+ > | Ordered Probit Model | > | Logistic Probability Model | > +---------------------------------------------+ > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Index function for probability > Constant 2.465362053 .24019735 10.264 .0000 > YR89 .5239024812 .82726851E-01 6.333 .0000 .39860445 > MALE -.7332997387 .79453652E-01 -9.229 .0000 .46489315 > WHITE -.3911594756 .11968218 -3.268 .0011 .87658090 > AGE -.2166550942E-01 .25416098E-02 -8.524 .0000 44.935456 > ED .6717275019E-01 .15977684E-01 4.204 .0000 12.218055 > PRST .6072712231E-02 .33018864E-02 1.839 .0659 39.585259 > Threshold parameters for index > Mu( 1) 1.834457988 .63291766E-01 28.984 .0000 > Mu( 2) 3.727215868 .82017394E-01 45.444 .0000 > > +------------------------------------------------------+ > | Marginal Effects for OrdLogit | > +----------+----------+----------+----------+----------+ > | Variable | WARM0=0 | WARM0=1 | WARM0=2 | WARM0=3 | > +----------+----------+----------+----------+----------+ > | ONE | -.2438 | -.3635 | .2739 | .3334 | > | YR89 | -.0518 | -.0773 | .0582 | .0708 | > | MALE | .0725 | .1081 | -.0815 | -.0992 | > | WHITE | .0387 | .0577 | -.0435 | -.0529 | > | AGE | .0021 | .0032 | -.0024 | -.0029 | > | ED | -.0066 | -.0099 | .0075 | .0091 | > | PRST | -.0006 | -.0009 | .0007 | .0008 | > +----------+----------+----------+----------+----------+ > */ . * RM4CLDV Ch 5: Ordinal Probit - verified 12/6/00 . use ordwarm2,clear (77 & 89 General Social Survey) . oprobit warm yr89 male white age ed prst, nolog Ordered probit estimates Number of obs = 2293 LR chi2(6) = 294.32 Prob > chi2 = 0.0000 Log likelihood = -2848.611 Pseudo R2 = 0.0491 ------------------------------------------------------------------------------ warm | Coef. Std. Err. z P>|z| [95% Conf. Interval] ---------+-------------------------------------------------------------------- yr89 | .3188147 .0468519 6.805 0.000 .2269867 .4106427 male | -.4170287 .0455459 -9.156 0.000 -.5062971 -.3277603 white | -.2265002 .0694773 -3.260 0.001 -.3626733 -.0903272 age | -.0122213 .0014427 -8.471 0.000 -.0150489 -.0093937 ed | .0387234 .0093241 4.153 0.000 .0204485 .0569982 prst | .003283 .001925 1.705 0.088 -.0004899 .0070559 ---------+-------------------------------------------------------------------- _cut1 | -1.428578 .1387742 (Ancillary parameters) _cut2 | -.3605589 .1369219 _cut3 | .7681637 .1370564 ------------------------------------------------------------------------------ . prchange oprobit: Changes in Predicted Probabilities for warm yr89 Avg|Chg| SD D A SA 0->1 .06230063 -.05840641 -.06619486 .04342893 .08117232 male Avg|Chg| SD D A SA 0->1 .08185343 .08099359 .08271328 -.06222743 -.10147941 white Avg|Chg| SD D A SA 0->1 .04386704 .03874894 .04898515 -.0271824 -.06055167 age Avg|Chg| SD D A SA Min->Max .16764789 .18767587 .1476199 -.14207253 -.19322326 -+1/2 .00241081 .00232452 .00249711 -.00180408 -.00301754 -+sd/2 .04038233 .03903589 .04172876 -.03013682 -.05062783 MargEfct .00964335 .00232452 .00249716 -.00180413 -.00301755 ed Avg|Chg| SD D A SA Min->Max .1504221 -.16398501 -.13685921 .12518799 .1756562 -+1/2 .0076383 -.00736549 -.00791112 .00571549 .00956112 -+sd/2 .02413012 -.02328731 -.02497292 .01803988 .03022036 MargEfct .03055511 -.00736526 -.00791229 .00571641 .00956114 prst Avg|Chg| SD D A SA Min->Max .04505957 -.04247186 -.04764727 .03203639 .05808274 -+1/2 .00064762 -.00062443 -.00067082 .00048462 .00081059 -+sd/2 .00938461 -.00904983 -.00971937 .00702184 .01174738 MargEfct .00259048 -.00062443 -.00067081 .00048464 .0008106 SD D A SA Pr(y|x) .11177309 .32895118 .39563131 .1636444 yr89 male white age ed prst x= .398604 .464893 .876581 44.9355 12.2181 39.5853 sd(x)= .489718 .498875 .328989 16.779 3.16083 14.4923 . /* > +---------------------------------------------+ > | Ordered Probit Model | > +---------------------------------------------+ > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Index function for probability > Constant 1.428578307 .13929995 10.255 .0000 > YR89 .3188147042 .48583695E-01 6.562 .0000 .39860445 > MALE -.4170287097 .46041475E-01 -9.058 .0000 .46489315 > WHITE -.2265002456 .70073897E-01 -3.232 .0012 .87658090 > AGE -.1222129917E-01 .14782079E-02 -8.268 .0000 44.935456 > ED .3872338654E-01 .92824800E-02 4.172 .0000 12.218055 > PRST .3282991562E-02 .19192003E-02 1.711 .0872 39.585259 > Threshold parameters for index > Mu( 1) 1.068019409 .34793659E-01 30.696 .0000 > Mu( 2) 2.196742035 .43462929E-01 50.543 .0000 > +------------------------------------------------------+ > | Marginal Effects for OrdProbt | > +----------+----------+----------+----------+----------+ > | Variable | WARM0=0 | WARM0=1 | WARM0=2 | WARM0=3 | > +----------+----------+----------+----------+----------+ > | ONE | -.2717 | -.2919 | .2109 | .3527 | > | YR89 | -.0606 | -.0651 | .0471 | .0787 | > | MALE | .0793 | .0852 | -.0616 | -.1030 | > | WHITE | .0431 | .0463 | -.0334 | -.0559 | > | AGE | .0023 | .0025 | -.0018 | -.0030 | > | ED | -.0074 | -.0079 | .0057 | .0096 | > | PRST | -.0006 | -.0007 | .0005 | .0008 | > +----------+----------+----------+----------+----------+ > */ . * RM4CLDV Ch 6: Nominal Outcomes . * LIMDEP would not run?? . . * RM4CLDV Ch 8: PRM - verified 12/6/00 . use couart2, clear (Academic Biochemists / S Long) . qui poisson art fem mar kid5 phd ment, nolog . prchange poisson: Changes in Predicted Rate for art min->max 0->1 -+1/2 -+sd/2 MargEfct fem -0.3591 -0.3591 -0.3624 -0.1804 -0.3616 mar 0.2440 0.2440 0.2502 0.1183 0.2500 kid5 -0.7512 -0.2978 -0.2981 -0.2279 -0.2977 phd 0.0794 0.0200 0.0206 0.0203 0.0206 ment 7.9124 0.0333 0.0411 0.3910 0.0411 exp(xb): 1.6101 fem mar kid5 phd ment x= .460109 .662295 .495082 3.10311 8.76721 sd(x)= .498679 .473186 .76488 .984249 9.48392 . /* > LIMDEP marginals differ slightly since they seem to be based on different p > redicted > value. I think SPOST is correct. --jsl > +---------------------------------------------+ > | Poisson Regression | > +---------------------------------------------+ > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Constant .3046168422 .10298215 2.958 .0031 > FEM -.2245942259 .54613757E-01 -4.112 .0000 .46010929 > MAR .1552433860 .61374690E-01 2.529 .0114 .66229508 > KID5 -.1848827020 .40127172E-01 -4.607 .0000 .49508197 > PHD .1282257535E-01 .26397193E-01 .486 .6271 3.1031093 > MENT .2554274773E-01 .20060778E-02 12.733 .0000 8.7672125 > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Constant .5156846879 .18742784 2.751 .0059 > FEM -.3802147066 .99397208E-01 -3.825 .0001 .46010929 > MAR .2628109348 .11170213 2.353 .0186 .66229508 > KID5 -.3129872196 .73031579E-01 -4.286 .0000 .49508197 > PHD .2170728880E-01 .48042973E-01 .452 .6514 3.1031093 > MENT .4324122001E-01 .36510678E-02 11.843 .0000 8.7672125 > */ . * RM4CLDV Ch 8: NBRM - verified 12/6/00 . use couart2, clear (Academic Biochemists / S Long) . qui nbreg art fem mar kid5 phd ment, nolog . prchange nbreg: Changes in Predicted Rate for art min->max 0->1 -+1/2 -+sd/2 MargEfct fem -0.3444 -0.3444 -0.3474 -0.1730 -0.3467 mar 0.2355 0.2355 0.2413 0.1141 0.2411 kid5 -0.7184 -0.2827 -0.2830 -0.2163 -0.2826 phd 0.0940 0.0235 0.0245 0.0241 0.0245 ment 10.4123 0.0366 0.0466 0.4433 0.0466 exp(xb): 1.6020 fem mar kid5 phd ment x= .460109 .662295 .495082 3.10311 8.76721 sd(x)= .498679 .473186 .76488 .984249 9.48392 . /* > LIMDEP marginals differ slightly since they seem to be based on different p > redicted > value. I think SPOST is correct. --jsl > +---------------------------------------------+ > | Negative Binomial Regression | > +---------------------------------------------+ > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Constant .2561440419 .14088343 1.818 .0690 > FEM -.2164184294 .76725070E-01 -2.821 .0048 .46010929 > MAR .1504894586 .84210399E-01 1.787 .0739 .66229508 > KID5 -.1764152459 .53828662E-01 -3.277 .0010 .49508197 > PHD .1527115579E-01 .36258477E-01 .421 .6736 3.1031093 > MENT .2908234254E-01 .31962782E-02 9.099 .0000 8.7672125 > Overdispersion parameter for negative binomial model > Alpha .4416204769 .52241406E-01 8.453 .0000 > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Constant .4367074542 .26508492 1.647 .0995 > FEM -.3689780977 .14436516 -2.556 .0106 .46010929 > MAR .2565738709 .15844949 1.619 .1054 .66229508 > KID5 -.3007755026 .10128350 -2.970 .0030 .49508197 > PHD .2603623929E-01 .68223606E-01 .382 .7027 3.1031093 > MENT .4958333474E-01 .60140867E-02 8.245 .0000 8.7672125 > */ . * RM4CLDV Ch 8: ZIP . use couart2, clear (Academic Biochemists / S Long) . qui zip art fem mar kid5 phd ment, inf(fem mar kid5 phd ment) vuong nolog . prchange zip: Changes in Predicted Rate for art min->max 0->1 -+1/2 -+sd/2 fem -0.3796 -0.3796 -0.3829 -0.1907 mar 0.2581 0.2581 0.2604 0.1232 kid5 -0.7736 -0.2954 -0.2956 -0.2261 phd -0.0419 -0.0110 -0.0108 -0.0106 ment 5.6902 0.0710 0.0625 0.5940 exp(xb): 1.7032 base x values for count equation: fem mar kid5 phd ment x= .460109 .662295 .495082 3.10311 8.76721 sd(x)= .498679 .473186 .76488 .984249 9.48392 base x values for binary equation: fem mar kid5 phd ment x= .460109 .662295 .495082 3.10311 8.76721 sd(x)= .498679 .473186 .76488 .984249 9.48392 . /* > +----------------------------------------------------------------------+ > | Zero Altered Poisson Regression Model | > | Logistic distribution used for splitting model. | > | Vuong statistic for testing ZIP vs. unaltered model is 5.9814 | > | Distributed as standard normal. A value greater than | > | +1.96 favors the zero altered Z.I.Poisson model. | > | A value less than -1.96 rejects the ZIP model. | > +----------------------------------------------------------------------+ > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Poisson/Negbin regression model > Constant .6408380343 .93871217E-01 6.827 .0000 > FEM -.2091445779 .53192738E-01 -3.932 .0001 .46010929 > MAR .1037509411 .57307146E-01 1.810 .0702 .66229508 > KID5 -.1433196644 .36169739E-01 -3.962 .0001 .49508197 > PHD -.6166062012E-02 .22554341E-01 -.273 .7846 3.1031093 > MENT .1809772537E-01 .14247000E-02 12.703 .0000 8.7672125 > Zero inflation model > Constant -.5770603212 .50202649 -1.149 .2504 > FEM .1097471967 .27754038 .395 .6925 .46010929 > MAR -.3540134918 .30942379 -1.144 .2526 .66229508 > KID5 .2171006204 .18853735 1.151 .2495 .49508197 > PHD .1272250778E-02 .13287880 .010 .9924 3.1031093 > MENT -.1341135368 .31257044E-01 -4.291 .0000 8.7672125 > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Constant 1.328163759 .21308395 6.233 .0000 > FEM -.4132860761 .11147948 -3.707 .0002 .46010929 > MAR .2819299084 .12056336 2.338 .0194 .66229508 > KID5 -.3196400023 .76448810E-01 -4.181 .0000 .49508197 > PHD -.1168056422E-01 .49914435E-01 -.234 .8150 3.1031093 > MENT .6775586547E-01 .10035190E-01 6.752 .0000 8.7672125 > */ . * RM4CLDV Ch 8: ZINB . use couart2, clear (Academic Biochemists / S Long) . qui zinb art fem mar kid5 phd ment, inf(fem mar kid5 phd ment) vuong nolog . prchange zinb: Changes in Predicted Rate for art min->max 0->1 -+1/2 -+sd/2 fem -0.3307 -0.3307 -0.3333 -0.1660 mar 0.1642 0.1642 0.1667 0.0788 kid5 -0.6716 -0.2584 -0.2586 -0.1977 phd -0.0045 -0.0012 -0.0012 -0.0012 ment 8.3177 0.2496 0.0425 0.4221 exp(xb): 1.7010 base x values for count equation: fem mar kid5 phd ment x= .460109 .662295 .495082 3.10311 8.76721 sd(x)= .498679 .473186 .76488 .984249 9.48392 base x values for binary equation: fem mar kid5 phd ment x= .460109 .662295 .495082 3.10311 8.76721 sd(x)= .498679 .473186 .76488 .984249 9.48392 . . /* > +----------------------------------------------------------------------+ > | Zero Altered Neg.Binomial Regression Model | > | Logistic distribution used for splitting model. | > | Vuong statistic for testing ZIP vs. unaltered model is 2.3185 | > | Distributed as standard normal. A value greater than | > | +1.96 favors the zero altered Z.I.Neg_Bin model. | > | A value less than -1.96 rejects the ZIP model. | > +----------------------------------------------------------------------+ > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Poisson/Negbin regression model > Constant .4167466020 .14757670 2.824 .0047 > FEM -.1955068333 .78739846E-01 -2.483 .0130 .46010929 > MAR .9758258391E-01 .85764108E-01 1.138 .2552 .66229508 > KID5 -.1517324502 .54075119E-01 -2.806 .0050 .49508197 > PHD -.7001514503E-03 .36482078E-01 -.019 .9847 3.1031093 > MENT .2478620306E-01 .31105703E-02 7.968 .0000 8.7672125 > Overdispersion parameter > Alpha .3766810687 .50499088E-01 7.459 .0000 > Zero inflation model > Constant -.1916860235 1.1904273 -.161 .8721 > FEM .6359326068 .75194896 .846 .3977 .46010929 > MAR -1.499469145 .82562789 -1.816 .0693 .66229508 > KID5 .6284274950 .45924528 1.368 .1712 .49508197 > PHD -.3771530810E-01 .32187945 -.117 .9067 3.1031093 > MENT -.8822934899 .37812504 -2.333 .0196 1.0000000 > +---------+--------------+----------------+--------+---------+----------+ > |Variable | Coefficient | Standard Error |b/St.Er.|P[|Z|>z] | Mean of X| > +---------+--------------+----------------+--------+---------+----------+ > Constant .6769194657 .31020676 2.182 .0291 > FEM -.2515104394 .17217705 -1.461 .1441 .46010929 > MAR .1321340645 .21811993 .606 .5447 .66229508 > KID5 -.2571164238 .12177700 -2.111 .0347 .49508197 > PHD .4732121566E-01 .73778012E-01 .641 .5213 3.1031093 > MENT .4233995753E-01 .90234066E-02 4.692 .0000 8.7672125 > Constant .1054079271E-01 .35024919E-01 .301 .7635 > FEM .6769194657 .31020676 2.182 .0291 .46010929 > MAR -.2515104394 .17217705 -1.461 .1441 .66229508 > KID5 .1321340645 .21811993 .606 .5447 .49508197 > PHD -.2571164238 .12177700 -2.111 .0347 3.1031093 > MENT .4732121566E-01 .73778012E-01 .641 .5213 8.7672125 > */ . . log close