U
    Kvf                     @   sn   d Z ddlZddlmZmZ ddlmZmZm	Z	 dd Z
dd Zed	krjddlZeed
ddg ed dS )z5
Created on Sun Apr 17 22:13:36 2011

@author: josef
    N)assert_assert_almost_equal)skewnorm	skewnorm2ACSkewT_genc            	      C   sJ  t dddddg} tdddd	d
gd}tt j|| ddd tdddd	d
gd}tt j|| ddd t dddddg}tdddd	d
gd}t t 	|| }t t 	|| |d  }d||||f }t
||dd tdddd	d
gd}t t 	|| }t t 	|| |d  }d||f }t
||d|d d S )Ng7R:X-煛E;Q63E?73?:6?r         
   vIh%<=ZrtolZatol        '5S>?A?zbC?gJz5z$maxabs=%15.13g, maxrel=%15.13g
%r
%r)decimalzmaxabs=%15.13g, maxrel=%15.13g)r   err_msg)nparrayr   pdfr   allcloser   cdfmaxabsr   )	pdf_rZpdf_snZpdf_sn2cdf_rZcdf_snZmaxabsZmaxrelmsgZcdf_sn2 r#   W/tmp/pip-unpacked-wheel-2v6byqio/statsmodels/sandbox/distributions/tests/test_extras.pytest_skewnorm   s4      r%   c                  C   s  t  } ddddddg}tddd	d
ddg}| |dd}tj||ddd tj||ddd tddddddg}| |dd}tj||ddd tj||ddd tddddddg}| |dd}ttj||ddd td d!d"d#d$d%g}| |dd}ttj||d&dd td'd(d)d*d+d,g}| |dd}ttj||ddd td-d.d/d0d1d2g}| |dd}ttj||ddd d S )3Nr   r   g      r   r   r   r   r   g,>r	   r
   r   i@B r   gư>r   g?r   g   0>r   r   r   gq7 %U>gO>g<vB?gXyK?gሃVB?g#\?   r   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   )Zskewtxr    Zpdf_str!   Zcdf_str#   r#   r$   
test_skewt3   sh                r(   __main__z-vvsz-xz--pdbDone)__doc__Znumpyr   Znumpy.testingr   r   Z(statsmodels.sandbox.distributions.extrasr   r   r   r%   r(   __name__Zpytestmain__file__printr#   r#   r#   r$   <module>   s   $7