U Kv f>@stdZddlZddlZddlmZddlmZddlm Z ddl m Z e Z ddd d d d d ddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;dd?d@dAdBg Created on Sun Jun 30 20:25:22 2013 Author: Josef Perktold N)assert_allclose) add_constant)Holder) TLinearModelzApr.1982zApr.1983zApr.1984zApr.1985zApr.1986zAug.1982zAug.1983zAug.1984zAug.1985zAug.1986zDec.1982zDec.1983zDec.1984zDec.1985zDec.1986zFeb.1284zFeb.1982zFeb.1983zFeb.1985zFeb.1986zJan.1982zJan.1983zJan.1984zJan.1985zJan.1986zJul.1982ZJuly1983ZJuly1984ZJuly1985ZJuly1986ZJune1982ZJune1983ZJune1984ZJune1985ZJune1986zMar.1982zMar.1983zMar.1984zMar.1985zMar.1986ZMay1982ZMay1983ZMay1984ZMay1985ZMay1986zNov.1982zNov.1983zNov.1984zNov.1985zNov.1986zOct.1982zOct.1983zOct.1984zOct.1985zOct.1986z Sept.1982z Sept.1983z Sept.1984z Sept.1985z Sept.1986gFxgr鷯gq= ףpg c?gH.!gZڊg\mgc]F?gʡEg(\µ?~jth?gvOjM?g0*gbX9ȶg1?gO@aӻ?g333333?gJY?g5;Nё?g:pΈg333333?g%ugX9vgݓZgX5;N?gjMStg%uggs?g٬?ggj+?gZӼㄿg'饿gUN@?g2U0*?gMJ?g-C6*?g"u?g+e?g5;NёgZd;O?g{Gz?gBig rh?g h"lxzg -?g_Qڋ?g?ܵ|g2%䃮g(\ſglxz,C?gw?ggsg ףp= ?g+e☿guV?gm{ggj+?gͪՖg_LegSt$?goŏ?glV}?gAf?g_vO?gΈ g,eXgȘ?gy&1l?gN@?gHPsgW2ı?g@߾g y)gj+ݓg6atolr rgh㈵>rtolrFrr"g;On?-C6 ?rHg{Gz?)res1res2rparamsloc_fitZ coefficientsbsetableZtvaluesZpvaluesZdofZdofsenpsqrtZ scale_fit fitted_valuesZmeanZllfZlogLik)selfrLrMZ scale_estrVL/tmp/pip-unpacked-wheel-2v6byqio/statsmodels/miscmodels/tests/test_tmodel.py test_basicIs***z!CheckTLinearModelMixin.test_basiccCs2|j}t|j|jdddt|j|jddddS)Ng333333?gMb`?rGg?gMbp?)rLrZbsejacrP)rUrLrVrVrWtest_bse[szCheckTLinearModelMixin.test_bsecCsH|j}|j}|}|jj|}t||jjddt||jjdddS)NgMb0?rJg>rE) rLrMZpredictmodelendogrrOrTZ residuals)rUrLrMZ fittedvaluesZresidrVrVrW test_fittedas  z"CheckTLinearModelMixin.test_fittedcCs`|j}|j}t|j|jddt|j|jddt|jj|jjddt|jj|jjdddS)Ng-C6?rErIrJg|=)rLresfrrNrPrZr[exog)rUrLr]rVrVrW test_formulams z#CheckTLinearModelMixin.test_formulacCsb|j}t|dd}|jjj\}}|j|||ks6t|j|dksHtt|j ||ks^tdS)Nk_extrarr ) rLgetattrrZr^shapeZdf_residAssertionErrorZdf_modellenrN)rUresr`ZnobsZk_varsrVrVrWtest_dfws  zCheckTLinearModelMixin.test_dfcCsX|j}|j}tt|j}||||||||||dSN rLr]rRZeyerdrNsummaryZt_testZf_testrUrLr]ZcontrrVrVrW test_smokes   z!CheckTLinearModelMixin.test_smokeN) __name__ __module__ __qualname__rXrYr\r_rfpytestmarksmokerkrVrVrVrWrCGs  rCc@seZdZeddZdS) TestTModelcCsvtj}ttj}t||}|jddd}tjdtjtjdd}|jddd}ddlm}||_ ||_ ||_ d |_ dS) NbfgsFmethodZdisp price ~ CRSPZpriceCRSP)datar ) res_t_dfestr) mm m_mariettarrxrfit from_formulaZresults_tmodelrzrMrLr]r`)clsr[r^modremodfr]rMrVrVrW setup_classs    zTestTModel.setup_classN)rlrmrn classmethodrrVrVrVrWrrsrrc@s(eZdZeddZejjddZdS)TestTModelFixedcCsjtj}ttj}t||dd}|jddd}tjdtjtjddd}|jddd}||_||_d |_ dS) Nr")fix_dfrsFrtrvrw)ryrr ) r{r|rrxrr}r~rLr]r`)rr[r^rrerr]rVrVrWrs  zTestTModelFixed.setup_classcCsX|j}|j}tt|j}||||||||||dSrgrhrjrVrVrWrks   zTestTModelFixed.test_smokeN) rlrmrnrrrorprqrkrVrVrVrWrs r)__doc__roZnumpyrRZ numpy.testingrZstatsmodels.tools.toolsrZstatsmodels.tools.testingrZstatsmodels.miscmodels.tmodelrr{Z date_labelarrayr|rxZam_candaterCrrrrVrVrVrWsv        H