
&]\c           @` ss   d  d l  m Z m Z m Z d  d l Z d  d l m Z m Z d  d l	 m
 Z
 m Z d   Z d   Z d   Z d S(   i    (   t   divisiont   print_functiont   absolute_importN(   t   assert_allcloset   assert_equal(   t   tukeylambda_variancet   tukeylambda_kurtosisc          C` s6  t  d  }  t |  t j d d d d t d  } t | d d d t  d  }  t |  d	 t j d d t d  } d t j d t j d	 d d d } t | | d d t  d  }  t |  d d d t d  } t | d d d t  d  }  t |  d d d t d  } t | d d d d S(   s/   Compare results with some known exact formulas.i    i   i   t   atolg-q=g333333?g|=g      ?i   g      @i   g      ?g333333i   Ng?gUUUUUU?gUUUUUU?(   R   R   t   npt   piR   (   t   vart   kurtt   desired(    (    sG   lib/python2.7/site-packages/scipy/stats/tests/test_tukeylambda_stats.pyt"   test_tukeylambda_stats_known_exact
   s"    
(c          C` sW  t  d d d d  }  t  d d d d  } d d d g d	 d
 d g d d d g d d d g d d d g d d d g d d d g d d d g d d d  g d! d" d# g d$ d% d& g g } xO | D]G \ } } } t |  } t | | |  t |  } t | | |   q Wt |   \ } } } t |  } t | | |  t |  } t | | |   d' S((   sA   Compare results with some values that were computed using mpmath.R   g|=t   rtoli    g-q=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?g333333?gڽP?g@᜿g      $@gds5O?gX\@g      4@g ?gc@N(   t   dictR   R   R   t   zip(   t   a10t   a12t   datat   lamt   var_expectedt   kurt_expectedR
   R   (    (    sG   lib/python2.7/site-packages/scipy/stats/tests/test_tukeylambda_stats.pyt   test_tukeylambda_stats_mpmath-   s.    c          C` sx   d d g }  t  |   } t | t j t j t j g   d d g }  t |   } t | t j t j t j g   d S(   s;   Test values of lambda outside the domains of the functions.g      g      g      пN(   R   R   R   t   arrayt   nant   infR   (   R   R
   R   (    (    sG   lib/python2.7/site-packages/scipy/stats/tests/test_tukeylambda_stats.pyt   test_tukeylambda_stats_invalidO   s    "(   t
   __future__R    R   R   t   numpyR   t   numpy.testingR   R   t   scipy.stats._tukeylambda_statsR   R   R   R   R   (    (    (    sG   lib/python2.7/site-packages/scipy/stats/tests/test_tukeylambda_stats.pyt   <module>   s   	#	"