
    M/PhZL                        d dl Zd dlmZmZ d dlmZ d dlmZm	Z	 d dl
mZ d dlmZ d dlmZmZmZmZmZmZmZmZ  G d d	          Z G d
 d          Z G d d          Z G d d          Z G d d          Z G d d          Z G d de          Z G d de          Z G d de          Z G d de          Z  G d d          Z! G d d          Z" G d  d!          Z#dS )"    N)assert_allcloseassert_equal)stats)poissonnbinom)PYTHON_IMPL_WASM)Bunch)genpoisson_ptruncatedpoissontruncatednegbin	zipoissonzinegbinzigenpoissonDiscretizedCountDiscretizedModelc                   ,    e Zd Zd Zd Zd Zd Zd ZdS )TestGenpoisson_pc                     t          j        dd          }t          j        dddd          }t          ||d           d S )N   r   V瞯<rtolr   pmfr
   r   selfpoisson_pmfgenpoisson_pmfs      m/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/statsmodels/distributions/tests/test_discrete.pytest_pmf_p1zTestGenpoisson_p.test_pmf_p1   D    k!Q''%)!Q155^%@@@@@@    c                     t          j        dd          }t          j        dddd          }t          ||d           d S )N   r   r   r   r   r   s      r   test_pmf_p2zTestGenpoisson_p.test_pmf_p2   r!   r"   c                     t          j        dd          }t          j        dddd          }t          ||d           d S )N
   r$   g}:   -q=r   r   )r   r   genpoisson_pmf_5s      r   test_pmf_p5zTestGenpoisson_p.test_pmf_p5#   sF    k"a(('+B5!<<%5EBBBBBBr"   c                     t          j        dd          }t          j        dddd          }t          ||d           d S )Nr(   r$   r   r   r   r   r   logpmfr
   r   r   s      r   test_logpmf_p1zTestGenpoisson_p.test_logpmf_p1(   D    nQ**%,Q1a88^%@@@@@@r"   c                     t          j        dd          }t          j        dddd          }t          ||d           d S )N   r   r   r$   r   r   r-   r   s      r   test_logpmf_p2zTestGenpoisson_p.test_logpmf_p2-   r0   r"   N)__name__
__module____qualname__r    r%   r+   r/   r3    r"   r   r   r      sp        A A A
A A A
C C C
A A A
A A A A Ar"   r   c                   *    e Zd ZdZd Zd Zd Zd ZdS )TestTruncatedPoisson-
    Test Truncated Poisson distribution
    c                     t          j        dd          t          j        dd          z  }t          j        ddd          }t	          ||d           d S Nr$   r   Hz>r   )r   r   sfr   r   r   r   tpoisson_pmfs      r   test_pmf_zeroz"TestTruncatedPoisson.test_pmf_zero7   sS    k!Q'''*Q*:*::'+Aq!44\======r"   c                     t          j        dd          t          j        t          j        dd                    z
  }t          j        ddd          }t          ||d           d S r<   r   r.   nplogr>   r   r   r   poisson_logpmftpoisson_logpmfs      r   test_logpmf_zeroz%TestTruncatedPoisson.test_logpmf_zero<   ]     1--wz!Q7G7G0H0HH*1!Q::dCCCCCCr"   c                     t          j        dd          dt          j        dd          z
  z  }t          j        ddd          }t	          ||d           d S )N   r2   r   r$   r=   r   )r   r   cdfr   r   r?   s      r   test_pmfzTestTruncatedPoisson.test_pmfA   sX    k!Q''1w{1a/@/@+@A'+Aq!44\======r"   c                     t          j        dd          t          j        t          j        dd                    z
  }t          j        ddd          }t          ||d           d S )NrL   r2   r$   r=   r   rC   rF   s      r   test_logpmfz TestTruncatedPoisson.test_logpmfF   rJ   r"   Nr4   r5   r6   __doc__rA   rI   rN   rP   r7   r"   r   r9   r9   3   sb         > > >
