
[c        
   @` s  d  d l  m Z m Z m Z m Z d  d l Z d  d l Z d  d l m	 Z	 m
 Z
 y d  d l Z Wn e k
 ry e Z n Xe Z y d  d l Z Wn e k
 r e Z n Xe Z d d l m Z d d l m Z d d l m Z d d	 l m Z d d
 l m Z d d l m Z d   Z d   Z e j  j! d  d    Z" d   Z# d   Z$ e j  j! d  d    Z% e j  j! d  d    Z& d   Z' d   Z( e j  j! d  d    Z) d   Z* e j  j+ d  d    Z, d   Z- d   Z. d   Z/ d   Z0 d   Z1 d    Z2 d!   Z3 d"   Z4 d#   Z5 e j  j! d  e j  j6 d$ d% d& d' d( g  d)     Z7 d*   Z8 e j  j! d  d+    Z9 d,   Z: d-   Z; e j  j! d  d.    Z< e j  j! d  d/    Z= e j  j! d  d0    Z> e j  j! d  d1    Z? e j  j! d  e j  j6 d2 d  d3 d d d4 d5 d6 g  d7     Z@ e j  j! d  d8    ZA e j  j! d9  d:    ZB e j  j! d  d;    ZC e j  j! d  d<    ZD e j  j! d=  d>    ZE d S(?   i    (   t   absolute_importt   divisiont   print_functiont   unicode_literalsN(   t   assert_equalt   assert_allclosei   (   t   funcsi   (   t   units(   t   range(   t   catch_warnings(   t   NUMPY_LT_1_10(   t   NumpyRNGContextc          C` s]  t  d  Kt j j d  }  t j |   } |  j d<  }  t j |  d d } t |  d k sj t  | j	 |  j	 k  s t  t j |  d d } t |  d k s t  | j	 |  j	 k  s t  t j
 d>  } t j g  t t |   D] } t | | d   ^ q  j d d	 d
  } t j |  } | d k sBt  t j | d d } t | d d d d d g d d d d d g d d d d d g d d d d d g g  t j | d d } t | d  d! d" d# d$ g d% d& d' d( d) g d* d+ d, d- d. g g  t j | d d/ } t | d0 d1 d2 d3 g d4 d5 d6 d7 g d8 d9 d: d; g g  Wd  QXd  S(?   Ni90  i'  i
   i  t   axisi   i    i   i   i   g     Xx@g     @j@g     l@g     @o@g     p@g      r@g     `s@g     t@g     u@g      w@g     `x@g     y@g     z@g      |@g     `}@g     ~@g     @g     @g     0@g     Ё@g     p@g     ;@g     @@@g     B@g     @E@g     G@g     _@g     `@g     0a@g     a@g     pb@g     pl@g     m@g     m@g     Pn@g     n@i   g      @g       @g      *@g      2@g      7@g      <@g     @@g      C@g     E@g      H@g     J@g      M@(   i
   i  i   i<   (   R   t   npt   randomt   randnR   t   median_absolute_deviationt   reshapet   lent   AssertionErrort   sizet   aranget   arrayR   t   sumR   (   t   randvart   madt   xt   it   a(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_median_absolute_deviation   s6    Kc          C` s  t  j j d d d g  }  t j |   d k s6 t  t  j j d d d g d d d d g }  t j |   d k s{ t  t j |   t j |  |  j  k s t  t  j j d d g d d g g d d d g d d g g }  t j |   t j |   d k st  t j |  j  d k s)t  t  j j	 t j |  d d j d d g  t  j j	 t j |  d d j d d g  d  S(	   Ni   i   i   i   t   maski    g      ?R   (
   R   t   maR   R   R   R   R   t   datat   testingt   assert_array_equal(   R   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt%   test_median_absolute_deviation_maskedA   s    * <	"	u   NUMPY_LT_1_10c          C` s   t  j d d d t  j g d d t  j d g g  }  t t j |  d t  j d d d d g  t  j j |   }  t j |   d k s t	  d  S(   Ni   i   i   i   i   t   funcR   (
   R   R   t   nanR   R   R   t	   nanmedianR   t   masked_invalidR   (   R   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt#   test_median_absolute_deviation_nans`   s    c          C` s   t  j d  t  j d  d  d   t  j t  j f }  t t j |  d d t  j d   t t j |  d t  j d d g  t  j d   d  S(	   Ni   i   i   R   i   i   (   i   i   i   (   i   i   (	   R   t   onesR   t   newaxisR   R   R   t   zerosR   (   R   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt,   test_median_absolute_deviation_multidim_axisk   s
    5	c          C` sk   t  j d d d g  t j }  t j |   } | j |  j k sF t  | j t j |  j  k sg t  d  S(   Ni   i   i   (	   R   R   t   ut   mR   R   t   unitR   t   value(   R   R   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt'   test_median_absolute_deviation_quantitys   s    u   not HAS_SCIPYc       
   C` s  d }  d d d g } x d d d g D] } t  j | |  d | d d	 } | d k | d k @j   sk t  t  j | |  d | d d
 } | d k | d k @j   s% t  q% Wd d d d d g } d }  d } t  j | |  d | d d
 } t j d d d d d g d d d d d g g  } t | | d d d d t j g  | D]$ } t  j | |  d | d d
 ^ qG j   } t | | d d d d t  j | |  d | d d } t j d d d d d  g d! d" d# d$ d% g g  } t | | d d d d t j g  | D]$ } t  j | |  d | d d ^ q
 j   } t | | d d d d t  j d d d d& } t | d  t  j d d d d& } t | d  t  j d' d( d d) d d& } t | d d d t j d*   t | d d d t j d*   d } d }  xG d& d	 d
 d g D]3 } t  j | |  d | } | j	 d+ k s!t  q!Wt j |  } xG d& d	 d
 d g D]3 } t  j | |  d | } | j	 d, k szt  qzWt j |   }  xG d& d	 d
 d g D]3 } t  j | |  d | } | j	 d- k st  qWt j d d d g  } xG d& d	 d
 d g D]3 } t  j | |  d | } | j	 d. k s5t  q5Wt j d d d g  }  xG d& d	 d
 d g D]3 } t  j | |  d | } | j	 d/ k st  qWd  S(0   Ni   i    i   g        g      ?g      ?t   conft   intervalu   wilsonu   jeffreysi   i   i   i   gffffff?gMb?gp=
ף?gn?gʡE?g㥛 ?gx&1?g#~j?gPn?gZd;O?t   atolgMbP?t   rtolu   flatgTN?g6Yѵ?gW?gW?g.4i?g+?gO@a?gOc*?ge1?u   waldi  i  g鲘?g     @@(   i   (   i   (   i   (   i   i   (   i   i   (
   R   t   binom_conf_intervalt   allR   R   R   R   t	   transposet   sqrtt   shape(   t   nt   kR2   t   rest   resultt   tablet   kvalR3   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_binom_conf_interval   sj    "&	7	7	""c          C` s  d }  t  j d d d  } t  j t |  d t  j } t j | | d |  \ } } } } | j |  f k su t  | j |  f k s t  | j |  f k s t  | j d |  f k s t  | d k j	   s t  t
 | (t j | | d |  \ } } } } | d k j	   st  d  S(	   Ni   g        g      $@id   t   dtypet   binsi   g      ?(   R   t   linspaceR)   R   t   boolR   t   binned_binom_proportionR:   R   R7   t   False(   t   nbinsR   t   successt   bin_ctrt   bin_hwt   pt   perr(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_binned_binom_proportion   s    c          C` s'  t  j d d d d d d  }  d |  k s0 t  t  j d d d d d d d  }  d |  k sc t  t  j d d d d d d  }  |  d k  s t  t  j d d d d d d  }  |  d k  s t  t  j d d d d d d  }  |  d k  s t  t  j d d d d d d  }  |  d k s#t  d  S(   Ni   i   i    g      @i   i   (   R   t   signal_to_noise_oir_ccdR   (   R>   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_signal_to_noise_oir_ccd   s    !c          C` s  t  j d d d d d d d d d	 d
 g
  }  t  j d d d d d d
 d d d d g
 d d d d d d d d d d g
 g  } t d   t | t j |  d   Wd  QXt d  E t  j t j |  d d t  j  } t t  j |   | d d Wd  QXd  S(   Ni   i   i   i   i   i   i   i   i	   i    i*   i'  t   bootfuncR4   g{Gz?(   R   R   R   R   R   t	   bootstrapt   meanR   (   t   bootarrt   answert
   bootresult(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_bootstrap   s    -'*$c          ` s  d d l  m   t d   t j d d d d d d	 d
 d d d g
 d d d d d	 d d d d	 d g
 g  j }  t j d  }   f d   } t j |  d d | } t | | d d Wd  QXt d   t j d d d d d d	 d
 d d d g
 d d d d d	 d d d d	 d g
 g  j }  t j d  }   f d   } t j |  d d | } t | | d d Wd  QXt d  l t j d d d f  }   f d   } t j |  d d | } | j	 d  k st
  t | | d d Wd  QXd  S(!   Ni    (   t	   spearmanri*   i   i   i   i   i   i   i   i   i	   gw/?gKq?c         ` s     |   d S(   Ni    (    (   R   (   RX   (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyRQ     s    RQ   R4   gMbP?g鷯?g#?c         ` s     |   d S(   Ni   (    (   R   (   RX   (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyRQ   *  s    gq?gAf?g2ı.n?gʡE?c         ` s
     |   S(   N(    (   R   (   RX   (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyRQ   7  s    (   gw/?gKq?(   g鷯?g#?(   gq?g鷯?(   gAf?g#?(   g2ı.n?gʡE?(   i   i   (   t   scipy.statsRX   R   R   R   t   TR   RR   R   R:   R   (   RT   RU   RQ   RV   (    (   RX   s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_bootstrap_multiple_outputs  s6    '-	'-				c       
   C` sN   t  d  < t j j d d d d	 }  t t j |   d d d Wd  QXd  S(
   Ni90  i   i   R   id   g       @R5   g?(   id   id   (   R   R   R   t   normalR   R   t   mad_std(   R    (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_mad_std@  s    u   not NUMPY_LT_1_10c       
   C` s   t  d   t j j d d d d }  t j j t j |   |   }  t j |   } t j	 |  sj t
  t j |  d <t j |  d t } t j	 |  s t
  t   ? t j |   } t j	 |  s t
  t j |  s t
  Wd  QXWd  QXd  S(	   Ni90  i   i   R   i
   t
   ignore_nan(   i
   i
   (   i   i   (   R   R   R   R\   R   t   masked_wheret   isnanR   R]   t   isscalarR   R%   t   TrueR	   (   R    t   rslt(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_mad_std_scalar_returnF  s    
c          C` s   t  d   t j j d d d d	 }  t j |  d
 <t   b } t j |  d t } t	 r | d } t
 | j  j d  s t  n t j |  s t  Wd  QXWd  QXd  S(   Ni90  i   i   R   i
   R_   i    u    Numpy versions <1.10 will return(   i
   i
   (   i   i   (   R   R   R   R\   R%   R	   R   R]   RG   R
   t   strt   messaget
   startswithR   Ra   (   R    t   warnsRd   t   w(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_mad_std_warnsY  s    
!c       
   C` s   t  d  w t j d d g  }  t j |  (t j j d d d d |  d d  d d  f <t t j |  d	 t	 d
 d d Wd  QXt
 s t j t j d d d d d t j g   s t  n  t t j d d d d d t j g d	 t	 d  d  S(   Ni90  if   i   i   R   id   i   iR_   g       @R5   g?i   i   gtV?(   id   id   (   R   R   t   emptyR%   R   R\   R   R   R]   Rc   R
   Ra   R   (   R    (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_mad_std_withnang  s    
1(6*c          C` s   t  j d d d d g d d d d g g  }  t  j d d d d g  } t  j d d g  } t t j |  d d	 |  t t j |  d d |  d  S(
   Ni   i   i   i   g#>ԍ@g.?g.?R   i    (   R   R   R   R   R]   (   R    t   result_axis0t   result_axis1(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_mad_std_with_axist  s    c          C` s   t  j d d d d t  j g d d d d t  j g g  }  t  j d d d d t  j g  } t  j d d g  } t t j |  d d	 d
 t |  t t j |  d d d
 t |  d  S(   Ni   i   i   i   g#>ԍ@g.?g.?R   i    R_   (   R   R   R%   R   R   R]   Rc   (   R    Rn   Ro   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_mad_std_with_axis_and_nan  s    "c          C` s   t  j d d d d t  j g d d d d t  j g g  }  t j |  d d d t } t  j j |  sm t  t  j j	 t  j
 |   |   }  t j |  d d d t } t  j j |  s t  d  S(   Ni   i   i   i   R   i    R_   (   R   R   R%   R   R]   Rc   R   t   isMaskedArrayR   R`   Ra   (   R    R>   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt)   test_mad_std_with_axis_and_nan_array_type  s    c          C` s>   d t  j d t  j d   }  t t j |  d d d d  S(   Ng       @g      ?R5   gư>(   R   R9   t   logR   R   t   gaussian_fwhm_to_sigma(   t   fwhm(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_gaussian_fwhm_to_sigma  s     c          C` sB   d d t  j d t  j d   }  t t j |  d d d d  S(   Ng      ?g       @R5   gư>(   R   R9   Rt   R   R   t   gaussian_sigma_to_fwhm(   t   sigma(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_gaussian_sigma_to_fwhm  s    $c           C` s   t  t j t j d  d  S(   Ng      ?(   R   R   Ru   Rx   (    (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt$   test_gaussian_sigma_to_fwhm_to_sigma  s    	c           C` s    t  t j d d d d  d  S(   Ni   R3   u   root-ni   i   (   i   i   (   R   R   t   poisson_conf_interval(    (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt    test_poisson_conf_interval_rootn  s    u   intervalu   root-n-0u   pearsonu   sherpagehrelsu   frequentist-confidencec         C` s;   d } t  t j | d d t j | d |  d d d  S(   Nid   R3   u   root-nR5   g{Gz?(   R   R   R|   (   R3   R;   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_poisson_conf_large  s    c          C` s   t  j d  }  t t j |  d d t j |  d d d d  d   d  d  d  f t  j |    t  j t  j t j |  d d   s t	  d  S(	   Ni   i   i   R3   u   root-n-0i    (   i   i   i   (   i    i    i    (
   R   R+   R   R   R|   t   Nonet	   ones_liket   anyRa   R   (   R;   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt#   test_poisson_conf_array_rootn0_zero  s
    =c          C` s   t  j d  }  t t j |  d d t j |  d	 d d d  d   d  d  d  f t  j |    t  j t  j t j |  d d   s t	  d  S(
   Ni   i   i   R3   u   frequentist-confidencei    u   root-n-0(   i   i   i   (   i    i    i    (
   R   R+   R   R   R|   R   R   R   Ra   R   (   R;   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt3   test_poisson_conf_array_frequentist_confidence_zero  s    =c          C` su   d d d g }  t  t j |  d d d d d g d d d g g  t j t j t j |  d d   sq t  d  S(   Ni    R3   u   root-n-0i   (   R   R   R|   R   R   Ra   R   (   R;   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt"   test_poisson_conf_list_rootn0_zero  s
    c          C` s   d t  j d
  }  t t j |  d d t j |  d d d d  d   d  d  d  f t  j |    d |  d <t  j t  j t j |  d d   s t	  d  S(   Ni   i   i   i   R3   u   root-n-0i    i   i   (   i   i   i   (   i    i    i    (   i   i   i   (
   R   R)   R   R   R|   R   R   R   Ra   R   (   R;   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_poisson_conf_array_rootn0  s    =
c          C` s   d t  j d
  }  t t j |  d d t j |  d d d d  d   d  d  d  f t  j |    d |  d <t  j t  j t j |  d d   s t	  d  S(   Ni   i   i   i   R3   u   frequentist-confidencei    i   i   (   i   i   i   (   i    i    i    (   i   i   i   (
   R   R)   R   R   R|   R   R   R   Ra   R   (   R;   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_poisson_conf_array_fc  s    =
c          C` s   t  j d' d( d) d* d+ d, d- d. d/ d0 d1 g  }  t t j |  d! d!  d f d" d# |  d! d!  d d!  f j d$ d% d& d% d! S(2   u6   Test intervals against those published in Gehrels 1986i    g~jt?i   g/$?gffffff
@i   g-?gZd;O@i   gZd;?gy&1@i   g㥛  @gx@i   gQ@g$ @i   g(\@g^I+#@i   gV@g
ףp=%@i   g!rh@gfffff'@i	   g|?5^:@gQ8*@i
   gMb@g
ףp=,@NR3   u   frequentist-confidenceR5   gMbP?R4   (   i    i    g~jt?(   i   g/$?gffffff
@(   i   g-?gZd;O@(   i   gZd;?gy&1@(   i   g㥛  @gx@(   i   gQ@g$ @(   i   g(\@g^I+#@(   i   gV@g
ףp=%@(   i   g!rh@gfffff'@(   i	   g|?5^:@gQ8*@(   i
   gMb@g
ףp=,@(   R   R   R   R   R|   RZ   (   t   nlh(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt0   test_poisson_conf_frequentist_confidence_gehrels  s    	c          C` s   t  j d' d( d) d* d+ d, d- d. d/ d0 d1 g  }  t t j |  d! d!  d f d" d d# d$ j |  d! d!  d d!  f d% d& d! S(2   u   Test intervals against those published in Gehrels 1986

    Note: I think there's a typo (transposition of digits) in Gehrels 1986,
    specifically for the two-sigma lower limit for 3 events; they claim
    0.569 but this function returns 0.59623...

    i    i   g$C@i   gZd;O?gZd@gq=
ףp?g1Zd@i   gOn?g!@i   g!rh?gHz$@i   gS?gp=
ף'@i   gMbX9@gq=
ףp*@i   gMb@gGz.-@i   gT㥛 @gzG/@i	   gv/@gQE1@i
   gK7A@gGz2@NRy   R3   u   frequentist-confidenceR5   g{Gz?(   i    i   i    g$C@(   i   i   gZd;O?gZd@(   i   i   gq=
ףp?g1Zd@(   i   i   gOn?g!@(   i   i   g!rh?gHz$@(   i   i   gS?gp=
ף'@(   i   i   gMbX9@gq=
ףp*@(   i   i   gMb@gGz.-@(   i   i   gT㥛 @gzG/@(   i	   i   gv/@gQE1@(   i
   i   gK7A@gGz2@(   R   R   R   R   R|   RZ   (   R   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt7   test_poisson_conf_frequentist_confidence_gehrels_2sigma   s    		"	c          C` s   t  j d( d) d* d+ d, d- d. d/ d0 d1 d2 g  }  t t j |  d! d!  d f d" d d# d$ j |  d! d!  d d!  f d% d& d' t d! S(3   u6   Test intervals against those published in Gehrels 1986i    i   g;On@i   g_vOV?g!@i   ghs?g=
ףp%@gA`"?g\(\)@i   g(\?gHz,@i   gMbX?gHz0@i   g?gQ1@i   gS?gQ+3@i   g%C @gq=
ף4@i	   g/$@gGz.6@i
   gy&1@gp=
ף7@NRy   R3   u   frequentist-confidenceR5   g{Gz?t   verbose(   i    i   i    g;On@(   i   i   g_vOV?g!@(   i   i   ghs?g=
ףp%@(   i   i   gA`"?g\(\)@(   i   i   g(\?gHz,@(   i   i   gMbX?gHz0@(   i   i   g?gQ1@(   i   i   gS?gQ+3@(   i   i   g%C @gq=
ף4@(   i	   i   g/$@gGz.6@(   i
   i   gy&1@gp=
ף7@(   R   R   R   R   R|   RZ   Rc   (   R   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt7   test_poisson_conf_frequentist_confidence_gehrels_3sigma  s    	"	u   ni   i
   i   id   c         C` s=   t  t j |  d d d t j |  d d d d d d  S(   NR3   u   sherpagehrelsi   u   frequentist-confidenceR5   g{Gz?(   R   R   R|   (   R;   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_poisson_conf_gehrels86/  s    c          C` s   t  t j d d d  d d d t j d d g d	 d
 d d g d d d g }  t  |  d d  d f d d d t  |  d d  d f d d d d S(   u   Test that the lower-level routine gives the snae number.

    Test numbers are from table1 1, 3 in
    Kraft, Burrows and Nousek in
    `ApJ 374, 344 (1991) <http://adsabs.harvard.edu/abs/1991ApJ...374..344K>`_
    g      @g      @gGz?i    gףp=
W%@R5   gMbP?g      @u   kraft-burrows-nousekt
   backgroundg       @t	   conflevelg?Ni   gQ?g{G!@g{Gzt?(   i    gףp=
W%@(   i    gףp=
W%@(   gQ?g{G!@(   R   R   t   _scipy_kraft_burrows_nousekR|   (   R2   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_scipy_poisson_limit8  s    #u   not HAS_MPMATHc           C` sH   t  t j d d d  d d d t  t j d d	 d
  d d d d  S(   Ng      @g       @g?gQ?g{G!@R5   g{Gzt?g      @g      @gGz?i    gףp=
W%@gMbP?(   gQ?g{G!@(   i    gףp=
W%@(   R   R   t   _mpmath_kraft_burrows_nousek(    (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_mpmath_poisson_limitI  s    c          C` s@  t  j t  # }  t j d d g d d d Wd  QXd t |  j  k sO t  t  j t  ) }  t j d d g d d d	 d
 g Wd  QXd t |  j  k s t  t  j t  ) }  t j d d g d d d	 d
 g Wd  QXd t |  j  k s t  t  j t   }  t j d d  Wd  QXd t |  j  k s<t  d  S(   Ni   i   u   root-nRy   i   u   Only sigma=1 supportedu   pearsonR   g      @g       @u   background not supportedu   sherpagehrelsR   u   conflevel not supportedi   u   foou   Invalid method(   t   pytestt   raisest
   ValueErrorR   R|   Rf   R0   R   (   t   e(    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   test_poisson_conf_value_errorsQ  s    "(c          C` s   t  j t  # }  t j d d d d d d Wd  QXd t |  j  k sO t  t  j t   }  t j d d d d Wd  QXd t |  j  k s t  t  j t  # }  t j d d d d	 d d
 Wd  QXd t |  j  k s t  d  S(   Ng      @u   kraft-burrows-nousekR   g      @R   ic   u   number between 0 and 1u   Set conflevel for methodg      gGz?u   Background must be(   R   R   R   R   R|   Rf   R0   R   (   R   (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt"   test_poisson_conf_kbn_value_errorse  s    u   HAS_SCIPY or HAS_MPMATHc           C` s9   t  j t  $ t j d d d d d d d Wd  QXd  S(   Ng      4@R3   u   kraft-burrows-nousekR   g      $@R   gffffff?(   R   R   t   ImportErrorR   R|   (    (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt!   test_poisson_limit_nodependenciesy  s    (F   t
   __future__R    R   R   R   R   t   numpyR   t   numpy.testingR   R   t   scipyR   RG   t	   HAS_SCIPYRc   t   mpmatht
   HAS_MPMATHt    R   R   R-   t   extern.six.movesR   t   tests.helperR	   t   utils.compatR
   t
   utils.miscR   R   R#   t   markt   skipifR(   R,   R1   RA   RN   RP   RW   R[   R^   t   xfailRe   Rk   Rm   Rp   Rq   Rs   Rw   Rz   R{   R}   t   parametrizeR~   R   R   R   R   R   R   R   R   R   R   R   R   R   R   (    (    (    s=   lib/python2.7/site-packages/astropy/stats/tests/test_funcs.pyt   <module>   sx   "

	"			O		3															
6