|
|
本帖最后由 Menuett 于 2013-12-22 15:59 编辑
0 H3 r4 T' H3 z煮酒正熟 发表于 2013-12-20 12:05 1 q# S$ p; m4 v ]2 V: u2 l1 b/ m
基本可以说是显著的。总的来说,在商界做统计学分析,95%信心水平是用得最多的,当95%上不显著时,都会去 ... % t* o- [$ E! e7 f+ ~. k/ A" f9 A
+ I3 W' j D: l+ H, y# ]! O
这个其实是一种binomial response,应该用Contigency Table或者Logisitic Regression(In case there are cofactors)来做。只记比率丢弃了Number of trial的信息(6841和1217个客户)。 & Z2 }# `$ u$ l: j* V% C# @
) d: {. \1 Q9 E v6 C
结果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)
3 `! j0 T8 V y9 O' l8 G( f9 K' z, b) D* ?: m5 X
R example:
7 b# s4 k, ~! |9 W7 P8 u& T/ p7 Z+ E* p6 ~5 [7 y8 }2 o1 U8 ?
> M<-as.table(rbind(c(1668,5173),c(287,930)))
: P3 O% x& x0 Z8 y# M> chisq.test(M)# D4 P& M# P( H2 Z, A
c2 G6 P* d8 U: A5 g- P3 L7 E
Pearson's Chi-squared test with Yates' continuity correction& M7 r6 {, a% O
, W/ B5 R) M" R. Cdata: M
. a0 |' f9 W! x# P8 c; YX-squared = 0.3175, df = 1, p-value = 0.5731
9 |$ m% `4 D) ? Y' l: v: H o& V: X8 \+ x$ `% a9 C& ]" G
Python example:7 Y1 U& }; h! b1 G* U# y# n' K$ b& S
" t: e& }4 ~6 [* {" ` H( ?5 X: M>>> from scipy import stats5 l- G' z/ O! S
>>> stats.chi2_contingency([[6841-5173,5173],[1217-930,930]])
& f5 @: b; I! B3 V5 J, B7 S(0.31748297614660292, 0.57312422493552839, 1, array([[ 1659.73628692, 5181.26371308],
9 ?) o, F+ u, s6 c3 o [ 295.26371308, 921.73628692]])) |
|