
    !gm              
       ~   d dl Zd dlmZ  edei ddddddd	d
dddddddd ddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?Z ej                  g d@dAdBdCdDdEej                  dFd dGdHdIdJdKdLej                  dFd d ej                  ej                  ej                  ej                  ej                  ej                  dFd dMdNdOdPdQdRej                  dFd dSdTdUdVdWdXej                  dFd dYdZd[d\d]d^ej                  dFd d_d`dadbdcddej                  dFd dedfdgdhdidjej                  dFd       j                  d
dk      Z	dlj                         Zg dmZ ej                  g dn      j                  d
d
      Zg dmZg dmZ edee	eeeeedoeZ edei ddddddd	ddddddddd ddddpddddqddrddsddtdduddvd d!d"d#d$d%dwdxd&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?Z ej                  g dydzd{d|d}d~ej                  dFd d ej                  ej                  ej                  ej                  ej                  ej                  dFd ddddddej                  dFd ddddddej                  dFd ddddddej                  dFd ddddddej                  dFd dddd ddej                  dFd       j                  ddk      Z	dlj                         Zg dZ ej                  g d      j                  dd      Zg dZg dZ edee	eeeeedoeZ edei ddddddd	d
dddddddd ddddddddrdddddej                  ddd d!d"d#d$d%d&d'd(dd*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?Z ej                  g ddddddej                  dFd ddddddej                  dFd d ej                  ej                  ej                  ej                  ej                  ej                  dFd ddddddej                  dFd ddddddej                  dFd ddddddej                  dFd ddddddej                  dFd ddÑdđdődƑdǑej                  dFd       j                  d
dk      Z	dlj                         Zg dmZ ej                  g dȢ      j                  d
d
      Zg dmZg dmZ edee	eeeeedoeZ edei ddrddddd	ddddddddd ddddɓddddddʓdd˓dej                  dd̓d d!d"d#d$d%dwdxd&d'd(dd*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?Z ej                  g d͑dΑdϑdБdёdґej                  dFd d ej                  ej                  ej                  ej                  ej                  ej                  dFd dӑdԑdՑd֑dבdؑej                  dFd dّdڑdۑdܑdݑdޑej                  dFd dߑdddddej                  dFd dߑdddddej                  dFd dddd ddej                  dFd       j                  ddk      Z	dlj                         Zg dZ ej                  g d      j                  dd      Zg dZg dZ edee	eeeeedoeZ edei ddddddd	d
dddddddd ddddddddrdddddej                  ddd d!d"d#d$d%d&d'd(dd*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?Z ej                  g ddddddej                  dFd ddddddej                  dFd d ej                  ej                  ej                  ej                  ej                  ej                  dFd dddddd ej                  dFd ddddddej                  dFd ddd	d
ddej                  dFd ddddddej                  dFd ddddddej                  dFd       j                  d
dk      Z	dlj                         Zg dmZ ej                  g d      j                  d
d
      Zg dmZg dmZ edee	eeeeedoeZ edei ddrddddd	ddddddddd dddddddddddddej                  ddd d!d"d#d$d%dwdxd&d'd(dd*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?Z ej                  g dddd d!d"ej                  dFd d ej                  ej                  ej                  ej                  ej                  ej                  dFd d#d$d%d&d'd(ej                  dFd d)d*d+d,d-d.ej                  dFd d/d0d1d2d3d4ej                  dFd d5d0d6d7d8d9ej                  dFd d:d;d<d d=d>ej                  dFd       j                  ddk      Z	dlj                         Zg dZ ej                  g d?      j                  dd      Zg dZg dZ edee	eeeeedoeZ edei ddrddddd	ddddddddd dddddddddd@ddAdej                  ddBd d!d"d#d$d%dwdxd&d'd(dd*d+dCdDd,dEd.d/d0d1d2dFd4d5d6d7d8d9d:d;d<d=d>d?Z ej                  g ddGdHdIdJdKej                  dFd d ej                  ej                  ej                  ej                  ej                  ej                  dFd d#dLdMdNdOdPej                  dFd d)dQdRdSdTdUej                  dFd d/dVdWdXdYdZej                  dFd d5dVd[d\d]d^ej                  dFd d:d_d`dadbdcej                  dFd       j                  ddk      Z	dlj                         Zg dZ ej                  g dd      j                  dd      Zg dZg dZ edee	eeeeedoeZy(f      N)ParamsTableTestBunchrank   N
   ic   k   k_eq   k_dv	convergedrc	k_autoCnsllg5@
