M/PhdZddlZddlmZddlmZmZmZm Z m Z m Z m Z ddl mZmZmZdZdZdZd Zed kr*eeeedSdS) z> Created on Wed Feb 19 12:39:49 2014 Author: Josef Perktold N)stats) SkewNorm_genskewnorm ACSkewT_gen NormExpan_gen pdf_moments ExpTransf_gen LogTransf_gen)mc2mvskmnc2mcmvsk2mncc :ttjddtjdtjddtjdz ttjddtjddtjddtjddz ttjddtjddtjddtjddz tjdd}td||td tjdtjd d}td||td tjd ttj ddtj dtj ddtj dz ttj ddtj ddtj ddtj ddz ttj d dtj d dtj d dtj d dz dS) Nrisizesample mean var: theoretical mean varg?) printrpdfrnormchirvsmeanvarcdfsf)rs t/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/statsmodels/sandbox/distributions/examples/ex_extras.py example_nr"s (,q  UZ^^A.. Qq0A0AEJNNSTDUDU0UVVV (,q   a 2 2HL44H4H59==YZ[\K]K]4]^^^ (,r% %)--!"4"4hl2e6L6Luy}}]^_`OaOa6abbb ,qc " " "C sxxzz37799555 (."3"3444 ,qc " " "C sxxzz37799555 (."3"3444 (,q  UZ^^A.. Qq0A0AEJNNSTDUDU0UVVV (,q   a 2 2HL44H4H59==YZ[\K]K]4]^^^ (+d4 %),,tA"6"6 D8N8NQVQZQ]Q]^bcdQeQe8efffffcrt}|ddd}td||td|ddtdt jdt|ddtjt jdd}t||dS) N rrrrrz t mean varr) rrrrrrtnpabs)skewtrs r! example_Tr*"s MME ))Bqc) " "C sxxzz37799555 %++b"3"3444 , b))*** %++b   &RT** + +C #((**cggii     r#cPt}|jdd}t|d}tj|dd}t dt ||d }t d t |t d t t j|t j|z t d t t|j t |j t|}t|}t d t |t dt |t|}t dt d|ddz |ddzt ||ddz |ddzgt ||ddz |ddzgdS)Nrdrsample)modezsample: mu,sig,sk,kurmvsk)momentsznormexpan: mu,sig,sk,kurzmvsk diff distribution - sampleznormexpan attributes mvskzcentral momentsznon-central momentsz pdf approximation from momentszpdf atrr)rrrrdescriberr'arrayr cntr0r r rr)rr normexpansmvskdmvskmncmcpdffns r!examples_normexpandr;.s~~H (,qc " " "Cc111I N3   #E !""" %LLL OOFO + +E $%%% %LLL +,,, "(5//BHUOO +,,, %&&& ')- !!! ). 5//C B  "III     #JJJ OOE ,--- (BqE!GBqE!G$$$ %%AqAq! " "### )--AqAq) * *+++++r#c<tdttjdd}t|dttjddt|ttjdt|dtdttj}t| dd ttj dd t| d d ttj d d td ttj }t| d ttj d dttj }t| d tj d d | d dddS)NzResults for lognormalrzLog transformed normal general)anamerrrzResults for expgammar%r/zResults for loglaplacegUUUUUU?) rr rrrlognormrr gamma_cdfloggammalaplace loglaplace) lognormalg loggammaexpg loglaplacegloglaplaceexpgs r!examples_transfrKOs !"""uzQ7WXXXJ *..   %-  Aa !!! *     %-  a !!! *..a. !!! !!! --L ,  Ab ! !""" %.  Qr " "### ,  Ab ! !""" %.  Qr " "###  "### ..K +  1   %   q # #$$$"5=11N .  a !!! 1!D!!!!!r#__main__)__doc__numpyr'scipyr(statsmodels.sandbox.distributions.extrasrrrrrr r statsmodels.stats.moment_helpersr r r r"r*r;rK__name__r#r!rTs9>>>>>>>>>>>>>>>>>>GFFFFFFFFF g g g !!!,,,B("("("V z IKKK IKKKO r#