ó î&]\c@`szddlmZmZmZddlmZddlZddlm Z m Z m Z m Z m Z ddlZddlmZd„Zd„Zejjd„ƒZejjd „ƒZd „Zd „Zd ejfd „ƒYZdejfd„ƒYZdejfd„ƒYZdejfd„ƒYZd„Zd„Zd„Z d„Z!d„Z"d„Z#d„Z$d„Z%dS(i(tdivisiontprint_functiontabsolute_import(tstatsN(tassert_almost_equaltassert_tassert_array_almost_equaltassert_array_almost_equal_nulptassert_allclose(traisesc C`sÃtjjdƒd}tjj|ƒ}|jƒ}|jddƒ}tj|ƒ}tjdddƒ}|j |ƒ}tj j |d|d |ƒ}|d|d }t tj ||d ƒ|d kƒ|j|tjƒ} |jtj |ƒ} t| d ddƒt| d ddƒt|j|tjƒ| ddƒt|jtj |ƒ| ddƒt|j|ƒ|d j ƒ|dd ƒt|j||d ƒ||j ƒ|dd ƒdS(NiîÀ…iôtddofiiùÿÿÿiiõtloctscaleiig{®Gáz„?gà?tdecimali (tnptrandomtseedtrandntmeantstdRt gaussian_kdetlinspacetevaluatetnormtpdfRtsumtintegrate_box_1dtinfRt integrate_boxt integrate_kdetintegrate_gaussian( t n_basesampletxntxnmeantxnstdtgkdetxstkdepdftnormpdft intervalltprob1tprob2((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyt test_kde_1d s* %"#c C`søtjjdƒd}tjj|ƒ}tjj|ƒ}tj|d|ƒ}tjtj||dd|ƒƒ}tj|d|ƒ}tj dddƒ}|j |ƒ}tj j |d|d |ƒ}|d |d } t tj||dƒ| d kƒ|j|tjƒ} |jtj |ƒ} t| d dd ƒt| d dd ƒt|j|tjƒ| ddƒt|jtj |ƒ| ddƒt|j|ƒ|djƒ| ddƒt|j||dƒ||jƒ| ddƒdS(NiîÀ…iôtweightsiiùÿÿÿiiõR R iig{®Gáz„?gà?R i (RRRRtrandtaveragetsqrtRRRRRRRRRRRRRR( RR twnR!R"R#R$R%R&R'R(R)((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_kde_1d_weighted*s,&%"#cC`s3tjjdƒd}tjddgƒ}tjddgddggƒ}tjj||d|ƒj}tj|ƒ}tjdd d …dd d …f\}}tj |j ƒ|j ƒgƒ}|j |ƒ}|j ddƒ}tjj tj||gƒd |d |ƒ} |j ƒd |j ƒd} ttj|| dƒ| ddkƒd} d} |j| |d g| | gƒ} |j| | g| |d gƒ}t| ddd ƒt|ddd ƒt|j|ƒ|djƒ| dddƒt|j||ƒ|| jƒ| dddƒdS(NiîÀ…iôgð?g@g@g@tsizeiùÿÿÿiy@@Rtcoviiig{®Gáz„?g}Ô%­I²Ôg}Ô%­I²Tgà?R (RRRtarraytmultivariate_normaltTRRtmgridtvstacktravelRtreshapeRtdstackRRRRRR(RRt covarianceR R#txtyt grid_coordsR%R&R'tsmalltlargeR(R)((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyt test_kde_2dJs.!+!-)""cC`sKtjjdƒd}tjddgƒ}tjddgddggƒ}tjj||d|ƒj}tjj|ƒ}tj|d|ƒ}tj d d d …d d d …f\}}tj |j ƒ|j ƒgƒ}|j |ƒ} | j ddƒ} tjjtj||gƒd |d |ƒ} |j ƒd|j ƒd} ttj| | dƒ| ddkƒd} d} |j| |dg| | gƒ}|j| | g| |dgƒ}t|dddƒt|dddƒt|j|ƒ| djƒ| dddƒt|j||ƒ| | jƒ| dddƒdS(NiîÀ…iôgð?g@g@g@R1R+iùÿÿÿiy@@RR2iiig{®Gáz„?g}Ô%­I²Ôg}Ô%­I²Tgà?R (RRRR3R4R5R,RRR6R7R8RR9RR:RRRRRR(RRR;R R/R#R<R=R>R%R&R'R?R@R(R)((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_kde_2d_weightedqs0!+!-)""c C`sæd„}tjjdƒd}tjj|ƒ}tj|ƒ}tj|d|ƒ}tj|d|jƒ}tjdddƒ}|j|ƒ}|j|ƒ}t ||ƒ|j|ƒ} t || ƒt t tj|ddƒdS( NcS`stj|jd|jdƒS(s*Same as default, just check that it works.gð¿i(Rtpowertntd(tkde_obj((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyt scotts_factoršsiîÀ…i2t bw_methodiùÿÿÿii3t wrongstring( RRRRRRtfactorRRRt assert_raisest ValueError( RGRR R#tgkde2tgkde3R$R%tkdepdf2tkdepdf3((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_kde_bandwidth_method™s   c C`sæd„}tjjdƒd}tjj|ƒ}tj|ƒ}tj|d|ƒ}tj|d|jƒ}tjdddƒ}|j|ƒ}|j|ƒ}t ||ƒ|j|ƒ} t || ƒt t tj|ddƒdS( NcS`stj|jd|jdƒS(s*Same as default, just check that it works.gð¿i(RRCtneffRE(RF((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyRG´siîÀ…i2RHiùÿÿÿii3RI( RRRRRRRJRRRRKRL( RGRR R#RMRNR$R%RORP((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyt"test_kde_bandwidth_method_weighted³s   t_kde_subclass1cB`seZd„ZRS(cC`sDtj|ƒ|_|jj\|_|_|j|_|jƒdS(N( Rt atleast_2dtdatasettshapeRERDRGtcovariance_factort_compute_covariance(tselfRV((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyt__init__Òs (t__name__t __module__R[(((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyRTÑst_kde_subclass2cB`seZd„ZRS(cC`s&|j|_tt|ƒj|ƒdS(N(RGRXtsuperR^R[(RZRV((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyR[Ús (R\R]R[(((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyR^Ùst_kde_subclass3cB`seZd„Zd„ZRS(cC`s ||_tjj||ƒdS(N(R;RRR[(RZRVR;((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyR[às cC`sOtjj|jƒ|_tjtjjdtj|jƒƒ|j|_ dS(Ni( RtlinalgtinvR;tinv_covR.tdettpiRDt _norm_factor(RZ((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyRYäs&(R\R]R[RY(((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyR`ßs t_kde_subclass4cB`seZd„ZRS(cC`sd|jƒS(Ngà?(tsilverman_factor(RZ((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyRXës(R\R]RX(((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyRgêsc`sctjdddddgdtƒ}tjddd d ƒ}tj|ƒ‰ˆ|ƒ}t|ƒ}||ƒ}t||d dƒt|ƒ}||ƒ}t||d dƒt |ˆj ƒ}||ƒ}t||d dƒt |ƒ} | |ƒ} d d dddg} t | | ddƒˆ} ‡fd†| _ | jƒ| |ƒ} t|| d dƒdS(Niùÿÿÿiûÿÿÿiiitdtypeiöÿÿÿi tnumi2tnulpgœn—,°?gÇ»e‡?°?gWÍ[3®?g¥aúà)¶?gKd8ú;´?R ic`sˆjS(N(RJ((tkde(s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyts(RR3tfloatRRRRTRR^R`R;RgRRXRY(tx1R$tystkde1ty1tkde2ty2tkde3ty3tkde4ty4t y_expectedtkde5ty5((Rls<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_gaussian_kde_subclassingïs,$          cC`s£tjdddddgdtƒ}tjddd dƒ}d d d d dg}tj|ƒ}|jddƒ|jddƒ||ƒ}t||ddƒdS(NiùÿÿÿiûÿÿÿiiiRiiöÿÿÿi Rjg>L\ž9™?gp4=¨?gùHŒ–꟫?g8(H MT«?gIJ™CR ‘?RHgà?tscottR i(RR3RnRRRt set_bandwidthR(RoR$RyRlRt((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyt$test_gaussian_kde_covariance_cachings$ cC`s¥tjdddddgdtƒ}tjddd d ƒ}tj|ƒ}|j|_|jƒ||ƒ}tj|d d ƒ}||ƒ}t ||d dƒdS(sÞUgly, but people may rely on this. See scipy pull request 123, specifically the linked ML thread "Width of the Gaussian in stats.kde". If it is necessary to break this later on, that is to be discussed on ML. iùÿÿÿiûÿÿÿiiiRiiöÿÿÿi Rji2RHt silvermanRkN( RR3RnRRRRhRXRYR(RoR$RlRrRsRt((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_gaussian_kde_monkeypatch%s$    cC`sPtjdƒ}tj|ƒ}dddddg}t||ƒ|ddƒdS(sRegression test for #1181.igB…×\AÁ?g­±ø ESÇ?gÔö f§úÈ?R iN(RtarangeRRR(RoRlRy((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_kde_integer_input:scC`sûtjjdƒd}tjj|ƒ}tj|ƒ}tjdddƒ}|j|ƒ}|j|ƒ}t ||ddƒtj |ƒ}|j |ƒ}t ||ddƒtj|ƒ}tj |j|ƒƒ}|j |ƒ}t ||ddƒdS(Nii2iñÿÿÿi iR ( RRRRRRRRRRtlogtlogpdf(RR R#R$Rtpdf2R…tlogpdf2((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_pdf_logpdfBsc C`stjjdƒd}tjj|ƒ}tjj|ƒ}tj|d|ƒ}tjdddƒ}|j|ƒ}|j |ƒ}t ||ddƒtj |ƒ}|j |ƒ}t ||ddƒtj|ƒ}tj |j|ƒƒ}|j |ƒ}t ||ddƒdS(Nii2R+iñÿÿÿi iR ( RRRRR,RRRRRRR„R…( RR R/R#R$RR†R…R‡((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_pdf_logpdf_weightedZs cC`sŽtjjdƒtjjddƒ}tjjdddgd|jƒ}|jƒ}tjtj |ƒd|ƒt ||ddd dƒdS( Ni90R1idgð?g$@R+tatolg›+¡†›„=trtol( RRRt lognormaltchoiceR1tcopyRRtlog10R(tvalsR+t orig_weights((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_weights_intactss $ cC`s«tjjdƒdddddg}ddd d d g}tj|d |ƒ}tj|d tj|ƒƒ}d ddg}t|j|ƒ|j|ƒddddƒdS(Ni90gš™™™™™É?g+@g5@gÀR@gÀX@iiiiiR+g333333Ó?i iXRŠg›+¡†›„=R‹(RRRRRtfloat64RR(tvaluesR+tpdf_itpdf_fR ((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyttest_weights_integer~s(&t __future__RRRtscipyRtnumpyRt numpy.testingRRRRRtpytestR RKR*R0tmarktslowRARBRQRSRRTR^R`RgR|RRRƒRˆR‰R’R—(((s<lib/python2.7/site-packages/scipy/stats/tests/test_kdeoth.pyts. (   '(    &