D D D
> > >
D D D D Dr"   r9   c                   >    e Zd Zd Zd Zd Zd Zd Zd Zd Z	d Z
d	S )
TestZIPoissonc                     t          j        dd          }t          j        ddd          }t          ||d           d S N   r$   r   r)   r   r   r   r   r   r   r   zipoisson_pmfs      r   rA   zTestZIPoisson.test_pmf_zeroN   A    k!Q''!aA..]??????r"   c                     t          j        dd          }t          j        ddd          }t          ||d           d S Nr(   r   r   r)   r   r   r.   r   r   r   rG   zipoisson_logpmfs      r   rI   zTestZIPoisson.test_logpmf_zeroS   sD     1--$+Aq!44(8uEEEEEEr"   c                     t          j        dd          }t          j        ddd          }t          ||dd           d S )Nr$   皙?皙?r   atolrX   rY   s      r   rN   zTestZIPoisson.test_pmfX   sC    k!Q''!aC00]DIIIIIIr"   c                     t          j        dd          }t          j        ddd          }t          ||dd           d S )N   rW   rb   rc   rd   r^   r_   s      r   rP   zTestZIPoisson.test_logpmf]   sF     1--$+Aq#66(8t$OOOOOOr"   c                     t          j        dd          }t          j        ddd          }t          ||d           d S rV   )r   rM   r   r   )r   poisson_cdfzipoisson_cdfs      r   test_cdf_zerozTestZIPoisson.test_cdf_zerob   r[   r"   c                     t          j        dd          }t          j        ddd          }t          ||d           d S r]   )r   ppfr   r   )r   poisson_ppfzipoisson_ppfs      r   test_ppf_zerozTestZIPoisson.test_ppf_zerog   r[   r"   c                 
   t          j        d          t          j        d          }}t          j        dd          }t          j        dd          }t	          ||d           t	          ||d           t          j        g d          }t          j        |          t          j        |          }}t          j        |d          }t          j        |d          }t	          ||d           t	          ||d           d S )N   r   绽|=r   r   r(   r'           )r   meanvarr   r   rD   array)r   poisson_meanpoisson_varzipoisson_meanzipoisson_varms         r   test_mean_varzTestZIPoisson.test_mean_varl   s    $+L$4$4gk"ook"A..!b!,,n5AAAA]????HZZZ  $+LOOW[^^k"1--!a--n5AAAA]??????r"   c                 J   t          j        dd          t          j        dd          }}t          j        ddd          }t          j        ddd          }t          j        ddd          }t          d|d           t          ||d           t          ||d           d S )Nr   rr   r$   r   rs   r   )r   momentr   r   )r   
poisson_m1
poisson_m2zip_m0zip_m1zip_m2s         r   test_momentszTestZIPoisson.test_momentsz   s    !(2!6!6q"8M8MJ
!!R++!!R++!!R++6....
F7777
F777777r"   N)r4   r5   r6   rA   rI   rN   rP   rk   rp   r~   r   r7   r"   r   rT   rT   L   s        @ @ @
F F F
J J J
P P P
@ @ @
@ @ @
@ @ @8 8 8 8 8r"   rT   c                   ,    e Zd Zd Zd Zd Zd Zd ZdS )TestZIGeneralizedPoissonc                     t          j        dddd          }t          j        ddddd          }t          ||d           d S )NrW   r$   r   r   r)   r   r
   r   r   r   r   gp_pmfzigp_pmfs      r   rA   z&TestZIGeneralizedPoisson.test_pmf_zero   sK    !!Q1--#Aq!Q22u555555r"   c                     t          j        dddd          }t          j        ddddd          }t          ||d           d S )Nrg   rW   r   r   r)   r   r
   r.   r   r   r   	gp_logpmfzigp_logpmfs      r   rI   z)TestZIGeneralizedPoisson.test_logpmf_zero   sK     '1a33	")!Q1a88	;U;;;;;;r"   c                     t          j        dddd          }t          j        ddddd          }t          ||dd           d S )NrW   r$   rb   rc   rd   r   r   s      r   rN   z!TestZIGeneralizedPoisson.test_pmf   sM    !!Q1--#Aq!Q44t$??????r"   c                     t          j        dddd          }t          j        ddddd          }t          ||dd           d S )Nr$   rW   r   rb   rc   rd   r   r   s      r   rP   z$TestZIGeneralizedPoisson.test_logpmf   sM     '1a33	")!Q1c::	;TEEEEEEr"   c                 &   t          j        g d          }t          j        |          t          j        |          }}t          j        |ddd          }t          j        |ddd          }t          ||d           t          ||d           d S )Nrt   r   r   ru   rs   r   )rD   rx   r   rv   rw   r   r   )r   r}   ry   rz   zigenpoisson_meanzigenpoisson_vars         r   r~   z&TestZIGeneralizedPoisson.test_mean_var   s     HZZZ  $+LOOW[^^k(-aAq99'+AsAq99&7eDDDD%5EBBBBBBr"   N)r4   r5   r6   rA   rI   rN   rP   r~   r7   r"   r   r   r      sj        6 6 6
< < <
@ @ @
F F F
C C C C Cr"   r   c                   \    e Zd Zd Zd Zd Zd Zd Zd Zd Z	d Z
d	 Zd
 Zd Zd Zd ZdS )	TestZiNBPc                     t          j        ddd          \  }}t          j        d||          }t          j        ddddd          }t	          ||dd           d S )N   rb   r$   d   {Gz?h㈵>rd   r   convert_paramsr   r   r   r   npnb_pmftnb_pmfs        r   r%   zTestZiNBP.test_pmf_p2   sa    &r3221CA&&,sBQ55d>>>>>>r"   c                     t          j        ddd          \  }}t          j        d||          }t          j        ddddd          }t	          ||dd           d S )Nr'   r   r$      r   rd   r   r   r   r.   r   r   r   r   	nb_logpmf
tnb_logpmfs        r   r3   zTestZiNBP.test_logpmf_p2   sa    &r1a001M#q!,,	_S"aD99
	:DtDDDDDDr"   c                     t          j        ddd          \  }}t          j        d||          }t          j        ddddd          }t	          ||dd           d S )Nr   rb   r$   r'   r   r)   rd   r   r   r   rM   r   r   r   r   
nbinom_cdfzinbinom_cdfs        r   test_cdf_p2zTestZiNBP.test_cdf_p2   sa    &r3221ZAq))
|BCA66
Lu5IIIIIIr"   c                     t          j        ddd          \  }}t          j        d||          }t          j        ddddd          }t	          ||dd           d S )Nr   r   r$   gHzG?r   r)   rd   r   r   r   rm   r   r   r   r   
nbinom_ppfzinbinom_ppfs        r   test_ppf_p2zTestZiNBP.test_ppf_p2   sa    &sAq111Za++
|D#q!Q77
Lu5IIIIIIr"   c                 0   t          j        ddd          \  }}t          j        ||          t          j        ||          }}t          j        dddd          }t          j        dddd          }t          ||d           t          ||d           d S Nrg   r   r$   r   rs   r   )r   r   r   rv   rw   r   )r   r   r   nbinom_mean
nbinom_var	zinb_meanzinb_vars          r   test_mran_var_p2zTestZiNBP.test_mran_var_p2   s    &q!Q//1"(+a"3"3VZ15E5EZM!Q1--	<1a++YU;;;;
H5999999r"   c                    t          j        ddd          \  }}t          j        d||          t          j        d||          }}t          j        ddddd          }t          j        ddddd          }t          j        ddddd          }t	          d|d           t	          ||d           t	          ||d           d S r   r   r   r   r   r   r   r   r   nb_m1nb_m2zinb_m0zinb_m1zinb_m2s           r   test_moments_p2zTestZiNBP.test_moments_p2       &q!Q//1}Q1--v}Q1/E/Eu/!Q1a00/!Q1a00/!Q1a007////wU3333wU333333r"   c                     t          j        ddd          \  }}t          j        d||          }t          j        ddddd          }t	          ||dz  d           d S )Nr   ?r$         ?r=   r   r   r   r   r   r   r   s        r   rN   zTestZiNBP.test_pmf   sc    &q#q111Jq!Q''	,q!S!S11	7Q;T::::::r"   c                     t          j        ddd          \  }}t          j        d||          }t          j        ddddd          }t	          ||dd           d S )Nr(   r   r$   g{Gzt?r   rd   r   r   s        r   rP   zTestZiNBP.test_logpmf   sa    &q!Q//1M!Q**	_Q1a77
	:DtDDDDDDr"   c                     t          j        ddd          \  }}t          j        d||          }t          j        ddddd          }t	          ||dd           d S )Nr   r   r$   r   r)   rd   r   r   s        r   test_cdfzTestZiNBP.test_cdf   sa    &q#q111Z1a((
