U
    Kvf                     @   s4   d Z ddlZddlmZmZ ddlmZ dd ZdS )a  
Test data from Lilliefors test for normality
An Analytic Approximation to the Distribution of Lilliefors's
Test Statistic for Normality
Author(s): Gerard E. Dallal and Leland WilkinsonSource: The American
Statistician, Vol. 40, No. 4 (Nov., 1986), pp. 294-296
    N)assert_almost_equalassert_allclose)	TableDistc                     sj  t ddddddgd d d } t dd	d
ddddddddddddddddddddgtt dd d!d"d#d$gd%dd&d'd(d)gd*d+d,d-d.d/gd0d1d2d3d4d5gd6d7d8d9d:d;gd<d=d0d>d?d@gdAdBdCdDd3dEgdFdGdHdIdJdKgdLdMdNdOd+dPgdQdRdSdTdUdVgdWdXdFdYdDdZgd[d\d]d^d_d`gdadbdcdddedfgdgdhd\didjdkgdldmdndRdTdogdpdqd[drdBdsgdtdudadvdwdxgdydzdtd{d]d|gd}d~ddtdndAgddddddgddddddgddddddgddddddggd d d d df d  t|  tddd tddd tddd ddd  krduks$n tddd  krDdksJn tddd  krjdkspn tddd  krdksn tddd  krdksn tddd  krdksn t fddt	t
D }t |dj}t |j }t|t | fddt	t
D }t |dj}t |  }t|t | tt t t fddt	t
D dj   fddt	t
D }t |djd }| d d d df k  stfddt	t
D }t |djd }| d d dd f k snt fddt	t
D }t |djd }|jd d k st fddt	t
D }t |djd }|jdd  k  st fddt	dt
D }t |djd }|jd d k sftd S )Ng?g333333?g?g?g{Gz?gMbP?               	   
                                       (   d   i  i  i/  iA  iZ  ix  i  i  i!  i?  iW  i  i  i  i  i)  iC  is  i     i  i  i0  i_  i        i	  i   iM  i        i  i=  im           i  i`              i#  iR              iE              i  i:           i1                 i(                 i              i           i           i
              i                       s   x               J   M   R   Y   h   z   %   '   )   -   4   =         #   *   g     @@gsh|??g'1Z?g(\?gRQ?gFx?c                    s*   g | ]"} | D ]} || qqS  prob.0icZcrit_lflfsizerg   J/tmp/pip-unpacked-wheel-2v6byqio/statsmodels/stats/tests/test_tabledist.py
<listcomp>=   s     
  z"test_tabledist.<locals>.<listcomp>c                    s(   g | ] } j D ]} || qqS rg   )alphacritrj   ro   rp   rg   rq   rr   C   s       c                    s(   g | ] } j D ]} || qqS rg   rs   crit3rj   ru   rg   rq   rr   J   s     c                    s4   g | ],} j d d d D ]} || qqS Nr   皙?rv   rj   ru   rg   rq   rr   M   s        c                    s4   g | ],} j d d d D ]} || qqS rx   rv   rj   ru   rg   rq   rr   R   s     c                    s6   g | ].} |d df D ]} |d | qqS Nr   g?rh   rj   rn   rg   rq   rr   W   s     c                    s6   g | ].} |d df D ]} |d | qqS )rz   Nry   rh   rj   rn   rg   rq   rr   \   s     c                    s6   g | ].} |d df D ]} || d qqS r{   rh   rj   rn   rg   rq   rr   b   s        )nparrayfloatr   r   ri   rt   rw   AssertionErrorrangelenZreshapeZn_alphaabsrs   r   Z
zeros_likeprintmaxall)rs   valsdeltarg   rn   rq   test_tabledist   s               &&&&&& 
"" r   )	__doc__Znumpyr}   Znumpy.testingr   r   Zstatsmodels.stats.tabledistr   r   rg   rg   rg   rq   <module>   s   