作者anovachen ( )
看板Statistics
标题Re: [问题] 费雪精确检定只限於2x2列联表？
时间Fri Oct 4 02:14:16 2013
考虑资料:
x y z
m 1 1 1
n 1 0 0
o 1 0 0
本来想用R练习看看，
可是查不到做3x3的Fisher exact test的语法，
所以改用SAS:
data test;
input V1 $ V2 $ count;
cards;
x m 1
x n 1
x o 1
y m 1
y n 0
y o 0
z m 1
z n 0
z o 0
;
run;
proc freq data=test;
weight count;
tables V1*V2/exact;
run;
Statistics for Table of V1 by V2
Statistic DF Value Prob
Chi-Square 4 2.2222 0.6950
Likelihood Ratio Chi-Square 4 2.9110 0.5728
Mantel-Haenszel Chi-Square 1 1.2656 0.2606
Phi Coefficient 0.6667
Contingency Coefficient 0.5547
Cramer's V 0.4714
WARNING: 100% of the cells have expected counts less
than 5. Chi-Square may not be a valid test.
Fisher's Exact Test
Table Probability (P) 0.3000
Pr <= P 1.0000
Sample Size = 5
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 111.255.26.96
1F：→ andrew43:r 内建的 fisher.test() 可以直接做 3*3 的 matrix 10/04 02:18
2F：→ andrew43:此例做出来也是 p = 1 10/04 02:20
3F：→ andrew43:fisher.test( matrix(c(1,1,1,1,0,0,1,0,0),3) ) 10/04 02:21
谢谢!!
---
p-value比较浅白的定义:虚无假设成立下抽取此样本及更极端样本的机率。
那这个范例中...我要怎麽手算出p-value阿QQ
2x2还比较好处理(维基百科也有范例)...
这个p-value=1要怎麽算出来??
※ 编辑: anovachen 来自: 111.255.26.96 (10/04 02:33)
4F：→ andrew43:查看看 Freeman-Halton extension 10/04 13:20
5F：→ andrew43:看到一个网页有列公式： 10/04 13:20
7F：→ yhliu:3x3, 固定行列和, 应有4个 cells 可变动, 要考虑的不是 10/04 20:44
8F：→ yhliu:bi-variate 吧? 10/04 20:44
感谢指正！
2x3才是bivariate，
我上次手算的答案应该是瞎猫碰上死耗子= ="
(因为第一列和第一行全都是1，其他格子都被固定住了，
但是如果(1,1)项改成其他数字，我的算法就不对了)
不负责任讲解:
3x3表格，参考该网页的Freeman-Halton extension:
a d g
b e h
c f i
p=[C(a+b+c,a)*C(b+c,b)]*[C(d+e+f,d)*C(e+f,e)]*[C(g+h+i,g)*C(h+i,h)]
/{n!/[(a+d+g)!(b+e+h)!(c+f+i)!]}
如果用高中所学的排列组合来理解，
分母就像"不尽相异物直线排列"，
而分子部分，固定第一列的数字之後，下面有两列可变动，
所以还要再算出第二列的可能排列方式。
e.g. 固定a+b+c，a的可能值有C(a+b+c,a)，
但是还要考虑a取完之後，b和c的变动情形，
所以要乘上C(b+c,b)，
而a,b都取完後，c已经无法再自由变动了，
所以C(c,c)=1，可忽略。
※ 编辑: anovachen 来自: 111.255.243.177 (10/08 07:25)