|Aq#q!44
Lu5IIIIIIr"   c                     t          j        ddd          \  }}t          j        d||          }t          j        ddddd          }t	          ||dd           d S )Nr(   r   gQ?r   r)   rd   r   r   s        r   test_ppfzTestZiNBP.test_ppf   sa    &q!Q//1Za++
|D!Q155
Lu5IIIIIIr"   c                    t          j        ddd          \  }}d\  }}t          ||dd           t          ||dd           t          j        ddd	          \  }}d
\  }}t          ||dd           t          ||dd           d S )N   g333333?r$   )g?gq\?r)   rd   rg   g(\?r   )gD@gYY?)r   r   r   )r   r   r   n_truep_trues        r   test_convertzTestZiNBP.test_convert   s    &r4331@6E::::6E::::&q$221>6E::::6E::::::r"   c                 f   dt          j        g d          fD ]}t          j        |dd          \  }}t	          j        ||          t	          j        ||          }}t          j        |ddd          }t          j        |ddd          }t          ||d           t          ||d           d S )N	   rt   r   r   rs   r   )rD   rx   r   r   r   rv   rw   r   )r   r}   r   r   r   r   r   r   s           r   r~   zTestZiNBP.test_mean_var   s    RXjjj))* 	> 	>A*1a33DAq&,k!Q&7&7Aq9I9IK aAq11I|Aq!Q//HK????Ju=====	> 	>r"   c                    t          j        ddd          \  }}t          j        d||          t          j        d||          }}t          j        ddddd          }t          j        ddddd          }t          j        ddddd          }t	          d|d           t	          ||d           t	          ||d           d S )Nr   r   r$   r   rs   r   r   r   s           r   r   zTestZiNBP.test_moments   r   r"   N)r4   r5   r6   r%   r3   r   r   r   r   rN   rP   r   r   r   r~   r   r7   r"   r   r   r      s        ? ? ?E E EJ J JJ J J: : :4 4 4; ; ;E E EJ J JJ J J	; 	; 	;> > >4 4 4 4 4r"   r   c                       e Zd Zd Zd ZdS )CheckDiscretizedc                 Z    |                                 }|                    dd           |S )Nr   )tolistinsert)r   paramsargss      r   r   zCheckDiscretized.convert_params	  s(    }}Br"   c           	      N   | j         }| j        }| j        }| j        }| j        }| j        }t          j                            d           t          ||          }|j        |k    sJ t          j
        d          } |j        |g|R  }	 |j        |g|R  }
