>h8\tSSKJr SSKrSSKrSSKJrJr SSKrSSK J r SSK J r J r Jr SSKJr SSKJr SrS rS rS rS rS rSrSr"SS5r"SS5r\S:XaPSSKrSSKJr \R@"SS/SQS9r!\"S\!S9RE5r#\"S\!S9RE5r$\"\#SS9r%gg))lrangeN) DataFrameIndex)stats)_has_intercept_intercept_idx_remove_intercept_patsy)summary2)OLScUcUR5$US:Xa UR$US:Xa UR$US:Xa UR$US:Xa UR$[ SU-5e)Nhc0hc1hc2hc3z robust options %s not understood) cov_paramscov_HC0cov_HC1cov_HC2cov_HC3 ValueError)modelrobusts jC:\Users\julio\OneDrive\Documentos\Trabajo\Ideas Frescas\venv\Lib\site-packages\statsmodels/stats/anova.py_get_covariancersk ~!! 5}} 5}} 5}} 5}};fDEEc 6URSS5nURSS5nURSS5nURSS5nU(aUR5nURRnURRnUR SnURR n URRRn URRn [U R5[U 5- S-n S U-n S S S X-/n[[R"U S 45US9nUS;a[!XXxXXX5 $US;a[#X XU U5$US;a[%X XU U5$US;a ['S5e[)S[+U5-5e)a Anova table for one fitted linear model. Parameters ---------- model : fitted linear model results instance A fitted linear model typ : int or str {1,2,3} or {"I","II","III"} Type of sum of squares to use. **kwargs** scale : float Estimate of variance, If None, will be estimated from the largest model. Default is None. test : str {"F", "Chisq", "Cp"} or None Test statistics to provide. Default is "F". Notes ----- Use of this function is discouraged. Use anova_lm instead. testFscaleNtyprrzPR(>%s)dfsum_sqmean_sqcolumnsr!I)II)III)IVzType IV not yet implementedzType %s not understood)getlowerrendogexogshape endog_namesdata design_info exog_nameslentermsrrnpzerosanova1_lm_singleanova2_lm_singleanova3_lm_singleNotImplementedErrorrstr)rkwargsrrr rr2r3nobs response_namer7r8n_rowspr_testnamestables r anova_singlerI#s. ::fc "D JJw %E **UA C ZZ$ 'F  KK  E ;;  D ::a=DKK++M++""..K''J+##$~k'BBQFF$G 8Y 6E bhh{+U ;E hd+ &g? ? F' &( ( F' &( ( !"?@@1CH<==rc [USS5n U cB[RRU5up[R"U R U5n [R "[UR5[UR545n URVs/sHoRU5PM nn[U5H unnSU UU4'M [R"XS-5n[U5nUU)n[R"UR5nUU)nUR5n[!US/-5Ul[R$U U)R'S5U4UR(USS/4'UR*UR,4UR(SSS/4'US:XaUSUS- UR*UR,- - XW'[.R0R3USUSUR,5XX'[R4[R44UR(SXx/4'USUS- US 'U$s snf) a Anova table for one fitted linear model. Parameters ---------- model : fitted linear model results instance A fitted linear model **kwargs** scale : float Estimate of variance, If None, will be estimated from the largest model. Default is None. test : str {"F", "Chisq", "Cp"} or None Test statistics to provide. Default is "F". Notes ----- Use of this function is discouraged. Use anova_lm instead. effectsNr!r*Residualr"r#rr$)getattrr;linalgqrdotTr<r9r: column_names term_namesslice enumeraterarraytolistrindexc_sumlocssrdf_residrfsfnan)rr2r3rCr7rHrErrFrrKqrarrnameslicesislice_r#idxrSrXs rr=r=_s .eY-Giill4 &&e$ ((C ))*C 0H0H,IJ KC2=2H2H I2H$%2HF If%&AvI&VVC! $F  %C SD\F+001JSD!J    E ,-EK)+sC4y}}Q/?/G)HEIIedH%%&-2YY-FEIIj8D/)* s{h%+5 ENN24 E#Jd $)NN413 *to-.Xt4E) L/JsI cPURSSn[U5nSSX4/n[[R"US45US9n[ US5n [ X5n /n /n [ U5GHupURU5n[URUR5n/n[UR5nUHn[UR5nURU5(dM0UU:XaM8URU5nUR[URUR55 UR[URUR55 M [R"UR R"R$S5Un[R"UR R"R$S5UnUR&(a[R("[R("UU 5UR*5nSSKJn UR1U5unnUR$SUR$S- n[R("USS2U*S24R*U5nOUnUR$SnUS :XabUR3UU S 9nUR4=UR6UR8U U4'nUR:UR6UR8U U4'UUR6UR8U S4'U R=UR5 U R=UR?55 GM [AU S /-5UlURB[RD"XR R"R$SS-/-5nXUS-URF-URH- nUUS'URFURH[RJ[RJ4UR6S SSX4/4'U$) aA Anova type II table for one fitted linear model. Parameters ---------- model : fitted linear model results instance A fitted linear model **kwargs** scale : float Estimate of variance, If None, will be estimated from the largest model. Default is None. test : str {"F", "Chisq", "Cp"} or None Test statistics to provide. Default is "F". Notes ----- Use of this function is discouraged. Use anova_lm instead. Type II Sum of Squares compares marginal contribution of terms. Thus, it is not particularly useful for models with significant interaction terms. Nr#r"r.r&r!r)rNrcov_prL)&r:r rr;r<rrUrTrstartstopsetfactorsissubsetextendeyerr3r4sizerPrQscipyrNrOf_testfvaluer[rXpvalueappendrdrilocargsortr\r]r`)rr7rErrFr terms_inforGrHcov robust_cov col_orderrXrftermcolsL1L2term_sett other_setcolLVLrN orth_compl_rbL12r^ test_valuer\s rr>r>sT2""1%J(4J tT +E bhh{+u =E % &C /JI EZ(  & DJJ * t||$AAIII  ++H 4I!''* &CHH56 &CHH56 VVEKK$$**1- .r 2 VVEKK$$**1- .r 2 77&&:.rtt4C $!99S>LJq bhhqk)A&&AqbcE*,,b1CC A 3; S 3A;<88 CEIIekk!nd* +j12EIIekk!ng- .+, %++a.$&'$ TYY[!O)R ,-EK JJrzz){{/?/?/E/Ea/H/J.K"KL ME +d #eii / >CE(O=BYYFFBFF=LEIIj8D$889 LrcU[U5- nURnSSX4/n[[R"US45US9n[ X5n /n /n [ U5GHupURU 5n[R"URRRS5UnUnURSnUS:XabURUU S9nUR=URURU U4'nUR URURU U4'UURURU S4'U R#U R%55 GM ['U S /-5UlXUS-UR(-UR*- nUUS'UR(UR*[R,[R,4URS SSX4/4'U$) Nr#r"r.r&r!rrrjrL)rr:rr;r<rrUrTrrrr3r4rurvr[rXrwrxrdrr\r]r`)rr7rErrFrr{rGrHr|r~rXrfrrrrrbr^rr\s rr?r?s n[))F""J tT +E bhh{+u =E % (CI EZ(  & VVEKK$$**1- .t 4 HHQK 3; S ,A;<88 CEIIekk!nd* +j12EIIekk!ng- .+, %++a.$&' TYY[!)" ,-EK +d #eii / >CE(O=BYYFFBFF=LEIIj8D$889 LrcURSS5n[U5S:XaUSn[U40UD6$US;a[S[ U5-5eURSS5nURSS 5n[U5nS U-nS S S SXG/n[ [ R"US45US9n U(dUSRnUV s/sHoRPM sn U S 'UV s/sHoRPM sn U S '[ R"U S R5*U RU RSS S 4'U S R5*U S'US:XamU SU S - U- U S'[R R#U SU S U S 5X'[ R$U RU SR'5U4'U $s sn fs sn f)a Anova table for one or more fitted linear models. Parameters ---------- args : fitted linear model results instance One or more fitted linear models scale : float Estimate of variance, If None, will be estimated from the largest model. Default is None. test : str {"F", "Chisq", "Cp"} or None Test statistics to provide. Default is "F". typ : str or int {"I","II","III"} or {1,2,3} The type of Anova test to perform. See notes. robust : {None, "hc0", "hc1", "hc2", "hc3"} Use heteroscedasticity-corrected coefficient covariance matrix. If robust covariance is desired, it is recommended to use `hc3`. Returns ------- anova : DataFrame When args is a single model, return is DataFrame with columns: sum_sq : float64 Sum of squares for model terms. df : float64 Degrees of freedom for model terms. F : float64 F statistic value for significance of adding model terms. PR(>F) : float64 P-value for significance of adding model terms. When args is multiple models, return is DataFrame with columns: df_resid : float64 Degrees of freedom of residuals in models. ssr : float64 Sum of squares of residuals in models. df_diff : float64 Degrees of freedom difference from previous model in args ss_dff : float64 Difference in ssr from previous model in args F : float64 F statistic comparing to previous model in args PR(>F): float64 P-value for significance comparing to previous model in args Notes ----- Model statistics are given in the order of args. Models must have been fit using the formula api. See Also -------- model_results.compare_f_test, model_results.compare_lm_test Examples -------- >>> import statsmodels.api as sm >>> from statsmodels.formula.api import ols >>> moore = sm.datasets.get_rdataset("Moore", "carData", cache=True) # load >>> data = moore.data >>> data = data.rename(columns={"partner.status" : ... "partner_status"}) # make name pythonic >>> moore_lm = ols('conformity ~ C(fcategory, Sum)*C(partner_status, Sum)', ... data=data).fit() >>> table = sm.stats.anova_lm(moore_lm, typ=2) # Type 2 Anova DataFrame >>> print(table) r r!rr(z6Multiple models only supported for type I. Got type %srrrNzPr(>%s)r]r\df_diffss_diffr&)r0r9rIrrArr;r<rr\r]diffvaluesr[rXrr^r_r`isnull) argsrBr rrrn_modelsrFrGrHmdls ranova_lmrsL **UA C 4yA~QE,V,, ('),S23 3 ::fc "D JJw %E4yH$G  9d DE bhh!}-u =E R'+,tGGt,E%L156#6E*-/WWU:5F5M5M-N,NEIIekk!"oy()e ))++E) s{9%i(885@c E#Ji0@$)*$5735&& %*##%w./ L-6s G $Gc\[R"S/U-5nUH nXnSX5'M U$)NTF)r;rV)reslices_to_excludenindrss r _not_slicers5 ((D6!8 C! L" Jrc[X$URS5nX5n[R"XSS2U4R U55nUR R U5n[ U5[ U5- nXx4$)a Residual sum of squares of OLS model excluding factors in `keys` Assumes x matrix is orthogonal Parameters ---------- y : array_like dependent variable x : array_like independent variables term_slices : a dict of slices term_slices[key] is a boolean array specifies the parameters associated with the factor `key` params : ndarray OLS solution of y = x * params keys : keys for term_slices factors to be excluded Returns ------- rss : float residual sum of squares df : int degrees of freedom r!N)rr4r;subtractrPrQr9) yx term_slicesparamskeysrparams1r\r]s r_ssr_reduced_modelrsi4 [ 3CkG ++a1c6w/ 0C %%))C.C1vG $H =rc8\rSrSrSrS SjrSrSrSrSr g) AnovaRMia: Repeated measures Anova using least squares regression The full model regression residual sum of squares is used to compare with the reduced model for calculating the within-subject effect sum of squares [1]. Currently, only fully balanced within-subject designs are supported. Calculation of between-subject effects and corrections for violation of sphericity are not yet implemented. Parameters ---------- data : DataFrame depvar : str The dependent variable in `data` subject : str Specify the subject id within : list[str] The within-subject factors between : list[str] The between-subject factors, this is not yet implemented aggregate_func : {None, 'mean', callable} If the data set contains more than a single observation per subject and cell of the specified model, this function will be used to aggregate the data before running the Anova. `None` (the default) will not perform any aggregation; 'mean' is s shortcut to `numpy.mean`. An exception will be raised if aggregation is required, but no aggregation function was specified. Returns ------- results : AnovaResults instance Raises ------ ValueError If the data need to be aggregated, but `aggregate_func` was not specified. Notes ----- This implementation currently only supports fully balanced designs. If the data contain more than one observation per subject and cell of the design, these observations need to be aggregated into a single observation before the Anova is calculated, either manually or by passing an aggregation function via the `aggregate_func` keyword argument. Note that if the input data set was not balanced before performing the aggregation, the implied heteroscedasticity of the data is ignored. References ---------- .. [*] Rutherford, Andrew. Anova and ANCOVA: a GLM approach. John Wiley & Sons, 2011. NcXlX lX@lSU;a [S5eXPlUb [ S5eX0lUS:Xa [RRUl OX`l URURU/U-S95(d+URbUR5 O Sn[U5eUR5 g)NCzSFactor name cannot be 'C'! This is in conflict with patsy's contrast function name.z)Between subject effect not yet supported!mean)subsetzThe data set contains more than one observation per subject and cell. Either aggregate the data manually, or pass the `aggregate_func` parameter.)r6depvarwithinrbetweenr@subjectpdSeriesraggregate_funcequalsdrop_duplicates _aggregate_check_data_balanced)selfr6rrrrrmsgs r__init__AnovaRM.__init__s   &=DE E  %'78 8 V #"$))..D "0 {{4//y67I/JKK"".!A!o% !!#rcURRUR/UR-SS9URR UR 5Ulg)NF)as_index)r6groupbyrrraggrrs rrAnovaRM._aggregatesQYYgt||nt{{:',..2kk;c$--. rcbSnURH,nU[URUR55-nM. 0n[ URR S5Hbn/nURH.nUR URURU5 M0 [U5nXS;a X5S-X5'M^SX5'Md Sn[U5U:wa [U5eUWnUHnXU:wdM [U5e URR SX-:a [S5eg)zraise if data is not balanced This raises a ValueError if the data is not balanced, and returns None if it is balance Return might change r!rzData is unbalanced.z9There are more than 1 element in a cell! Missing factors?N) rr9r6uniqueranger4rxrytupler) r factor_levelswi cell_countrXkeyr error_messagecounts rrAnovaRM._check_data_balanceds' ++B S2!5!5!78 8M 499??1-.EC{{ 499S>..u56#*C ",/A"5 "# /. z?m +]+ +3C3' // 99??1  5 5)* * 6rc URURRnURVs/sHnSU-PM nnSUR-nX4/-n[ R "SRU5URS9nURRnUHIn[R"S/URS-5n SXU'[R"U 5Xx'MK SRU5/n [XzURS5n USS2U 4n[X5n U R5n U R URS:a [#S 5eU HnUR%U5 M UH nXxU Xx'M U R&n U R(nU R*n/S Qn[,R."[R0"S 5US 9nUGHnURU;dMUS :wdM[3XX}U/5unnUU- nUU- U- nUSRUSS5:Xd US-U-U;aX- nUnO$[3XX}US-U-/5unnUU- nUU- U- nUU- n[4R6R9UUU5nUR;SS5R;SS5nUUR<US4'UUR<US4'UUR<US4'UUR<US4'GM [?U5$s snf)zVestimate the model and compute the Anova table Returns ------- AnovaResults instance z C(%s, Sum)*r6Fr!T:Nz$Independent variables are collinear.)F ValueNum DFDen DFPr > F)rr.r& InterceptrzC(z, Sum)rrrr) r6rrrrpatsydmatrixjoinr7term_name_slicesr;rVr4rr fitrankrpoprr]r\rrr<rrr^r_replacer[ AnovaResults)rrrfrrrorrrr term_excluderresultsrr]r\r' anova_tablessr1 df_resid1df1msmmsedf2rprs rr AnovaRM.fit&s IIdkk " ) )-1KK8Kq,"K8-9$ MM#((7+$)) <mm44 C((E71771:-.C$(CC !!xx}K )* AGGAJ? afIA ))+ :: "CD DA OOA C*/4K ##kk;ll288F#3WE C||3&3++="4+u#6i(*czS(388GCRL11sW,K?.C"C&8ksW,-'/OD)$h.C#:,C#IGGJJq#s+{{4,44XrB34 i025 h/25 h/23 h//2K((m9sL)rrr6rrr)NNN) __name__ __module__ __qualname____firstlineno____doc__rrrr__static_attributes__rrrrs(5nDH $$<0 *B@)rrc*\rSrSrSrSrSrSrSrg)riizD Anova results class Attributes ---------- anova_table : DataFrame cXlgNr)rrs rrAnovaResults.__init__qs&rc>UR5R5$r)summary__str__rs rrAnovaResults.__str__ts||~%%''rc[R"5nURS5 URUR5 U$)zLcreate summary results Returns ------- summary : summary2.Summary instance Anova)r Summary add_titleadd_dfr)rsumms rrAnovaResults.summaryws5! w D$$% rrN) rrrrrrrrrrrrrris'( rr__main__)olsz moore.csvr!)partner_status conformity fcategoryfscore)skiprowsrGz5conformity ~ C(fcategory, Sum)*C(partner_status, Sum)rz#conformity ~ C(partner_status, Sum)r*)r )&statsmodels.compat.pythonrnumpyr;pandasrrrrrtr statsmodels.formula.formulatoolsrrr statsmodels.iolibr #statsmodels.regression.linear_modelr rrIr=r>r?rrrrrrstatsmodels.formula.apirread_csvmoorermoore_lmmooreBrHrrrrs,#  '3 F"9>x4nVp%NgTD})})@8 z+ OOK!#9 :EJ # 6U C G G IF X1 %E/r