
p7]c           @   s   d  Z  d d l Z d d l m Z m Z d d l m Z m Z m	 Z	 d   Z
 d   Z e d k r d d l Z e j e d d	 d
 g  d GHn  d S(   s5   
Created on Sun Apr 17 22:13:36 2011

@author: josef
iN(   t   assert_t   assert_almost_equal(   t   skewnormt	   skewnorm2t   ACSkewT_genc    	      C   s  t  j d d d d d g  }  t j d d d d	 d
 g d  } t t  j | |  d d d d  t j d d d d	 d
 g d  } t t  j | |  d d d d  t  j d d d d d g  } t j d d d d	 d
 g d  } t  j t  j	 | |   } t  j t  j	 | |  | d  } d | | | | f } t
 | | d d t j d d d d	 d
 g d  } t  j t  j	 | |   } t  j t  j	 | |  | d  } d | | f } t
 | | d d d | d  S(   Ngg7R:X-gE;gQ63E?g73?g:6?iii    i   i   i
   t   rtolgvIh%<=t   atolg        g'5S>?gA?gzbC?gJz5s$   maxabs=%15.13g, maxrel=%15.13g
%r
%rt   decimals   maxabs=%15.13g, maxrel=%15.13gt   err_msg(   t   npt   arrayR   t   pdfR    t   allcloseR   t   cdft   maxt   absR   (	   t   pdf_rt   pdf_snt   pdf_sn2t   cdf_rt   cdf_snt   maxabst   maxrelt   msgt   cdf_sn2(    (    sR   lib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/test_extras.pyt   test_skewnorm   s*    !"!"!$
!$c          C   sa  t    }  d d d d d d g } t j d d d	 d
 d d g  } |  j | d d  } t j | | d d d d t j | | d d d d t j d d d d d d g  } |  j | d d  } t j | | d d d d t j | | d d d d t j d d d d d d g  } |  j | d d  } t t j | | d d d d   t j d! d" d# d$ d% d& g  } |  j | d d  } t t j | | d d' d d  t j d( d) d* d+ d, d- g  } |  j | d d  } t t j | | d d d d   t j d. d/ d0 d1 d2 d3 g  } |  j | d d  } t t j | | d d d d   d  S(4   Niig      i    i   i   gg7R:X-gE;g,>gQ63E?g73?g:6?i@B i
   R   R   gư>g?g        g   0>g'5S>?gA?gzbC?gq7 %U>gO>g<vB?gXyK?gሃVB?g#\?i   gvIh%<=g}:g$sw>gτ">gFN#;?g 5S>?ggK`?gf?g|=g?+v59?g+`Y?gcS#y?gm0_?gVizE?gK;??g.GI?gtY?gGv)gi?g5S>?gWos?gN?rt?(   R   R	   R
   R   R   R   R    (   t   skewtt   xR   t   pdf_stR   t   cdf_st(    (    sR   lib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/test_extras.pyt
   test_skewt3   sD    	"""t   __main__s   -vvss   -xs   --pdbt   Done(   t   __doc__t   numpyR	   t   numpy.testingR    R   t(   statsmodels.sandbox.distributions.extrasR   R   R   R   R   t   __name__t   pytestt   maint   __file__(    (    (    sR   lib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/test_extras.pyt   <module>   s   	$	7	