ó î&]\c@`sÆddlmZmZmZddlZddlmZmZm Z m Z m Z ddl m ZddlmZddlmZmZmZd„Zd„Zd „Zd „Zd „Zd „ZdS( i(tdivisiontprint_functiontabsolute_importN(t assert_equaltassert_array_equaltassert_array_almost_equaltassert_approx_equaltassert_allclose(traises(txlogy(tmarginst expected_freqtchi2_contingencycC`s tjdgƒ}t|ƒ}tt|ƒdƒ|d}t|tjdgƒƒtjdggƒ}t|ƒ\}}tjdggƒ}tjdggƒ}t||ƒt||ƒtjdƒjddƒ}t|ƒ\}}tjdgdggƒ}tjddd dd d ggƒ}t||ƒt||ƒtjd ƒjdd dƒ}t|ƒ\}}}tjdggdgggƒ}tjdgdgdgggƒ}tjddddgggƒ}t||ƒt||ƒt||ƒdS(Niii iiii3ii iiiiiiBiÒi<i\i|iHiN(tnptarrayR RtlenRtarangetreshape(tatmtm0tm1t expected0t expected1tm2t expected2((sAlib/python2.7/site-packages/scipy/stats/tests/test_contingency.pyt test_margins s2    $  !$!  cC`sýttdgƒtjdgƒƒtjddgddggddgddggddgddgggƒ}t|ƒ}t|tj|ƒƒtjdddgdddggƒ}t|ƒ}tjdddgd d d ggƒ}t||ƒdS( Nigð?iii ig(@g0@g2@g8@(RR R Rt ones_likeR(tobservedtetcorrect((sAlib/python2.7/site-packages/scipy/stats/tests/test_contingency.pyttest_expected_freq+s"N ' 'cC`sátjddgddggƒ}t|dtƒ\}}}}t|dƒt|dƒt|dƒt||ƒtjdddgƒ}t|dtƒ\}}}}t|dƒt|dƒt|dƒt||ƒdS(Niit correctionggð?ii(R RR tFalseRR(tobstchi2tptdoftexpected((sAlib/python2.7/site-packages/scipy/stats/tests/test_contingency.pyttest_chi2_contingency_trivial8s!       cC`sd}d}tjdddgdddgdd d gddd ggd d dgdddgdddgd ddgggƒ}t|ƒ\}}}}t|dddƒt|ddd ƒt|dƒd}d}tjddgddggddgddgggddgddggd dgddggggƒ}t|ƒ\}}}}t|ddd ƒt|d dd ƒt|dƒdS(!Ns  # Data vector. data <- c( 12, 34, 23, 4, 47, 11, 35, 31, 11, 34, 10, 18, 12, 32, 9, 18, 13, 19, 12, 12, 14, 9, 33, 25 ) # Create factor tags:r=rows, c=columns, t=tiers r <- factor(gl(4, 2*3, 2*3*4, labels=c("r1", "r2", "r3", "r4"))) c <- factor(gl(3, 1, 2*3*4, labels=c("c1", "c2", "c3"))) t <- factor(gl(2, 3, 2*3*4, labels=c("t1", "t2"))) # 3-way Chi squared test of independence s = summary(xtabs(data~r+c+t)) print(s) sÏ Call: xtabs(formula = data ~ r + c + t) Number of cases in table: 478 Number of factors: 3 Test for independence of all factors: Chisq = 102.17, df = 17, p-value = 3.514e-14 i i"ii#ii i i iii/i ii ii!ig{®GáŠY@t significantig"„j,5È#=isb # Data vector. data <- c( # 12, 17, 11, 16, # 11, 12, 15, 16, # 23, 15, 30, 22, # 14, 17, 15, 16 ) # Create factor tags:r=rows, c=columns, d=depths(?), t=tiers r <- factor(gl(2, 2, 2*2*2*2, labels=c("r1", "r2"))) c <- factor(gl(2, 1, 2*2*2*2, labels=c("c1", "c2"))) d <- factor(gl(2, 4, 2*2*2*2, labels=c("d1", "d2"))) t <- factor(gl(2, 8, 2*2*2*2, labels=c("t1", "t2"))) # 4-way Chi squared test of independence s = summary(xtabs(data~r+c+d+t)) print(s) sÏ Call: xtabs(formula = data ~ r + c + d + t) Number of cases in table: 262 Number of factors: 4 Test for independence of all factors: Chisq = 8.758, df = 11, p-value = 0.6442 iiiigj¼t“„!@g“©‚QIä?(R RR RR(tRcodetRoutputR"R#R$R%R&((sAlib/python2.7/site-packages/scipy/stats/tests/test_contingency.pyttest_chi2_contingency_RLs<             cC`sCtjddgddggƒ}t|dddtƒ\}}}}t|dt|||ƒjƒƒt|dddtƒ\}}}}|tjdd gd dggƒ}t|dt|||ƒjƒƒtjd d d gd d d ggƒ}t|ddƒ\}}}}t|dt|||ƒjƒƒdS( Nii<iZtlambda_slog-likelihoodR igà¿gà?i i (R RR R!RR tsumtTrue(tctgR$R%Rtc_corr((sAlib/python2.7/site-packages/scipy/stats/tests/test_contingency.pyttest_chi2_contingency_g«s!$$$%$'cC`s…tjddgddggƒ}ttt|ƒtjddgddggƒ}ttt|ƒtjdƒ}ttt|ƒdS(Niÿÿÿÿi iiii(ii(R Rt assert_raisest ValueErrorR tempty(R"((sAlib/python2.7/site-packages/scipy/stats/tests/test_contingency.pyttest_chi2_contingency_bad_args¹s !!(t __future__RRRtnumpyR t numpy.testingRRRRRtpytestRR3t scipy.specialR tscipy.stats.contingencyR R R RRR'R+R2R6(((sAlib/python2.7/site-packages/scipy/stats/tests/test_contingency.pyts (   _