作者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)