t          j        |
          }t          |	d t          |                   |d            |j        |g|R  }t          |d t          |
          dz
           |
dd          d            |j        |g|R  }t          ||	d            |j        |g|R  }t          ||d            |j        |g|R  }t          |d|z
  d           d} |j        |d|i}t          |          |k    sJ |                                dk    sJ t+          ||	          }|                    |
          }|                    |j        d          }	|                     |j                  }t          j         |j        t          j
        d          g|R             }t          |	|d           t          j         |j        t          j
        d          g|R            }t          |	|dd           t          j        |                    t8                              }t          |          }|dk    r9d}||dz
  xx         ||d                                          z  cc<   |d |         }|                    |j        d|          }	|	|dz
  xx         d|	d |                                         z
  z  cc<   t=          j        ||	d |         |z            }|j         dk    sJ |!                    |j                  }d}|                    |          }tD          r|                    t          j#                  }t          j        |          }|                    |j        d|          }	t          |          }|	|dz
  xx         d|	d |                                         z
  z  cc<   t=          j        ||	d |         |z            }|j         dk    sJ |$                    |                    t          j
        dd          dz                       }t          j%        g d          }tM          ||           t          j'        |                    t          j
        dd                    dz
  d          }	|$                    |	          }t          j
        dd          }tM          ||           |$                    |                    t          j
        d                              }t          j
        dd          }tM          ||           |(                    d|                    t          j
        dd          dz             z
            }t          j%        g d          }tM          ||           d S )N
 r(   gvIh%<=r   r   i  sizegMbP?distrstart_paramsprobswhich   r   rd   r'   )r   k_maxr   i  )r   r   gư>)g      g      ?g       @g      @g      @g      @r   ))d_offsetddistrparamgparamdshapesr   rD   randomseedr   arange_pmfrM   diffr   len_cdfr   r>   rvsrw   r   fitpredictr   r   bincountastypeintsumr   	chisquarepvalue	get_distrr   int32rm   rx   r   maximumisf)r   r   r   r   r   r   r   dpxir   cdf1p1rM   p2cdf2r>   nobsxxmodresr   freqktchi2dfrnobs_rvsr   qq1s                                r   
test_basiczCheckDiscretized.test_basic  s   =(
	vfh//yF""""Yq\\BGB    vz"&v&&&WT]])CGG)bu5555bgb"6"""Oc$ii!mO,d122hUCCCC RVB    AE****rvb"6"""c....RU2AG%0000RVV'$''2ww$vvxx%r,,,gg<g00KK
'K22""3:..gifi	"55556662E**** WZVZ	"5555662E6666{299S>>**IIr66AQKKK48<<>>)KKK8DKK
'K;;	!a%A"1"		O#aedl33|d""""
 mmCJ''gg8g$$
  	'**RX&&C{3KK
'KBBII	!a%A"1"		O#aeh&677|d"""" GGCGGBIb!,,t34455X44455QJswwryQ//0047;;GGAJJYr1QGGCGGBIaLL))**Yq!__QGGA	"a 0 04 788899X44455Qr"   N)r4   r5   r6   r   r  r7   r"   r   r   r     s7          
Z Z Z Z Zr"   r   c                   $    e Zd Zed             ZdS )TestDiscretizedGammac                 n    d| _         t          j        | _        d| _        d| _        d| _        d| _        d S )Nr   )r(   r   r   )r(   r   za, s)r   r   )r   r   gammar   r   r   r   r   clss    r   setup_classz TestDiscretizedGamma.setup_classm  s5    [
 


#r"   Nr4   r5   r6   classmethodr  r7   r"   r   r  r  k  s-        $ $ [$ $ $r"   r  c                   $    e Zd Zed             ZdS )TestDiscretizedExponentialc                 n    d| _         t          j        | _        d| _        d| _        d| _        d| _        d S )Nr   )r   r(   )r(   sr   )r   r   exponr   r   r   r   r   r  s    r   r  z&TestDiscretizedExponential.setup_classz  s5    [



r"   Nr   r7   r"   r   r#  r#  x  s-        ! ! [! ! !r"   r#  c                   $    e Zd Zed             ZdS )TestDiscretizedLomaxc                 n    d| _         t          j        | _        d| _        d| _        d| _        d| _        d S )Nr   )r$   r         ?)r$   r*  zc, s)r   r   )r   r   lomaxr   r   r   r   r   r  s    r   r  z TestDiscretizedLomax.setup_class  s5    [
 


%r"   Nr   r7   r"   r   r(  r(    s-        & & [& & &r"   r(  c                   $    e Zd Zed             ZdS )TestDiscretizedBurr12c                 n    d| _         t          j        | _        d| _        d| _        d| _        d| _        d S )Nr   )r$   r   r   r*  )r$   r   r*  zc, d, s)r   r   r   )r   r   burr12r   r   r   r   r   r  s    r   r  z!TestDiscretizedBurr12.setup_class  s5    \
#
 

(r"   Nr   r7   r"   r   r-  r-    s-        ) ) [) ) )r"   r-  c                       e Zd Zd ZdS )TestDiscretizedGammaExc           
      $   g d}t          j        t          j        d          |          }t          ddgddddd	g d
          }t	          t
          j                  }t          ||          }|                    ddg          }t          |          }t          |j        |j        d           t          |j        |j        d           t          |j        |j        d           t          |j        |j                   |                    |j        d          }|d t          |j                           }	|	dxx         d|	                                z
  z  cc<   t          |	|z  |j        d           t          t          j        |          |	|z                                  d           t          j        ||	|z  t          |j                            }
