!gdZddlZddlmZddlmZmZmZm Z m Z m Z m Z ddl mZmZmZdZdZdZd Zed k(reeeeyy) 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 ttjddtjjdtjddtjjdz ttjddtj jddtjddtj jddz ttjddtj jddtjddtj jddz tj dd}td|j|jtd tjdtj d d}td|j|jtd tjd ttjddtjjdtjddtjjdz ttjddtj jddtjddtj jddz ttjd dtj jd dtjd dtj jd dz y) Nrisizesample mean var: theoretical mean varg?) printrpdfrnormchirvsmeanvarcdfsf)rs q/var/www/dash_apps/app1/venv/lib/python3.12/site-packages/statsmodels/sandbox/distributions/examples/ex_extras.py example_nr"s! (,,q UZZ^^A. Qq0AEJJNNSTDU0UV (,,q  a 2HLL44H599==YZ[\K]4]^ (,,r% %))--!"4hll2e6Luyy}}]^_`Oa6ab ,,qc "C sxxz37795 (.."34 ,,qc "C sxxz37795 (.."34 (,,q UZZ^^A. Qq0AEJJNNSTDU0UV (,,q  a 2HLL44H599==YZ[\K]4]^ (++d4 %)),,tA"6 D8NQVQZQZQ]Q]^bcdQe8efc t}|jddd}td|j|j td|j ddtdt j j dt|j ddtjt j jdd}t|j|j y) N rrrrrz t mean varr) rrrrrrtnpabs)skewtrs r! example_Tr*"s ME ))Bqc) "C sxxz37795 %++b"34 , b)* %++b  &&RT* +C #((*cggi r#cpt}|jdd}t|d}tj|dd}t dt ||jd }t d t |t d t t j|t j|z t d t t|jt |jt|}t|}t d t |t dt |t|}t dt d|ddz |ddzt ||ddz |ddzgt |j|ddz |ddzgy)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;.sb~H (,,qc "Cc1I NN3  #E !" %L OOFO +E $% %L +, "((5/BHHUO +, %& ')-- ! ).. 5/C B  "I   #J OE ,- (BqE!GBqE!G$ %AqAq! "# )--AqAq) *+r#ctdttjdd}t|j dttj j ddt|jttj jdt|j dtdttj}t|jdd ttjj dd t|jd d ttjj d d td ttj}t|jd ttjj d dttj}t|jd tjj d d |jd ddy)NzResults for lognormalrzLog transformed normal general)anamerrrzResults for expgammar%r/zResults for loglaplacegUUUUUU?) rr rrrlognormrr gamma_cdfloggammalaplace loglaplace) lognormalg loggammaexpg loglaplacegloglaplaceexpgs r!examples_transfrKOs !"uzzQ7WXJ *..  %--  Aa ! *    %--  a ! *..a. ! ! -L ,  Ab !" %..  Qr "# ,  Ab !" %..  Qr "#  "# .K +  1  %    q #$"5==1N .  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!rTsf>>>GF g !,B("V z K K r#