
    IR-e.                         d dl Zd dlZd dlmZmZ d dlmZmZ d dl	m
Z
 d Zej                            e
 d          d             Zej                            e
 d          d	             Zd
 ZdS )    N)assert_allcloseassert_equal)jackknife_resamplingjackknife_stats)	HAS_SCIPYc                      t          j        g d          } t          j        g dg dg dg dg          }t          |t          |                      d S )N)            )r
   r   r   )r	   r   r   )r	   r
   r   )r	   r
   r   )nparrayr   r   dataanswers     Blib/python3.11/site-packages/astropy/stats/tests/test_jackknife.pytest_jackknife_resamplingr      s\    8LLL!!DXyyy)))YYY			BCCF-d3344444    zrequires scipy)reasonc                      t          j        d          } d}t          |t          | t           j                  dd         d           d S )N)	s            i  i  i     i  )gz,C'p@g        g^2!I@r   r   -C6?atol)r   r   r   r   meanr   s     r   test_jackknife_statsr       sI     8ABBD&FFOD"'::1Q3?dKKKKKKr   c                  F   t          j        g d          } t          j        | d          } | d d df         } d}d }t          |t	          | |d          dd         d	
           t          j        d          }t          |t	          | |d          d         d	
           d S )N)40   *   $   !            '   r#   &   r#   r$   r&      r#   r%      r&   )   +   -   "      r,         r   r+   r%   r0      r(   r0   r-   r          r*         /      r-   r-   r7   r4      r1      r4   r-   (   )r
   r	   )g_Qr\@gnlg9D6@c                     t          j        | t          j        |           z
  | t          j        |           z
  z            t          |           z  S )N)r   sumr   len)xs    r   mle_varz3test_jackknife_stats_conf_interval.<locals>.mle_var2   s:    vq271::~!bgajj.9::SVVCCr   gffffff?r   r   r   r   )gޓZQ@gKc@)r   r   reshaper   r   )r   r   rC   s      r   "test_jackknife_stats_conf_intervalrE      s    
 8	
 	
 	
 D :dG$$D1:D -FD D D FOD'4@@1EDQQQQ X+,,FFOD'4@@C$OOOOOOr   c                      t          j        t                    5  t          t	          j        d          t          j        d           d d d            d S # 1 swxY w Y   d S )Nr
   r#   )confidence_level)pytestraises
ValueErrorr   r   aranger    r   r   test_jackknife_stats_exceptionsrM   <   s    	z	"	" D D	!bgCCCCD D D D D D D D D D D D D D D D D Ds   /AAA)numpyr   rH   numpy.testingr   r   astropy.stats.jackknifer   r   "astropy.utils.compat.optional_depsr   r   markskipifr    rE   rM   rL   r   r   <module>rT      s         7 7 7 7 7 7 7 7 I I I I I I I I 8 8 8 8 8 85 5 5 	M*:;;L L <;L 	M*:;;P P <;P>D D D D Dr   