
    M/Ph                     4    d Z ddlZddlmZmZ ddlmZ 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                  f   t          j        g d          d d d         } t          j        g dt                    t          j        g dg dg dg dg dg d	g d
g dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg          d d d d df         dz  t          |           t	                              dd          d           t	                              dd          d           t	                              dd          d           d                    d d          cxk    rd!k    sn J d                    d d          cxk    rdk    sn J d                    d"d          cxk    rdk    sn J d                    d"d          cxk    rdk    sn J d                    d#d          cxk    rdk    sn J d                    d#d          cxk    rdk    sn J fd$t          t                              D             }t          j        |          
                    dj                  }t          j        |          j        z
  }t          |t          j        |                     fd%t          t                              D             }t          j        |          
                    dj                  }t          j        |z
            }t          |t          j        |                     t!          t          j        t          j        t          j        fd&t          t                              D                       
                    dj                  z
                                 fd't          t                              D             }t          j        |          
                    dj        d(z
            }|d d d df         k                                     sJ fd)t          t                              D             }t          j        |          
                    dj        d(z
            }|d d d(d f         k                                    sJ fd*t          t                              D             }t          j        |          
                    dj        d(z
            }|j        d d         k                                    sJ fd+t          t                              D             }t          j        |          
                    dj        d(z
            }|j        d(d          k                                     sJ fd,t          d-t                              D             }t          j        |          
                    dj        d(z
            }|j        d d         k                                    sJ d S ).N)g?333333?g?g?g{Gz?gMbP?)               	   
                                       (   d   i  i  )/  iA  iZ  ix  i  i  )i!  r   i?  iW  i  i  )i    i)  iC  is  i  )   i  i  0  i_  i  )      i	  i   iM  i  )      r    i  i=  im  )           r!   i`  )            i#  iR  )            r   iE  )            i  i:  )      r*      r)   i1  )               i(  )               i  )      r:         i  )         r3   r5   i  )      r9      r'   i
  )      r>         i  )      rM      r;      )         rM   rI   r&   )s   x               )J   M   R   Y   h   z   )%   '   )   -   4   =   )r         r   #   *   g     @@gsh|??r   r   g'1Z?g(\?rN   gRQ?gFx?c                 ^    g | ])}|         D ]}                     ||                   *S  prob.0iccrit_lflfsizes      f/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/statsmodels/stats/tests/test_tabledist.py
<listcomp>z"test_tabledist.<locals>.<listcomp>=   s<    NNNA71:NNaBGGAtAwNNNN    c                 \    g | ](}j         D ]}                    ||                   )S ro   )alphacritrs   rt   ru   rw   rx   s      ry   rz   z"test_tabledist.<locals>.<listcomp>C   s:    LLLA28LLaBGGAtAwLLLLr{   c                 \    g | ](}j         D ]}                    ||                   )S ro   r}   crit3r   s      ry   rz   z"test_tabledist.<locals>.<listcomp>J   sM     	 	 	!		 	A!T!W		 	 	 	 	r{   c                 r    g | ]3}j         d d         dz  D ]}                    ||                   4S Nr   皙?r   r   s      ry   rz   z"test_tabledist.<locals>.<listcomp>M   s^     * * *QXcrc]S(* * HHQQ   * * * *r{      c                 r    g | ]3}j         d d         dz  D ]}                    ||                   4S r   r   r   s      ry   rz   z"test_tabledist.<locals>.<listcomp>R   sZ     ! ! !QHSbSMC! !qBHHQQ   ! ! ! !r{   c                 l    g | ]0}|d df         D ]!}                     |dz  |                   "1S Nr   g?rp   rr   s      ry   rz   z"test_tabledist.<locals>.<listcomp>W   s_     & & &!QV_& & GGAGT!W%% & & & &r{   c                 l    g | ]0}|d df         D ]!}                     |dz  |                   "1S )r   Nr   rp   rr   s      ry   rz   z"test_tabledist.<locals>.<listcomp>\   s[       !AqrrEN ABGGAGT!W%%    r{   c                 l    g | ]0}|d df         D ]!}                     ||         dz            "1S r   rp   rr   s      ry   rz   z"test_tabledist.<locals>.<listcomp>b   s_     & & &!QV_& & GGAtAw}%% & & & &r{      )nparrayfloatr   r   rq   r~   r   rangelenreshapen_alphaabsr}   r   
zeros_likeprintmaxall)r}   valsdeltarv   rw   rx   s      @@@ry   test_tabledistr      s3   H88899$$B$?E8 D D DEJL LD h666666666666666666666666666666666666666666666666666666666666222000000-2 3 3, 4511ddd7-<, ?D-DG0 
5$	(	(Br**D111b))5111r**E2222774$$++++++++++288D"%%----------2774$$,,,,,,,,,,288D"%%----------2775"%%----------288E2&&....$...... ONNNNNs4yy)9)9NNND8D>>!!"bj11DF4LL28#EE2=//000 MLLLLs4yy)9)9LLLD8D>>!!"bj11DF4'>""EE2=//000 
"&	 	 	 	 	uSYY'7'7 	 	 	 GB
++g6 7 7 8 8 9 9 9* * * * *%D		*:*: * * *D8D>>!!"bj1n55D7111crc6?"'')))))! ! ! ! !%D		*:*: ! ! !D8D>>!!"bj1n55D7111abb5>!&&(((((& & & & & &uSYY/?/? & & &D8D>>!!"bj1n55D28CRC= %%'''''     uSYY/?/?   D8D>>!!"bj1n55D28ABB<$$&&&&&& & & & & &uQD		/B/B & & &D8D>>!!"bj1n55D28CRC= %%'''''''r{   )	__doc__numpyr   numpy.testingr   r   statsmodels.stats.tabledistr   r   ro   r{   ry   <module>r      se         > > > > > > > > 1 1 1 1 1 1W( W( W( W( W(r{   