t          |
j        dd           t          |
j        dd           |                                 t           j                            d           |                                }t          |j        |d         d           t          |j        |d         d           d S )N).   L      r   r   r(   ggaO;5@gO=?g^Iogg;%8?r   gz6>W?g^Iw@)g=
ףp=G@gGznR@gzG;@g      @gQ?)r   llfchi2df_modelr   aicr   r   r   r   r   r   g~jtx?)re   r   r   r   gQ?rs   )ddofg.c؛C?r=   gYrL?r   rc   )rD   repeatr   r	   r   r   r  r   r   r   r   r   r6  r9  r   r8  r   r   r   r  	statisticr  summaryr   r   	bootstrapbse)r   r  yres1r  r  r  r  r   probs_truncres_chi2	res_bootss               r   test_allzTestDiscretizedGammaEx.test_all  sc   !!!IbillD))X&2224 4 4 ek**q+++ggAq6g**1vv
DKd;;;;55555555S\4=111CJg66,S__,-B1{0000d*DJTBBBBt{T'9&>&>&@&@uMMMM?4t);(+CJ9 9 9 	*JTBBBB$???? 	
	vMMOO	
IaLt<<<<1D999999r"   Nr4   r5   r6   rE  r7   r"   r   r1  r1    s#        ,: ,: ,: ,: ,:r"   r1  c                       e Zd Zd ZdS )TestGeometricc                    d}dt          j        d|z
            z  }t          j        |d          } t	          t          j                  |          }t          j        d          }|                    |          }|                    |          }t          ||d           |	                    |          }|	                    |          }	t          |	|d           |
                    |          }
