ó áp7]c@slddlZddlmZmZddlmZmZmZm Z d„Z d„Z d„Z d„Z dS(iÿÿÿÿN(t assert_equaltassert_(t_calc_nodewise_rowt_calc_nodewise_weightt_calc_approx_inv_covtRegularizedInvCovariancecCsKtjjdƒtjjddƒ}t|ddƒ}t|jd ƒdS( NiA¤tsizei2iig{®Gáz„?i(i2i(i(tnptrandomtseedtnormalRRtshape(tXtghat((sRlib/python2.7/site-packages/statsmodels/stats/tests/test_regularized_covariance.pyttest_calc_nodewise_rowscCsftjjdƒtjjddƒ}tjjddƒ}t||ddƒ}tt|tƒƒdS( NiA¤Ri2iiig{®Gáz„?(i2i(RRR R RRt isinstancetfloat(R R tthat((sRlib/python2.7/site-packages/statsmodels/stats/tests/test_regularized_covariance.pyttest_calc_nodewise_weights cCs¾tjjdƒtjjddƒ}g}g}xUtdƒD]G}t||dƒ}t|||dƒ}|j|ƒ|j|ƒq>Wttj |ƒtj |ƒƒ}t |j dƒdS(NiA¤Ri2ig{®Gáz„?(i2i(ii( RRR R trangeRRtappendRtarrayRR (R tghat_ltthat_ltiR Rt theta_hat((sRlib/python2.7/site-packages/statsmodels/stats/tests/test_regularized_covariance.pyttest_calc_approx_inv_covs !cCsÐtjjdƒtjjdd ƒ}tjjtj|jƒƒ}td|ƒ}|j ƒtjj |j ƒ|ƒ}t |dkƒ|j ddƒt tj |j ƒd kƒtj |d kƒkƒdS( NiA¤Ri2itexoggš™™™™™¹?talphagà?i(i2i(RRR R tlinalgtinvtcovtTRtfittnormtapprox_inv_covRtsum(R Rtregcovtdiff((sRlib/python2.7/site-packages/statsmodels/stats/tests/test_regularized_covariance.pyttest_fit(s (tnumpyRt numpy.testingRRt(statsmodels.stats.regularized_covarianceRRRRRRRR'(((sRlib/python2.7/site-packages/statsmodels/stats/tests/test_regularized_covariance.pyts "