U
    Ãmœd®  ã                   @   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d€d�d‚dƒd„d…gd†d‡dˆd‰dŠd‹gdŒd�dŽd�d�d‘gdd’d“dd”d•gg¡d d …d d d…f d– ‰ t| ˆˆ ƒ‰tˆ d—d¡dƒ tˆ dd¡d—ƒ tˆ dd¡d—ƒ d˜ˆ d™d¡  k�rduk�s$n t‚d˜ˆ d™d¡  k�rDd—k�sJn t‚d˜ˆ dšd¡  k�rjd—k�spn t‚d˜ˆ dšd¡  k�r�d—k�s–n t‚d˜ˆ d›d¡  k�r¶d—k�s¼n t‚d˜ˆ d›d¡  k�rÜd—k�sân t‚‡ ‡‡fdœd�„t	t
ˆƒƒD ƒ}t  |¡ dˆj¡}t  |¡ˆj }t|t  |¡ƒ ‡‡fdžd�„t	t
ˆƒƒD ƒ}t  |¡ dˆj¡}t  |ˆ  ¡}t|t  |¡ƒ tt  t  t  ‡‡fdŸd�„t	t
ˆƒƒD ƒ¡ dˆj¡ˆ  ¡¡ƒ ‡‡fd d�„t	t
ˆƒƒD ƒ}t  |¡ dˆjd¡ ¡}|ˆ d d …d d…f k  ¡ �st‚‡‡fd¢d�„t	t
ˆƒƒD ƒ}t  |¡ dˆjd¡ ¡}|ˆ d d …d¡d …f k ¡ �snt‚‡ ‡‡fd£d�„t	t
ˆƒƒD ƒ}t  |¡ dˆjd¡ ¡}|ˆjd d… k ¡ �sÀt‚‡ ‡‡fd¤d�„t	t
ˆƒƒD ƒ}t  |¡ dˆjd¡ ¡}|ˆjd¡d … k  ¡ �st‚‡ ‡‡fd¥d�„t	d¦t
ˆƒƒD ƒ}t  |¡ dˆjd¡ ¡}|ˆjd d… k ¡ �sft‚d S )§Ngš™™™™™É?g333333Ã?gš™™™™™¹?gš™™™™™©?g{®Gáz„?gü©ñÒMbP?éÿÿÿÿé   é   é   é   é   é	   é
   é   é   é   é   é   é   é   é   é   é   é   é   é(   é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È?gF¶óýÔxÉ?c                    s*   g | ]"}ˆ | D ]}ˆ  |ˆ| ¡‘qqS © ©Úprob©Ú.0ÚiÚc©Zcrit_lfÚlfÚsizerg   ú_/home/sam/Atlas/atlas_env/lib/python3.8/site-packages/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__Únumpyr}   Znumpy.testingr   r   Zstatsmodels.stats.tabledistr   r‰   rg   rg   rg   rq   Ú<module>   s   