|
                    |          }t          ||
d           |                    |          }|                    |          }t          ||           |                    |dz
            }|                    |dz
            }t          ||           |                    |dz             }|                    |dz             }t          ||           |                    d	          }|                    d	          }t          |d           |                    |
          }|                    |
          }t          ||           |                    |
dz
            }|                    |
dz
            }t          ||           |                    |
dz             }|                    |
dz             }t          ||           |                    d	          }|                    d	          }t          ||           |                    d          }|                    d          }t          |d           d S )
Ng333333?r   r   )locr2   rs   r   g:0yE>r   )rD   rE   r   geomr   r&  r   r   r   rM   r>   rm   r   r  )r   p_geomscale_dexpondgeodpgr  pmf1r   r	  rM   sf1r>   ppf1rm   isf1r  s                   r   rE  zTestGeometric.test_all  s   BF1V8,,,z&b)))+u{++L99Yq\\xx||ggbkkT....xx||ggbkkT....ggbkkVVBZZCe,,,,xx~~ggdmmS$xxt$$ggdTk""S$xxt$$ggdTk""S$xx{{ggajjS" xx}}ggcllS$xxd
##ggcDj!!S$xxd
##ggcDj!!S$xx{{ggajjS$xx{{ggajjS"r"   NrF  r7   r"   r   rH  rH    s#        - - - - -r"   rH  c                   *    e Zd ZdZd Zd Zd Zd ZdS )TestTruncatedNBPr:   c                     t          j        ddd          \  }}t          j        d||          t          j        d||          z  }t          j        ddddd          }t          ||d           d S )Nr(   rb   r$   r   r   r   r   )r   r   r   r   r>   r   r   s        r   rA   zTestTruncatedNBP.test_pmf_zero  ss    -aa881Aq!$$vyAq'9'99!%aCA66d333333r"   c                    t          j        ddd          \  }}t          j        d||          t	          j        t          j        d||                    z
  }t          j        ddddd          }t          ||dd           d S )Nr(   r   r$   r   r   rd   )r   r   r   r.   rD   rE   r>   r   r   s        r   rI   z!TestTruncatedNBP.test_logpmf_zero  s    -aA661M!Q**RVFIaA4F4F-G-GG	$+Aq!Q::
	:DtDDDDDDr"   c                 6   t          j        ddd          \  }}t          j        d||          t          j        d||          z  }t          j        ddddd          }t          ||d           t          j        ddddd          }t          |d           d S )Nr$   r   r2   r(   r=   r   r   )r   r   r   r   r>   r   r   r   s        r   rN   zTestTruncatedNBP.test_pmf  s    -aa881Jq!Q''&)Aq!*<*<<	!%aCA66	76666!%aCA66Wa     r"   c                 f   t          j        ddd          \  }}t          j        d||          t	          j        t          j        d||                    z
  }t          j        ddddd          }t          ||d           t          j        ddddd          }t	          j        |          sJ d S )Nr(   rb   r$   r2   r=   r   )	r   r   r   r.   rD   rE   r>   r   isneginfr   s        r   rP   zTestTruncatedNBP.test_logpmf  s    -aa881M!Q**RVFIaA4F4F-G-GG	$+Aq#q!<<
	:D9999$+Aq#q!<<
{:&&&&&&&r"   NrQ   r7   r"   r   rU  rU    s]         4 4 4E E E! ! !' ' ' ' 'r"   rU  )$numpyrD   numpy.testingr   r   scipyr   scipy.statsr   r   statsmodels.compat.pythonr   statsmodels.tools.toolsr	   "statsmodels.distributions.discreter
   r   r   r   r   r   r   r   r   r9   rT   r   r   r   r  r#  r(  r-  r1  rH  rU  r7   r"   r   <module>rb     s'       7 7 7 7 7 7 7 7       ' ' ' ' ' ' ' ' 6 6 6 6 6 6 ) ) ) ) ) )	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	A A A A A A A A:D D D D D D D D258 58 58 58 58 58 58 58pC C C C C C C C@`4 `4 `4 `4 `4 `4 `4 `4Fa a a a a a a aH
$ 
$ 
$ 
$ 
$+ 
$ 
$ 
$
! 
! 
! 
! 
!!1 
! 
! 
!
& 
& 
& 
& 
&+ 
& 
& 
&
) 
) 
) 
) 
), 
) 
) 
)/: /: /: /: /: /: /: /:d/ / / / / / / /d!' !' !' !' !' !' !' !' !' !'r"   