
    VfdJ                        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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)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      Mlib/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.   r/   r!   N)__name__
__module____qualname__r   r$   r*   r.   r2    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_zero8   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   r1   r   r#   r<   r   )r   r   cdfr
   r   r>   s      r   test_pmfzTestTruncatedPoisson.test_pmfB   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 )NrK   r1   r#   r<   r   rB   rE   s      r   test_logpmfz TestTruncatedPoisson.test_logpmfG   rI   r!   Nr3   r4   r5   __doc__r@   rH   rM   rO   r6   r!   r   r8   r8   4   sb         > > >
D D D
> > >
D D D D Dr!   r8   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   r@   zTestZIPoisson.test_pmf_zeroO   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   rF   zipoisson_logpmfs      r   rH   zTestZIPoisson.test_logpmf_zeroT   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   atolrW   rX   s      r   rM   zTestZIPoisson.test_pmfY   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   rV   ra   rb   rc   r]   r^   s      r   rO   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 rU   )r   rL   r   r   )r   poisson_cdfzipoisson_cdfs      r   test_cdf_zerozTestZIPoisson.test_cdf_zeroc   rZ   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_zeroh   rZ   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   rC   array)r   poisson_meanpoisson_varzipoisson_meanzipoisson_varms         r   test_mean_varzTestZIPoisson.test_mean_varm   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   rq   r#   r   rr   r   )r   momentr   r   )r   
poisson_m1
poisson_m2zip_m0zip_m1zip_m2s         r   test_momentszTestZIPoisson.test_moments{   s    !(2!6!6q"8M8MJ
!!R++!!R++!!R++6....
F7777
F777777r!   N)r3   r4   r5   r@   rH   rM   rO   rj   ro   r}   r   r6   r!   r   rS   rS   M   s        @ @ @
F F F
J J J
P P P
@ @ @
@ @ @
@ @ @8 8 8 8 8r!   rS   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 )NrV   r#   r   r   r(   r   r	   r   r   r   r   gp_pmfzigp_pmfs      r   r@   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 )Nrf   rV   r   r   r(   r   r	   r-   r   r   r   	gp_logpmfzigp_logpmfs      r   rH   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 )NrV   r#   ra   rb   rc   r   r   s      r   rM   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#   rV   r   ra   rb   rc   r   r   s      r   rO   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 )Nrs   r   r   rt   rr   r   )rC   rw   r   ru   rv   r   r   )r   r|   rx   ry   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)r3   r4   r5   r@   rH   rM   rO   r}   r6   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   ra   r#   d   {Gz?h㈵>rc   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   rc   r   r   r   r-   r   r   r   r   	nb_logpmf
tnb_logpmfs        r   r2   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   ra   r#   r&   r   r(   rc   r   r   r   rL   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(   rc   r   r   r   rl   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 Nrf   r   r#   r   rr   r   )r   r   r   ru   rv   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   rM   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   rc   r   r   s        r   rO   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(   rc   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(   rc   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(   rc   rf   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	   rs   r   r   rr   r   )rC   rw   r   r   r   ru   rv   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   rr   r   r   r   s           r   r   zTestZiNBP.test_moments   r   r!   N)r3   r4   r5   r$   r2   r   r   r   r   rM   rO   r   r   r   r}   r   r6   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           	         | 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}|                    |          }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          }tI          ||           t          j%        |                    t          j
        dd                    dz
  d          }	|"                    |	          }t          j
        dd          }tI          ||           |"                    |                    t          j
        d                              }t          j
        dd          }tI          ||           |&                    d|                    t          j
        dd          dz             z
            }t          j#        g d          }tI          ||           d S )N
 r'   gvIh%<=r   r   i  sizegMbP?distrstart_paramsprobswhich   r   rc   r&   )r   k_maxr   i  )r   r   gư>)g      g      ?g       @g      @g      @g      @r   )'d_offsetddistrparamgparamdshapesr   rC   randomseedr   arange_pmfrL   diffr   len_cdfr   r=   rvsrv   r   fitpredictr   r   bincountastypeintsumr   	chisquarepvalue	get_distrrl   rw   r   maximumisf)r   r   r   r   r   r   r   dpxir   cdf1p1rL   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$${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)r3   r4   r5   r   r  r6   r!   r   r   r     s7          
T T T T T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_classh  s5    [
 


#r!   Nr3   r4   r5   classmethodr  r6   r!   r   r  r  f  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_classu  s5    [



r!   Nr  r6   r!   r   r!  r!  s  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  r6   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  r6   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?)rd   r   r   r   gQ?rr   )ddofg.c؛C?r<   gYrL?r   rb   )rC   repeatr   r   r   r   r  r   r   r   r   r   r4  r7  r   r6  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r3   r4   r5   rC  r6   r!   r   r/  r/    s#        ,: ,: ,: ,: ,:r!   r/  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   )locr1   rr   r   g:0yE>r   )rC   rD   r   geomr   r$  r   r   r   rL   r=   rl   r   r  )r   p_geomscale_dexpondgeodpgr  pmf1r   r  rL   sf1r=   ppf1rl   isf1r  s                   r   rC  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!   NrD  r6   r!   r   rF  rF    s#        - - - - -r!   rF  c                   *    e Zd ZdZd Zd Zd Zd ZdS )TestTruncatedNBPr9   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'   ra   r#   r   r   r   r   )r   r   r   r   r=   r   r   s        r   r@   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   rc   )r   r   r   r-   rC   rD   r=   r   r   s        r   rH   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   r1   r'   r<   r   r   )r   r   r   r   r=   r   r   r   s        r   rM   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'   ra   r#   r1   r<   r   )	r   r   r   r-   rC   rD   r=   r   isneginfr   s        r   rO   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!   NrP   r6   r!   r   rS  rS    s]         4 4 4E E E! ! !' ' ' ' 'r!   rS  )"numpyrC   numpy.testingr   r   scipyr   scipy.statsr   r   statsmodels.tools.toolsr   "statsmodels.distributions.discreter	   r
   r   r   r   r   r   r   r   r8   rS   r   r   r   r  r!  r&  r+  r/  rF  rS  r6   r!   r   <module>r_     s       7 7 7 7 7 7 7 7       ' ' ' ' ' ' ' ' ) ) ) ) ) )	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	 	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[ [ [ [ [ [ [ [|
$ 
$ 
$ 
$ 
$+ 
$ 
$ 
$
! 
! 
! 
! 
!!1 
! 
! 
!
& 
& 
& 
& 
&+ 
& 
& 
&
) 
) 
) 
) 
), 
) 
) 
)/: /: /: /: /: /: /: /:d/ / / / / / / /d!' !' !' !' !' !' !' !' !' !'r!   