p7]c@sdZddlmZddlZddlmZddlmZddl m Z iidd 6d 6Z d Z d Z d ZdZdZdZdZiZgZgZgZgZdZdZedkregeeD]k\ZZejeddedddkrejeddreeddeddf^qZgeeD]k\ZZejeddedddkrsejeddrseeddeddf^qsZ geeD]k\ZZejeddedddkrejeddreeddeddf^qZ!geeD]k\ZZejeddedddkroejeddroeeddeddf^qoZ"ddl#m$Z$e%edkr$dGHe$eddddgGHnd GHe$e ddddgGHd!GHe$e!ddddgGHd"GHe$e"ddddgGHndS(#sscript to test expect and moments in distributions.stats method not written as a test, prints results, renamed to prevent test runner from running it i(t iteritemsN(tstats(t expect_v2i(tdistcontitubtncfcCsR|\}}}}tj||d}tj||dd}||||fS(s>convert central moments to mean, variance, skew, kurtosis g?g@g@(tnptdivide(targstmctmc2tmc3tmc4tskewtkurt((sTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pytmc2mvsksc Cs|\}}}}|}|||}|d|||d}|d||d||||d}t||||fS(s>convert central moments to mean, variance, skew, kurtosis iii(R( Rtmnctmnc2tmnc3tmnc4R R R R ((sTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pytmnc2mvsks *cCsdS(Ng?((tx((sTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pytmom_nc0'scCs|S(N((R((sTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pytmom_nc1*scCs||S(N((R((sTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pytmom_nc2-scCs |||S(N((R((sTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pytmom_nc30scCstj|dS(Ni(Rtpower(R((sTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pytmom_nc43sc sx{tD]r\}}tt|j|tdd\}}}}tj||||g}tj|}|jdkrtj |n|GHj }fd}t j |i} |t d|| } |d|| } |td|| } |td|| } |td|| }|d|dddd | }y(t| | | |f\}}}}WnYd | | | ||fGHtjgd\}}}}|jd krtj |qnXtj||||g}||||}||dd | fGHtj|d krTtjtj|d krTtj |n tj || ||||gt|QtRtlocitscaleit exceptionigMbP?(RtgetattrRtdictRtarraytisfinitetsumt distnonfinitetappendtexpectt specialcasestgetRRRRRtnantsizetdistextmaxtabstdistlowtdistoktres(tdistnametdistargstmtvtstktsttmaskR,t special_kwdstmnc0tmnc1RRRtmnc1_lctmetvetsetketemtdiff((RsTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pytcheck_cont_basic<s>' (3 c Cs"ddlm}d}|tjjddd d |d d tjdg}tjd dg}tjddg}tjdddg}d}|j}||tjddgj}||tjddgj} tjdddg} }t||| | f} dS(stest for incorrect kurtosis of nct D. Hogben, R. S. Pinkham, M. B. Wilk: The Moments of the Non-Central t-DistributionAuthor(s): Biometrika, Vol. 48, No. 3/4 (Dec., 1961), pp. 465-468 i(tassert_almost_equalgV`V?gso5?g*s?ga?i iRRtdecimalig8WV?g?gEd~?gl%?gu*'ȥ?gmn?g@iiiN(gV`V?gso5?g*s?ga?( t numpy.testingRJRtnctRR'titemR)R( RJt mvsk_10_1tc1tc2tc3tc4tnctmc1R R R tmvsk_nc((sTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pyt nct_kurt_bugqs( ##t__main__iigư>g{Gz?i(t SimpleTables Mean difference at least 1e-6theadersR7t diststatsR,s" Variance difference at least 1e-2s Skew difference at least 1e-2s" Kurtosis difference at least 1e-2(&t__doc__tstatsmodels.compat.pythonRtnumpyRtscipyRt)statsmodels.sandbox.distributions.sppatchRt distparamsRR-RRRRRRRR6R1R4R5R*RIRWt__name__R<R:R3R(tmean_tvar_R Rtstatsmodels.iolibRYtlen(((sTlib/python2.7/site-packages/statsmodels/sandbox/distributions/tests/check_moments.pytsP        5  eeee