U
    Ãmœdž  ã                   @   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 )Nçíg7ÍR:X-ç…›�EûÔ;çQ6Ô3EˆÙ?ç73ŽåøÞ?çÔ:è6´¤»?éþÿÿÿéÿÿÿÿr   é   é   é
   ç‚vIhÂ%<=©ZrtolZatolç        ç'5 ·S> ?çú¦A¢—Øå?çzbÓC‹î?g¸ÔJzî�5z$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#   úl/home/sam/Atlas/atlas_env/lib/python3.8/site-packages/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£¢<váÔB?g÷XˆÛyKØ?gáˆƒVBÜ?g#çÐ\À©À?é   r   gÙ}ÚõÐò¾:gô$s„ø®w>gÏ„õÓ"Æ>gFNÜ#;?g 5 ·S> ?göÂgK‹`ä?g°fšé¼ì?g»½×Ùß|Û=g?+v5�Õ9?gÄ+`ó°Y?gcœßS#y?gƒÈÉm0_Ô?gVÄizEÔ?g»K ý;?À?g.Gã±ÝI?g­ÞætÑY?gG”v)g¡i?g5 ·S> ?gWosºèà?gN?rÝt”æ?)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__Ú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