兩個率或多個率的比較

兩個率或多個率的比較第1頁，課件共25頁，創(chuàng)作于2023年2月率的抽樣誤差與可信區(qū)間

一、率的抽樣誤差與標(biāo)準(zhǔn)誤二、總體率的可信區(qū)間第2頁，課件共25頁，創(chuàng)作于2023年2月一、率的抽樣誤差與標(biāo)準(zhǔn)誤

樣本率(p)和總體率(π)的差異稱為率的抽樣誤差(samplingerrorofrate)，用率的標(biāo)準(zhǔn)誤（standarderrorofrate）度量。如果總體率π未知，用樣本率p估計第3頁，課件共25頁，創(chuàng)作于2023年2月標(biāo)準(zhǔn)誤的計算第4頁，課件共25頁，創(chuàng)作于2023年2月二、總體率的可信區(qū)間

總體率的可信區(qū)間(confidenceintervalofrate)：根據(jù)樣本率推算總體率可能所在的范圍

第5頁，課件共25頁，創(chuàng)作于2023年2月率的統(tǒng)計學(xué)推斷

一、樣本率與總體率比較u檢驗二、兩個樣本率的比較u檢驗第6頁，課件共25頁，創(chuàng)作于2023年2月一、樣本率與總體率比較的u檢驗u檢驗的條件：np和n（1-p）均大于5時第7頁，課件共25頁，創(chuàng)作于2023年2月二、兩個獨(dú)立樣本率比較的u檢驗表5-1兩種療法的心血管病病死率比較療法死亡生存合計病死率(%)鹽酸苯乙雙胍26(X1)178204(n1)12.75(p1)安慰劑2(X2)6264(n2)3.13(p2)合計2824026810.45(pc)u檢驗的條件：n1p1和n1(1-p1)與n2p2和n2(1-p2)均>5第8頁，課件共25頁，創(chuàng)作于2023年2月卡方檢驗

χ2檢驗(Chi-squaretest)是現(xiàn)代統(tǒng)計學(xué)的創(chuàng)始人之一，英國人K.Pearson（1857-1936）于1900年提出的一種具有廣泛用途的統(tǒng)計方法，可用于兩個或多個率間的比較，計數(shù)資料的關(guān)聯(lián)度分析，擬合優(yōu)度檢驗等等。本章僅限于介紹兩個和多個率或構(gòu)成比比較的χ2檢驗。第9頁，課件共25頁，創(chuàng)作于2023年2月一、卡方檢驗的基本思想(1)療法死亡生存合計病死率(%)鹽酸苯乙雙胍26(a)178(b)204(a+b)12.75(p1)安慰劑2(c)62(d)64(c+d)3.13(p2)合計28(a+c.)240(b+d.)268(a+b+c+d=n)10.45(pc)表5-1兩種療法的心血管病病死率的比較2×2表或四格表(fourfoldtable)實際頻數(shù)A

(actualfrequency)(a、b、c、d)的理論頻數(shù)T(theoreticalfrequency)（H0:π1=π2=π）：a的理論頻數(shù)＝(a+b)×pc=(a+b)×[(a+c.)/n]=nRnC/n=21.3b的理論頻數(shù)＝(a+b)×(1-pc)=(a+b)×[(b+d.)/n]=nRnC/n=182.7c的理論頻數(shù)＝(c+d)×pc=(c+d)×[(a+c)/n]=nRnC/n=6.7d的理論頻數(shù)＝(c+d)×(1-pc)=(c+d)×[(b+d.)/n]=nRnC/n=57.3第10頁，課件共25頁，創(chuàng)作于2023年2月一、卡方檢驗的基本思想(2)各種情形下，理論與實際偏離的總和即為卡方值（chi-squarevalue）,它服從自由度為ν的卡方分布。第11頁，課件共25頁，創(chuàng)作于2023年2月3.847.8112.59P＝0.05的臨界值χ2分布（chi-squaredistribution）第12頁，課件共25頁，創(chuàng)作于2023年2月χ2檢驗的基本公式上述基本公式由Pearson提出，因此軟件上常稱這種檢驗為Peareson卡方檢驗，下面將要介紹的其他卡方檢驗公式都是在此基礎(chǔ)上發(fā)展起來的。它不僅適用于四格表資料，也適用于其它的“行×列表”。第13頁，課件共25頁，創(chuàng)作于2023年2月二、四格表專用公式（1）

為了不計算理論頻數(shù)T,可由基本公式推導(dǎo)出，直接由各格子的實際頻數(shù)（a、b、c、d）計算卡方值的公式：第14頁，課件共25頁，創(chuàng)作于2023年2月二、四格表專用公式（2）2(1)～u2＝2.19492＝4.82（n>40，所有T5時）第15頁，課件共25頁，創(chuàng)作于2023年2月三、連續(xù)性校正公式（1）

χ2分布是一連續(xù)型分布，而行×列表資料屬離散型分布，對其進(jìn)行校正稱為連續(xù)性校正(correctionforcontinuity),又稱Yates校正（Yates'correction）。⑴當(dāng)n≥40，而1≤T＜5時，用連續(xù)性校正公式⑵當(dāng)n＜40或T＜1時，用Fisher精確檢驗(Fisherexacttest)校正公式：第16頁，課件共25頁，創(chuàng)作于2023年2月三、連續(xù)性校正公式（2）因為1＜T＜5，且n＞40時，所以應(yīng)用連續(xù)性校正χ2檢驗第17頁，課件共25頁，創(chuàng)作于2023年2月四、配對四格表資料的χ2檢驗第18頁，課件共25頁，創(chuàng)作于2023年2月配對四格表資料的χ2檢驗也稱McNemar檢驗（McNemar'stest）H0：b，c來自同一個實驗總體（兩種劑量的毒性無差異）；H1：b，c來自不同的實驗總體（兩種劑量的毒性有差別）；α=0.05。第19頁，課件共25頁，創(chuàng)作于2023年2月五、行×列（R×C）表資料的χ2檢驗第20頁，課件共25頁，創(chuàng)作于2023年2月R×C表的χ2檢驗通用公式第21頁，課件共25頁，創(chuàng)作于2023年2月幾種R×C表的檢驗假設(shè)H0第22頁，課件共25頁，創(chuàng)作于2023年2月R×C表的計算舉例第23頁，課件共25頁，創(chuàng)作于2023年2月R×C表χ2檢驗的應(yīng)用注意事項1.對R×C表，若較多格子（1/5）的理論頻數(shù)小于5或有一個格子的理論頻數(shù)小于1，則易犯第一類錯誤。 出現(xiàn)某些格子中理論頻數(shù)過小時怎么辦？（1）增大樣本含量（最好?。?）刪去該格所在的行或列（丟失信息?。?）根據(jù)專業(yè)知識將該格所在行或列與別的行或列合并。（丟失信息！甚至出假象）

第24頁，課件共25頁，創(chuàng)作于2023年2