[問題] 卡方同質性檢定 達顯著後 事後比較

看板Statistics作者 (唷唷~)時間13年前 (2011/07/03 19:17), 編輯推噓0(0推 0噓 29→)
留言29則, 2人參與, 最新討論串1/1
是新手發問，不符合板規請告知，本篇會自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
→
07/03 19:25, , 1F
既然是在整體性卡方檢定顯著後才做所謂 "事後比較", 就如
07/03 19:25, 1F
→
07/03 19:26, , 2F
ANOVA 中整體 F 檢定顯著才做細部比較, 那麼, 沒必要做什麼
07/03 19:26, 2F
→
07/03 19:27, , 3F
Bonferroni 修正! 就像 ANOVA 中整體 F 檢定顯著後做的比較
07/03 19:27, 3F
→
07/03 19:28, , 4F
用普通 t 臨界值即可. 在同質性檢定, k 群體卡方檢定顯著後
07/03 19:28, 4F
→
07/03 19:29, , 5F
才做兩兩比較, 用普通兩群體比較之檢定程序即可.
07/03 19:29, 5F
→
07/03 19:30, , 6F
整體的 type I error rate 已經由整體卡方檢定保證了.
07/03 19:30, 6F
※ 編輯: asukaakimoto 來自: 220.136.177.39 (07/03 19:34)
→
07/03 19:35, , 7F
噢噢謝謝樓上，那1.的部分有要修正嗎？還是用原來的
07/03 19:35, 7F
→
07/03 19:36, , 8F
式子呢？謝謝^^
07/03 19:36, 8F
→
07/03 19:40, , 9F
btw,說要用Bonferroni好像是因為那些文章又作了
07/03 19:40, 9F
→
07/03 19:41, , 10F
列聯表的subdivide，想請問這是否有關呢？然後好像就
07/03 19:41, 10F
→
07/03 19:42, , 11F
有要把p值乘上幾倍(@@不好意思因為對原理不是很懂)
07/03 19:42, 11F
→
07/03 19:42, , 12F
所以也不知道p值乘上的倍數該怎麼計@@ 謝謝回答^^
07/03 19:42, 12F
※ 編輯: asukaakimoto 來自: 220.136.177.39 (07/03 19:57)
→
07/04 09:23, , 13F
補習班教的那個公式是仿 ANOVA 之 Scheffe 方法的. 那是
07/04 09:23, 13F
→
07/04 09:24, , 14F
信賴區間公式, 也可以用該方法做檢定, 也就是兩比例差檢定
07/04 09:24, 14F
→
07/04 09:24, , 15F
之 z 臨界值改用卡方臨界值.
07/04 09:24, 15F
→
07/04 09:25, , 16F
列聯表的分割, 特別是分割/合併成一些 2 by 2 的表, 就像
07/04 09:25, 16F
→
07/04 09:26, , 17F
ANOVA 中之對比(contrast).
07/04 09:26, 17F
→
07/04 09:27, , 18F
如果先用整體卡方檢定, 顯著了才做那些細部比較, 就不需用
07/04 09:27, 18F
→
07/04 09:28, , 19F
那仿 Scheffe 方法或 Bonferroni 方法; 如果沒有先用卡方檢
07/04 09:28, 19F
→
07/04 09:29, , 20F
定控制整體型一誤機率, 就應使用 Bonferroni 校正或仿
07/04 09:29, 20F
→
07/04 09:29, , 21F
Scheffe 公式.
07/04 09:29, 21F
→
07/04 09:31, , 22F
如果沒有先用卡方檢定做整體錯誤率控制, 而且要比較的
07/04 09:31, 22F
→
07/04 09:31, , 23F
sub-tables 或 pairs of populations 是看過資料後才決定的,
07/04 09:31, 23F
→
07/04 09:32, , 24F
那麼用 Bonferroni 方法仍是錯的!
07/04 09:32, 24F
→
07/04 16:20, , 25F
噢噢，謝謝，所以原來的公式 "[ Pi*(1-Pi)/ni) + Pj*
07/04 16:20, 25F
→
07/04 16:20, , 26F
*(1-Pj)/nj ]^(1/2)"的公式是對的囉^^
07/04 16:20, 26F
→
07/04 16:25, , 27F
所以如果要做這種事後比較，只要做原來的兩比例t檢定
07/04 16:25, 27F
→
07/04 16:28, , 28F
那臨界t值的自由度該如何取呢?也是照列連表的方式嗎?
07/04 16:28, 28F
→
07/04 16:28, , 29F
非常謝謝^^
07/04 16:28, 29F
‣ 更多 asukaakimoto 的文章
文章代碼(AID): #1E44-oDv (Statistics)