k_eq_modelll_0g		udf_m   chi2g`V]Ã@pgue$r2_pg*%?cmdlinez'poisson deaths lnpyears smokes i.agecatcmdpoissonpredictpoisso_p	estat_cmdpoisson_estatgofpoiss_gchi2typeLRopt	moptimizevceoimtitlezPoisson regressionuserpoiss_lfcrittypezlog likelihood	ml_methode2singularHmethodzm-marquardt	techniquenrwhichmaxdepvardeaths
propertieszb Vgri7?gsY?gmNۮ?gY}?ggmD?gS 1\?ghYy0?g)Ճb+?goڸ?gyzP?g%gQ=M큔@gAc{O?g+Jyd?g8vp@gh6k>ghƎ?gÒ{W?g|6!A@g^(0>`?gv@g'd>g.?g  =
@gieD@gcTd?g
M%@gIX5o?gaݸ&N?g@@giBgH@gK?gbf	?gFԪ?gg	YWg@gPRg>C@gv$gI?g4c?0g[@	   z!b se z pvalue ll ul df crit eform)lnpyearssmokes	1b.agecat2.agecat3.agecat4.agecat5.agecat_cons)@g|?疺¾r   Ծ|?N?/?.?*kPG?Y)Ir@   g	E!q?r   `Pߊÿa̹0ڿQl\iLnMaR@r   r   r   r   r   r   r   r   rA   rF   r   gLi˰?Y̶?cjO=?!K?mwrB   rG   r   rK   gwAI)?[%o?Մ`97?/ʧrC   rH   r   rL   rO   gW6?g88uh?籛[LUrD   rI   r   rM   rP   rR   g'>@3:!rE   rJ   r   rN   rQ   rS   rT   gD}~*B@)params_tableparams_table_colnamesparams_table_rownamescovcov_colnamescov_rownamesgx@gc%	~   gtTz׌@gdoK!g?z0poisson deaths smokes i.agecat, exposure(pyears)offsetz
ln(pyears)g#=?g{|?gM<j
@g$uxO?gnjq?g
Ǥ[?g,r!~?gf%?g5~<ul@g)=g=BL3?g.O?gWNg!@g
`?g@N+,@gs2s6g(W#@gE@g N)
@g7?g@Ve!2@gE10g@/@gt[#@g%+̙@gC#?gI?3@gwq ].g*:
@g ?N@g*cgE?gz;Dgu  g%B,)r9   r:   r;   r<   r=   r>   r?   )1gjN2?r   `0*,D>M>T/L	
ӫ=	5݂r   r   r   r   r   r   r   r]   r   gJjdL}?^*?YU*?P&yt$?w+r^   r   rb   gZmH?=d} 0?je'?獇[ar_   r   rc   rf   gr4|?tC'?
0BLgr`   r   rd   rg   ri   gZ?_Hra   r   re   rh   rj   rk   gӢ?   gEƻ@gDlnI|@gi	<,z7poisson deaths lnpyears smokes i.agecat, constraints(1)Waldg> [?g6i?g@gWu>grC?gBy0?g""?g?g7Ѫs@g+co>gݪ.$?gZ̼^?g tn?g*lW?gA7Q=@gb=g3.?gSw9,?gCj@gޗG?gtt&@g_!1H9g F?g$hR@g3P>@g-?g"@gq VF<g4S@gGg
@g0Nqzg8?gqER5	gەz[?g$ǒgHHӦ)@gF:i?ѣj:r   㶒	l=r?0?Eo?rq   `mRĿrn   g],p?r   
|c 撿`	$/ru   S#?r   r   r   r   r   r   r   r   ro   rs   r   gqD؇?+?矾Q<?rx   $F=rp   rt   r   rw   gH?X?rz   Ϩ4BJ¿rq   ru   r   rx   rz   J]~?r|   $,Կrq   ru   r   rx   rz   r|   r|   r}   rr   rv   r   ry   r{   r}   r}   g&" -?g>A9Cg~}@(@gS.V#z?poisson deaths smokes i.agecat, exposure(pyears) constraints(1)gDFa?gC!fx?gx[tS	@gV?g(G?g˴"?g$?g_x/?gr8Xo@g=g3붢?g =\Q?gzN@g-w?g>,YM,@g`br6g{%@gss x@g6p@g@?g,DM3@g_,.g(A8@gc"@gg<#?g7Dg`QuQ g Z#)1g\{38?r   y0wLD.#)Mt1Hr   ̾̂r   r   r   r   r   r   r   r~   r   g<[}?=r*?+V(?r   ;r   r   r   geએH?!Cɘ-?r   EӼ1`r   r   r   r   .-=/?r   :pMxr   r   r   r   r   r   r   r   r   r   r   r   r   g6΢?gG~;@g(>du@gq5D1(z7poisson deaths lnpyears smokes i.agecat, constraints(2)g܎?g A?g*l1+!@gf ,<g8?g[K*?gsq#k?g@!?g$B]?g%?gPb_a>Կgąx?g{M&q?gi
h?gZ @gMm=gSRO7?gqK?gif~@g-UQ1m?gnuq*@gw7gN}@g~2{+	@gDaʲ@g}SM?gZ ص'@g5}P9gW#@g\܎@g"0eY@g)l+@g!YGe7gW#@g\܎@gjh"g%b?gܿ#g_ -=g5'g%ԩ)@g>?〰6r   íus?c43?E꽆7?r   .*:Ŀr   g`p?r   fjT~=r   Dn?r   r   r   r   r   r   r   r   r   r   r   gMA+?C(נ?AA8?r   Q,沿r   r   r   r   gf7ӥ?IP?r   '8Y¿r   r   r   r   r   XRI?r   8R*	տr   r   r   r   r   r   r   r   r   r   r   r   r   r   r   gG(&?gIDAgUS@gt2Ir'z?poisson deaths smokes i.agecat, exposure(pyears) constraints(2)g*@?g|usx?gvk9
@gΚ}5J?gŜ!?g6\C>?gT[X?g)˥ ?g^k@gЀv =g _?g[=?gPJ@g)R?gtK!,@gti?sEt6g="@g1d@gPTQ
@gtZ?gIQۊK2@gkn$?b0gY@gG^#@gPTQ@g"+]5@g
c+gY@g\@g/g'v6b/?guDgbv.J gN0)1g:?r   N
VD^cMKMkLFr   PJk܂r   r   r   r   r   r   r   r   r   gmA~C}?*?--'?r   -r   r   r   g>YH?QqCF,?r   5br   r   r   r   qvB?r   ,r   r   r   r   r   r   r   r   r   r   r   r   r   g432Ԣ?gi2,4@g,< /&zJpoisson deaths smokes i.agecat,exposure(pyears) constraints(2) vce(robust)vcetypeRobustrobustzlog pseudolikelihoodge?gQ@g)WIa?gAˣ?g\p?g$]?gr@g8b'=goW?g9D?gN־?g+3VS%)@g䘈㤏8gdɑ@gp=K@gn:O?gXR-@gC/|G5g@g&@h@gl\"1@gTk2g
@g`@g^?gYp?g<Yd3g g)z<)1g[^&?r   ۶[vd?tϧMp?E's?r   +r   r   r   r   r   r   r   r   r   g2?甑z?@)X|?r   Dr   r   r   g' KZ?]?r   v2r   r   r   r   -0d?r   OMr   r   r   r   r   r   r   r   r   r   r   r   r   g.|? )numpynpstatsmodels.tools.testingr   dictestarraynanreshaperU   splitrV   rW   rX   rY   rZ   results_noexposure_noconstraintresults_exposure_noconstraintresults_noexposure_constraintresults_exposure_constraintresults_noexposure_constraint2results_exposure_constraint2#results_exposure_constraint2_robustr       {/var/www/dash_apps/app1/venv/lib/python3.12/site-packages/statsmodels/discrete/tests/results/results_poisson_constrained.py<module>r      ss!    :
 """ " 	"
 " " " " " !" " #" " "" "  "!"" =#"$ %"& '"( %)"* +", -". /"0 1"2 &3"4 5"6 %7"8 9": );"< ="> ?"@ A"B C"H rxx %(7:I&(*1@  	 * -< 	 &	 )8	 :<	
 
 %&
 78
 :< FF FF FF FF FF  ./ 2A  & )8 <J  VV . ?@  % (7 :I  & (*vv 1@  	 * -<  & )8 :<  %& )8 <J  & (8 ;J  FF!  !  ./!  2A!" #" &#" )8#" :J#$ %$ FF%$ -%$ >?% & wq!}' * <AAC F  bhh  " wq!}# &== #7 #)&;&;))# #   ### # 	#
 # # # # #  # # ## # "# #  "!#" F##$ %#& '#( %)#* +#, -#. /#0 1#2 3#4 &5#6 7#8 %9#: ;#< )=#> ?#@ A#B C#D E#J rxx 00%0(70:I00%0')vv00?0 0 0 0 "$0 FF	0 FF	0 FF	0 FF	0
 0
 %&0
 )80
 ;J0 0 %0 (70 :I0 FF0 0 ./0 2A0 0 &0 )80 ;J0 0 VV0 .0 ?@0 0 %0 (70 :I0 0 &0 (*vv0 1@0 0 	0 *0 .<0 0 &0 )80 :<0 0 %&0 )80 ;J0 0 &'0 )90 ;K0  FF!0  !0  ./!0 1  29A! $ <AAC :  bhh   GAqM 11 !5 !)&;&;))! !   !!! ! 	!
 ! ! ! ! ! !! ! ! "! ! !  M!!" #!$ %!& %'!( )!* +!, -!. /!0 &1!2 3!4 %5!6 7!8 )9!: ;!< =!> ?!@ A!F rxx %(7:I%')vv0?  	 * -< 	 &	 )8	 :<	
 
 %&
 78
 :< FF FF FF FF FF  ./ 2A  & )8 ;J  FF - >?  & ,8 ;J  $ &(VV />   * -<  & )8 :<  %& +8 ;J  & )8 ;J  FF!  !  ./!  1A!" #" &#" )8#" :J#$ %$ VV%$ .%$ ?@% & wq!}' * <AAC F  bhh  " wq!}# &== !5 !)&;&;))! !   ### # 	#
 # # # # # !# # # "# # # 6!#$ %#& '#( %)#* +#, -#. /#0 1#2 3#4 &5#6 7#8 %9#: ;#< )=#> ?#@ A#B C#D E#J rxx 00%0)70:I00%0')vv00?0 0 0 0 "$0 FF	0 FF	0 FF	0 FF	0
 0
 %&0
 )80
 ;J0 0 &0 )80 ;J0 FF0 0 ./0 2A0 0 &0 )80 ;J0 0 VV0 .0 ?@0 0 &0 )80 ;J0 0 &0 (*vv0 1@0 0 	0 *0 -<0 0 &0 )80 :<0 0 %&0 (80 ;J0 0 &'0 )90 ;K0  FF!0  !0  ./!0 1  29A! $ <AAC :  bhh   gam 11 3 )&;&;))    !!! ! 	!
 ! ! ! ! ! !! ! ! !! ! !  M!!" #!$ %!& %'!( )!* +!, -!. /!0 &1!2 3!4 %5!6 7!8 )9!: ;!< =!> ?!@ A!F rxx %(7:I%')vv0?  	 * -< 	 &	 )8	 :<	
 
 %&
 78
 :< FF FF FF FF FF  ./ 2A  & )8 ;J  VV . ?@  & )8 ;J  & (*vv 1@  
 * -<  & )8 :<  %& *8 ;J  & )8 ;J  FF!  !  ./!  1A!" #" &#" )8#" :J#$ %$ ff%$  /%$ @A% & wq!}' * <AAC F  bhh  " wq!}# &== "6 ")&;&;))" "   ### # 	#
 # # # # #  # # # "# # # 6!#$ %#& '#( %)#* +#, -#. /#0 1#2 3#4 &5#6 7#8 %9#: ;#< )=#> ?#@ A#B C#D E#J rxx 00%0(70:I00%0')vv00?0 0 0 0 "$0 FF	0 FF	0 FF	0 FF	0
 0
 %&0
 )80
 ;J0 0 &0 )80 ;J0 FF0 0 ./0 2A0 0 &0 )80 ;J0 0 FF0 -0 >?0 0 &0 )80 ;J0 0 &0 (*vv0 1@0 0 	0 *0 -<0 0 &0 )80 :<0 0 %&0 (80 ;J0 0 &'0 )90 ;K0  FF!0  !0  ./!0 1  29A! $ <AAC :  bhh   gam 11  4  )&;&;))      $$$ $ 	$
 $ $ $ $ $  $ $ $ "$ $ $ B!$$ %$& '$( %)$* +$, -$. /$0 1$2 3$4 5$6 &7$8 9$: +;$< =$> )?$@ A$B C$D E$F G$L rxx 00%0(70:I00$0&(ff0/>0 0 0 0 "$0 FF	0 FF	0 FF	0 FF	0
 0
 %&0
 )80
 ;J0 0 &0 )80 ;J0 FF0 0 ./0 2A0 0 &0 )80 <J0 0 VV0 .0 ?@0 0 &0 )80 ;J0 0 &0 (*vv0 1@0 0 	0 *0 -<0 0 &0 )80 :<0 0 %&0 (80 <J0 0 '0 )90 ;K0  FF!0  !0  ./!0 1  29A! $ <AAC :  bhh   gam 11 '; ')&;&;))' ' #r   