. * prcounts: version 1.8.4 . * verified: 3/24/01 . which prcounts c:\stata6\ado\stbplus\p\prcounts.ado *! version 1.6.5 1/10/01 changes in labels for zi models . . use couart2.dta, clear (Academic Biochemists / S Long) . . ** CHECK RESULTS AGAINST XPOST . . *poisson . qui poisson art fem mar kid5 phd ment . prcounts pois, max(8) plot . su pois* Variable | Obs Mean Std. Dev. Min Max ---------+----------------------------------------------------- poisrate | 915 1.692896 .6685824 .8883344 9.627207 poispr0 | 915 .2092071 .0794247 .0000659 .4113403 poispr1 | 915 .3098447 .0634931 .0006345 .3678775 poispr2 | 915 .242096 .0311473 .0030544 .2706704 poispr3 | 915 .1346656 .0415861 .0098018 .2240418 poispr4 | 915 .0611696 .0383808 .0106732 .1951233 poispr5 | 915 .0249554 .0287183 .0018963 .1742638 poispr6 | 915 .0099346 .0201179 .0002808 .1603728 poispr7 | 915 .0041384 .0137756 .0000356 .1428533 poispr8 | 915 .001877 .0094055 3.96e-06 .1206255 poiscu0 | 915 .2092071 .0794247 .0000659 .4113403 poiscu1 | 915 .5190518 .1395755 .0007004 .7767481 poiscu2 | 915 .7611477 .1407294 .0037549 .9390502 poiscu3 | 915 .8958133 .1126566 .0135567 .9871097 poiscu4 | 915 .956983 .0824803 .0371477 .9977829 poiscu5 | 915 .9819384 .0589296 .0825709 .9996792 poiscu6 | 915 .991873 .0423403 .155454 .9999599 poiscu7 | 915 .9960114 .0310561 .2556911 .9999956 poiscu8 | 915 .9978884 .023188 .3763166 .9999995 poisprgt | 915 .0021116 .023188 4.77e-07 .6236834 poisval | 9 4 2.738613 0 8 poisobeq | 9 .1101396 .1153559 .0010929 .3005464 poispreq | 9 .1108765 .1174511 .001877 .3098447 poisoble | 9 .8150577 .2373893 .3005464 .9912568 poisprle | 9 .8122127 .2760109 .2092071 .9978884 . . list fem mar kid5 phd ment poisp* in 1 Observation 1 fem Men mar Married kid5 0 phd 2.52 ment 7 poispr0 .1414034 poispr1 .2766047 poispr2 .2705385 poispr3 .1764036 poispr4 .0862674 poispr5 .0337502 poispr6 .0110033 poispr7 .0030749 poispr8 .0007519 poisprgt .0002021 poispreq .2092071 poisprle .2092071 . . *0=0.141 (probs from XPost) . *1=0.277 . *2=0.271 . *3=0.176 . *4=0.086 . *5=0.034 . *6=0.011 . *7=0.003 . *8=0.001 . . *nbreg . qui nbreg art fem mar kid5 phd ment . prcounts nbre, max(8) plot . su nbre* Variable | Obs Mean Std. Dev. Min Max ---------+----------------------------------------------------- nbrerate | 915 1.704929 .7964249 .866364 12.14286 nbrepr0 | 915 .3035957 .0781645 .015145 .4801816 nbrepr1 | 915 .2722666 .037721 .0289042 .3030041 nbrepr2 | 915 .1800492 .0152423 .0397624 .1938989 nbrepr3 | 915 .1063388 .0206307 .0476373 .1439137 nbrepr4 | 915 .0598104 .0214007 .0194697 .1147368 nbrepr5 | 915 .0330554 .0184871 .0067502 .0955359 nbrepr6 | 915 .0182944 .0147986 .0022616 .0817175 nbrepr7 | 915 .0102735 .0115145 .0007389 .0716509 nbrepr8 | 915 .0059105 .0089094 .0002368 .0636456 nbrecu0 | 915 .3035957 .0781645 .015145 .4801816 nbrecu1 | 915 .5758622 .11392 .0440492 .7810719 nbrecu2 | 915 .7559114 .1142171 .0838115 .9169757 nbrecu3 | 915 .8622502 .0996432 .1314489 .9704345 nbrecu4 | 915 .9220607 .0821358 .1842904 .9899042 nbrecu5 | 915 .955116 .0663873 .2400891 .9966544 nbrecu6 | 915 .9734105 .053543 .2970283 .998916 nbrecu7 | 915 .9836839 .0434288 .3536867 .9996549 nbrecu8 | 915 .9895945 .0355319 .4089875 .9998918 nbreprgt | 915 .0104055 .0355319 .0001082 .5910125 nbreval | 9 4 2.738613 0 8 nbreobeq | 9 .1101396 .1153559 .0010929 .3005464 nbrepreq | 9 .1099549 .1153036 .0059105 .3035957 nbreoble | 9 .8150577 .2373893 .3005464 .9912568 nbreprle | 9 .8134983 .2344011 .3035957 .9895945 . . list fem mar kid5 phd ment nbrep* in 1 Observation 1 fem Men mar Married kid5 0 phd 2.52 ment 7 nbrepr0 .2499009 nbrepr1 .2591395 nbrepr2 .1936959 nbrepr3 .1260871 nbrepr4 .075993 nbrepr5 .0436011 nbrepr6 .0241746 nbrepr7 .0130704 nbrepr8 .0069315 nbreprgt .007406 nbrepreq .3035957 nbreprle .3035957 . . *0=0.250 (probs from XPost) . *1=0.259 . *2=0.194 . *3=0.126 . *4=0.076 . *5=0.044 . *6=0.024 . *7=0.013 . *8=0.007 . . *zip . qui zip art fem mar kid5 phd ment, inf(fem mar kid5 phd ment) . prcounts zip, max(8) plot . su zip* Variable | Obs Mean Std. Dev. Min Max ---------+----------------------------------------------------- ziprate | 915 1.693208 .6477186 .7477353 7.270344 zipall0 | 915 .1850598 .1060459 .0000161 .4318461 zippr0 | 915 .2985684 .1280144 .0007119 .5815108 zippr1 | 915 .2148732 .0356353 .0050586 .2906894 zippr2 | 915 .2103102 .030343 .0183893 .2613746 zippr3 | 915 .1426864 .0400169 .0445663 .2237694 zippr4 | 915 .0761093 .0380685 .0176175 .195136 zippr5 | 915 .0344642 .028516 .0044308 .1754144 zippr6 | 915 .0140484 .018923 .0009286 .1483315 zippr7 | 915 .0054281 .0119081 .0001668 .1482433 zippr8 | 915 .0020825 .0074178 .0000262 .1347246 zipcu0 | 915 .2985684 .1280144 .0007119 .5815108 zipcu1 | 915 .5134416 .1467138 .0057705 .7870966 zipcu2 | 915 .7237518 .1357939 .0241598 .9207867 zipcu3 | 915 .8664382 .104481 .0687261 .9768258 zipcu4 | 915 .9425476 .0715917 .1497305 .9944432 zipcu5 | 915 .9770118 .0462554 .2675184 .9988741 zipcu6 | 915 .9910602 .0292408 .410247 .9998027 zipcu7 | 915 .9964883 .0183864 .5584903 .9999695 zipcu8 | 915 .9985708 .0114556 .693215 .9999958 zipprgt | 915 .0014292 .0114556 4.23e-06 .306785 zipval | 9 4 2.738613 0 8 zipobeq | 9 .1101396 .1153559 .0010929 .3005464 zippreq | 9 .1109523 .1097276 .0020825 .2985684 zipoble | 9 .8150577 .2373893 .3005464 .9912568 zipprle | 9 .8119865 .2521604 .2985684 .9985708 . . list fem mar kid5 phd ment zipall0 ziprate in 1 Observation 1 fem Men mar Married kid5 0 phd 2.52 ment 7 zipall0 .13393 ziprate 2.037953 . . *CHECK ziprate=2.0379527 (-predict-) . *CHECK zipall0=.13393 . . *zinb . qui zinb art fem mar kid5 phd ment, inf(fem mar kid5 phd ment) . prcounts zinb, max(8) plot . su zinb* Variable | Obs Mean Std. Dev. Min Max ---------+----------------------------------------------------- zinbrate | 915 1.697666 .7165823 .468232 9.678316 zinball0 | 915 .0585168 .1215135 1.01e-30 .6024709 zinbpr0 | 915 .3119487 .1249968 .0169493 .7454953 zinbpr1 | 915 .2556639 .0418208 .0353106 .3098438 zinbpr2 | 915 .1807039 .0286549 .0506364 .2016944 zinbpr3 | 915 .1106089 .0262181 .0310986 .1509873 zinbpr4 | 915 .0630791 .0232455 .0121561 .1210988 zinbpr5 | 915 .0347421 .0188791 .0044736 .1011692 zinbpr6 | 915 .0188861 .0145442 .0015781 .0870389 zinbpr7 | 915 .010288 .0109449 .0005395 .0756084 zinbpr8 | 915 .0056812 .00819 .00018 .0657227 zinbcu0 | 915 .3119487 .1249968 .0169493 .7454953 zinbcu1 | 915 .5676126 .1256956 .0522599 .8773996 zinbcu2 | 915 .7483166 .1129196 .1028963 .9498874 zinbcu3 | 915 .8589255 .0933162 .1645512 .980986 zinbcu4 | 915 .9220046 .0737594 .2329503 .9931421 zinbcu5 | 915 .9567467 .0572472 .3043903 .9976156 zinbcu6 | 915 .9756327 .0443001 .375914 .9991937 zinbcu7 | 915 .9859208 .034455 .4453102 .9997332 zinbcu8 | 915 .991602 .0270324 .5110329 .9999132 zinbprgt | 915 .008398 .0270324 .0000868 .4889671 zinbval | 9 4 2.738613 0 8 zinbobeq | 9 .1101396 .1153559 .0010929 .3005464 zinbpreq | 9 .110178 .1140133 .0056812 .3119487 zinboble | 9 .8150577 .2373893 .3005464 .9912568 zinbprle | 9 .81319 .2340819 .3119487 .991602 . . capture drop checkvar . predict checkvar (option n assumed; predicted number of events) . di checkvar[1] 1.9852005 . capture drop checkvar . predict checkvar, p . di checkvar[1] .00034825 . . list fem mar kid5 phd ment zinball0 zinbrate in 1 Observation 1 fem Men mar Married kid5 0 phd 2.52 ment 7 zinball0 .0003482 zinbrate 1.985201 . . *CHECK zinbrate=1.9852005 (-predict-) . *CHECK zipall0=.00034825 . . . log close