作者asukaakimoto (唷唷~)
看板Statistics
标题[问题] 卡方同质性检定 达显着後 事後比较
时间Sun Jul 3 19:17:02 2011
是新手发问，不符合板规请告知，本篇会自d ^^"
想要请教是否有卡方 同质性检定 达显着後，事後检定的公式
之前google了很久都没有找到
之前补习班老师有给一个公式是
(Pi-Pj) +/- (X^2)^(1/2)*[ Pi*(1-Pi)/ni) + Pj*(1-Pj)/nj ]^(1/2)
X^2是 自由度为(R-1)*(C-1)，1-alpha 的临界值
(不知道会不会打得太乱看不太懂，等下会把公式用word打一次短网址上传^^")
http://www.hotimg.com/image/j3yNSu4
自己觉得公式有点奇怪，如果有人看过这个公式
不知是否可以帮忙回答下面的问题，谢谢^^
1."[ Pi*(1-Pi)/ni) + Pj*(1-Pj)/nj ]^(1/2)":
感觉这个式子应该是在估计信赖区间，而不是用在检定上@@
应该也要跟其他的事後检定一样，有虚无假设
H0:Pi=Pj，所以里面的式子应该要写成对共同母体比例的混何估计式
P=(Pi*ni+Pj*nj)/(ni+nj)
然後把原式改写成[P(1-P)(1/ni+1/nj)]^(1/2)这样@@
http://www.hotimg.com/image/FVzT7rF
2.有看到很多网路上的文章有写说应该要用Bonferroni修正
所以求卡方临界值时所用的alpha是维持原来的吗@@
该如何修正呢？
谢谢回答^^
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 220.136.177.39
1F：→ yhliu:既然是在整体性卡方检定显着後才做所谓 "事後比较", 就如 07/03 19:25
2F：→ yhliu:ANOVA 中整体 F 检定显着才做细部比较, 那麽, 没必要做什麽 07/03 19:26
3F：→ yhliu:Bonferroni 修正! 就像 ANOVA 中整体 F 检定显着後做的比较 07/03 19:27
4F：→ yhliu:用普通 t 临界值即可. 在同质性检定, k 群体卡方检定显着後 07/03 19:28
5F：→ yhliu:才做两两比较, 用普通两群体比较之检定程序即可. 07/03 19:29
6F：→ yhliu:整体的 type I error rate 已经由整体卡方检定保证了. 07/03 19:30
※ 编辑: asukaakimoto 来自: 220.136.177.39 (07/03 19:34)
7F：→ asukaakimoto:噢噢谢谢楼上，那1.的部分有要修正吗？还是用原来的 07/03 19:35
8F：→ asukaakimoto:式子呢？谢谢^^ 07/03 19:36
9F：→ asukaakimoto:btw,说要用Bonferroni好像是因为那些文章又作了 07/03 19:40
10F：→ asukaakimoto:列联表的subdivide，想请问这是否有关呢？然後好像就 07/03 19:41
11F：→ asukaakimoto:有要把p值乘上几倍(@@不好意思因为对原理不是很懂) 07/03 19:42
12F：→ asukaakimoto:所以也不知道p值乘上的倍数该怎麽计@@ 谢谢回答^^ 07/03 19:42
※ 编辑: asukaakimoto 来自: 220.136.177.39 (07/03 19:57)
13F：→ yhliu:补习班教的那个公式是仿 ANOVA 之 Scheffe 方法的. 那是 07/04 09:23
14F：→ yhliu:信赖区间公式, 也可以用该方法做检定, 也就是两比例差检定 07/04 09:24
15F：→ yhliu:之 z 临界值改用卡方临界值. 07/04 09:24
16F：→ yhliu:列联表的分割, 特别是分割/合并成一些 2 by 2 的表, 就像 07/04 09:25
17F：→ yhliu:ANOVA 中之对比(contrast). 07/04 09:26
18F：→ yhliu:如果先用整体卡方检定, 显着了才做那些细部比较, 就不需用 07/04 09:27
19F：→ yhliu:那仿 Scheffe 方法或 Bonferroni 方法; 如果没有先用卡方检 07/04 09:28
20F：→ yhliu:定控制整体型一误机率, 就应使用 Bonferroni 校正或仿 07/04 09:29
21F：→ yhliu:Scheffe 公式. 07/04 09:29
22F：→ yhliu:如果没有先用卡方检定做整体错误率控制, 而且要比较的 07/04 09:31
23F：→ yhliu:sub-tables 或 pairs of populations 是看过资料後才决定的, 07/04 09:31
24F：→ yhliu:那麽用 Bonferroni 方法仍是错的! 07/04 09:32
25F：→ asukaakimoto:噢噢，谢谢，所以原来的公式 "[ Pi*(1-Pi)/ni) + Pj* 07/04 16:20
26F：→ asukaakimoto:*(1-Pj)/nj ]^(1/2)"的公式是对的罗^^ 07/04 16:20
27F：→ asukaakimoto:所以如果要做这种事後比较，只要做原来的两比例t检定 07/04 16:25
28F：→ asukaakimoto:那临界t值的自由度该如何取呢?也是照列连表的方式吗? 07/04 16:28
29F：→ asukaakimoto:非常谢谢^^ 07/04 16:28