
    (ph                     >    S SK rS SKJrJr  S SKJrJr  S rS r	S r
g)    N)assert_allcloseassert_equal)tukeylambda_variancetukeylambda_kurtosisc                     [        S5      n [        U [        R                  S-  S-  SS9  [	        S5      n[        USSS9  [        S5      n [        U S	[        R                  -
  SS9  [	        S5      nS
[        R                  S-  -
  [        R                  S	-  S-
  S-  -  S-
  n[        XSS9  [        S5      n [        U SSS9  [	        S5      n[        USSS9  [        S5      n [        U SSS9  [	        S5      n[        USSS9  g)z/Compare results with some known exact formulas.r         -q=)atol333333?绽|=g      ?   g?   gUUUUUU?g333333gUUUUUU?N)r   r   nppir   )varkurtdesireds      [/var/www/html/venv/lib/python3.13/site-packages/scipy/stats/tests/test_tukeylambda_stats.py"test_tukeylambda_stats_known_exactr      s     q
!CCAE2"DD#E* s
#CCRUU/$DbeeAg~"%%'A+!11A5GD. q
!CCu-"DD$U+ q
!CC."DD$U+    c                  X   [        SSS9n [        SSS9n/ SQ/ SQ/ SQ/ SQ/ S	Q/ S
Q/ SQ/ SQ/ SQ/ SQ/ SQ/nU H4  u  p4n[        U5      n[        Xd40 UD6  [        U5      n[        Xu40 U D6  M6     [	        U6 u  p4n[        U5      n[        Xd40 UD6  [        U5      n[        Xu40 U D6  g)zACompare results with some values that were computed using mpmath.r   r   )r   rtolr
   )gg;@g2V#H@)gQIg9@gk+\)@)ggxuh~@g@)gMbPg٤{i
@ghYn?)gMbP?gӦNFb:
@g4\>?)g      ?gt@g<n\?)g?gU+2C@gӤB?)gQI?gGI8@gx@As?)r   gڽP?g@᜿)g      $@gds5O?gX\@)g      4@g ?gc@N)dictr   r   r   zip)a10a12datalamvar_expectedkurt_expectedr   r   s           r   test_tukeylambda_stats_mpmathr"   +   s    
E
"C
E
"C 	9;9:9<9;9;;D -1(="3'1S1#C(3s3	 -1 (+Dz$C}
s
#CC--$DD/3/r   c                  8   SS/n [        U 5      n[        U[        R                  " [        R                  [        R
                  /5      5        SS/n [        U 5      n[        U[        R                  " [        R                  [        R
                  /5      5        g)z;Test values of lambda outside the domains of the functions.g      g      g      пN)r   r   r   arraynaninfr   )r   r   r   s      r   test_tukeylambda_stats_invalidr'   M   sh    ,C
s
#Cbhh/01-C$Drxx 012r   )numpyr   numpy.testingr   r   scipy.stats._tukeylambda_statsr   r   r   r"   r'    r   r   <module>r,      s"     7B ,F0D3r   