|
|
本帖最后由 Menuett 于 2013-12-22 15:59 编辑 " L* P% g8 N* A& J" \( }6 f
煮酒正熟 发表于 2013-12-20 12:05 ![]()
& l3 f) O, o0 x基本可以说是显著的。总的来说,在商界做统计学分析,95%信心水平是用得最多的,当95%上不显著时,都会去 ...
5 j' M6 C u% E4 X: s/ u5 s7 r% ~: Q2 i. B$ c# O4 J9 p
这个其实是一种binomial response,应该用Contigency Table或者Logisitic Regression(In case there are cofactors)来做。只记比率丢弃了Number of trial的信息(6841和1217个客户)。
! S' p6 N$ h, l$ Z/ f. v- w0 z" F) R) j! c! \5 Z
结果p=0.5731。 远远不显著。要在alpha level 0.05的水平上检验出76.42%和75.62%的区别,即使实验组和对照组各自样本大小相同,各自尚需44735个样本(At power level 80%)。see: Statistical Methods for Rates and Proportions by Joseph L. Fleiss (1981)! q1 P9 ?; ?" z# Q# f- k* T8 h+ c
9 T" t& g: h# d. a3 x$ QR example:
9 N/ Q$ f/ }: a$ H! m" E6 M
G6 A3 @" |+ ? o7 g0 v4 B* a> M<-as.table(rbind(c(1668,5173),c(287,930)))4 G3 k$ A3 _3 ^% U4 V" M& K+ w& M" p
> chisq.test(M)- F+ d$ u6 C! ~8 r& z: `
, l% o( v/ W" @ Pearson's Chi-squared test with Yates' continuity correction
# u( I# V* h) f, \" J5 G5 \ i& T' D$ e, {1 H& X& U7 A# j4 k
data: M, P' a/ J |3 ?+ W
X-squared = 0.3175, df = 1, p-value = 0.57310 n- v T5 W- x& A7 }# w# Z$ K7 ^8 s
. E0 \; a+ a' ]) ]
Python example:; G! S* K( O* Y& N E
3 D& X! V$ p5 R, {
>>> from scipy import stats4 q2 p6 ~. ] U* l. g* p
>>> stats.chi2_contingency([[6841-5173,5173],[1217-930,930]])( B N% ?7 `9 z7 F* \+ h; c0 S
(0.31748297614660292, 0.57312422493552839, 1, array([[ 1659.73628692, 5181.26371308],
# _7 r6 f9 J* ~6 ^& } [ 295.26371308, 921.73628692]])) |
|