K層四格表OR值的一致性檢驗(yàn)方法模擬比較研究.pdf

1、研究背景:
　　 流行病學(xué)研究中，為了排除混雜因素對(duì)疾病與暴露之間關(guān)聯(lián)性的影響，常常需要將資料進(jìn)行分層分析。當(dāng)分組變量和結(jié)局變量均為二分類，混雜變量為K個(gè)水平時(shí)，觀察到的數(shù)據(jù)可以K×2×2表的形式展現(xiàn)。但是，在評(píng)估疾病與暴露因素之間是否有關(guān)聯(lián)以及在估計(jì)公共比數(shù)比(oddsratio，OR)之前，需要對(duì)各層之間的比數(shù)比進(jìn)行一致性檢驗(yàn)。
　　 關(guān)于K×2×2表資料的一致性檢驗(yàn)，眾多學(xué)者提出了不同方法，根據(jù)現(xiàn)已公開發(fā)表的文獻(xiàn)可

2、以概括為11種。本研究按求解公共OR值的計(jì)算方法類型以及數(shù)據(jù)結(jié)構(gòu)特點(diǎn)，將這11種一致性檢驗(yàn)方法劃分為如下6大類:
　　 ①M(fèi)antel-Haenszel方法(簡(jiǎn)稱M-H):基于Mantel和Haenszel(1959)求解公共OR值的方法，主要有Breslow-Day法和Breslow-Day-Tarone法兩種，后者是對(duì)前者方法的校正。
　　 ②Zelen精確概率方法(簡(jiǎn)稱exact):不用先求解公共OR值，這種一致性

3、檢驗(yàn)方法是在Fisher精確概率方法基礎(chǔ)上的一種拓展方法。
　　 ③近似非條件最大似然(AU，Asymptotic Unconditional)方法:基于近似非條件最大似然估計(jì)方法求解公共OR值，有Wald法、Score法及似然比檢驗(yàn)法三種方法。
　　 ④加權(quán)最小二乘(WLS，Weighted Least Square)方法:基于加權(quán)最小二乘方法求解公共OR值的一致性檢驗(yàn)方法。
　　 ⑤近似條件最大似然(AC，A

4、symptotic Conditional)方法:基于超幾何分布用近似條件最大似然估計(jì)方法求解公共OR值，有Zelen近似法和Score法兩種方法。
　　 ⑥稀疏數(shù)據(jù)(SD，sparsedata)檢驗(yàn)方法:基于近似條件最大似然估計(jì)方法求解公共OR值，并且結(jié)合中心極限定理或混合線性模型理論，主要有Liang和Self(1985)提出的T4和T5兩種方法。
　　 在劃分為上述6類的基礎(chǔ)上，作者期望以此為切入點(diǎn)能夠從不同于前人

5、的角度對(duì)K×2×2表資料一致性檢驗(yàn)方法進(jìn)行更加獨(dú)特的模擬比較，以做出更為合理細(xì)致的解釋。
　　 研究目的:
　　 本研究旨在通過Monte Carlo方法，模擬比較不同一致性檢驗(yàn)方法的檢驗(yàn)效能1—β和I型錯(cuò)誤α，以期望回答如下幾個(gè)問題:
　　 1.針對(duì)不同的數(shù)據(jù)類型，哪種檢驗(yàn)方法更適合于稀疏數(shù)據(jù)，哪種方法對(duì)于非稀疏數(shù)據(jù)更值得采用?
　　 2.Zelen精確概率法優(yōu)劣性如何，是否比其他近似方法更為“精確”?

6、
　　 3.基于超幾何分布而求解得到公共OR值的檢驗(yàn)方法是否比其他方法更值得提倡?
　　 4.廣泛應(yīng)用于Meta分析異質(zhì)性檢驗(yàn)的WLS方法，在默認(rèn)的檢驗(yàn)水準(zhǔn)0.05時(shí)，其檢驗(yàn)效能以及I型錯(cuò)誤與其他檢驗(yàn)方法比較如何?
　　 5.Meta分析異質(zhì)性檢驗(yàn)時(shí)，我們常常將檢驗(yàn)水準(zhǔn)設(shè)置為0.10，因此當(dāng)將檢驗(yàn)水準(zhǔn)設(shè)定為0.10時(shí)，這11種檢驗(yàn)方法的檢驗(yàn)效能及I型錯(cuò)誤如何，模擬比較結(jié)果與檢驗(yàn)水準(zhǔn)為0.05時(shí)比較是否發(fā)生較大改變

7、?
　　 方法:
　　 模擬過程考慮4個(gè)參數(shù)，即第i層樣本量Ni、ORl值、對(duì)照組暴露概率π0i和層數(shù)K，每個(gè)參數(shù)分別取不同值以此構(gòu)成多種參數(shù)組合，每種組合模擬運(yùn)算1000次。借助SAS9.2(R)ranbin隨機(jī)數(shù)函數(shù)在給定種子數(shù)、樣本量Ni、對(duì)照組暴露概率π0i情況下產(chǎn)生第i層對(duì)照組暴露數(shù);病例組暴露概率由π1i=OR(l)π0i/ORiπ0i+(1-π0i)計(jì)算而得，同理可以產(chǎn)生第i層病例組暴露數(shù)。當(dāng)進(jìn)行I型錯(cuò)誤α

8、比較時(shí)，各層間ORi值恒相等，先用隨機(jī)數(shù)函數(shù)rannor產(chǎn)生均數(shù)為logOR，方差為0的隨機(jī)數(shù)，然后再將logORi指數(shù)化;當(dāng)進(jìn)行檢驗(yàn)效能1-β比較時(shí)，K層ORl值來(lái)自不同的總體，先用隨機(jī)數(shù)函數(shù)rannor產(chǎn)生均數(shù)為logOR，方差為1的隨機(jī)數(shù)，然后再將logORl指數(shù)化。不管是比較(I)型錯(cuò)誤α還是檢驗(yàn)效能1-β時(shí)，對(duì)照組暴露概率π0i由兩種方式產(chǎn)生，一種是為固定相等的常數(shù)，還有一種是由均勻分布產(chǎn)生。另外，鑒于一致性檢驗(yàn)方法常用于Me

9、ta分析中，我們?cè)谶M(jìn)行Meta分析的效應(yīng)量異質(zhì)性檢驗(yàn)時(shí)，檢驗(yàn)水準(zhǔn)常取0.10，因此本研究在進(jìn)行模擬比較時(shí)，設(shè)定兩個(gè)檢驗(yàn)水準(zhǔn)，分別為0.05和0.10。
　　 結(jié)果:
　　 在控制I型錯(cuò)誤方面:(1)針對(duì)不同的數(shù)據(jù)類型:近似非條件Wald法更適合于稀疏數(shù)據(jù);非稀疏數(shù)據(jù)類型，Breslow-Day法、Breslow-Day-Tarone法、近似非條件Score法以及近似非條件似然比檢驗(yàn)法較為合適。(2)Zelen精確概率法過

10、于保守，并且隨著樣本量的增加以及層數(shù)逐漸增大時(shí)，精確概率法的運(yùn)算量以及所需的運(yùn)算時(shí)間也隨之迅速增加。(3)基于超幾何分布求解得到公共OR值的檢驗(yàn)方法，在控制I型錯(cuò)誤方面相對(duì)差于其他檢驗(yàn)方法。(4)廣泛應(yīng)用于Meta分析的WLS法，不管檢驗(yàn)水準(zhǔn)取值為0.05還是0.10時(shí)，與M-H以及AU方法相比均過于保守。(5)大多數(shù)情況下，檢驗(yàn)水準(zhǔn)設(shè)定為0.05和0.10時(shí)，這11種方法間的比較并不因?yàn)闄z驗(yàn)水準(zhǔn)的不同而發(fā)生較大的改變。但是在某些特定情

11、形下，檢驗(yàn)水準(zhǔn)設(shè)定為0.10甚至稍顯更為合適。
　　 從檢驗(yàn)效能方面來(lái)看:(1)這11種方法均不適用于稀疏數(shù)據(jù)情形，當(dāng)樣本量較大層數(shù)也較大時(shí)，M-H方法、AU方法以及AC方法均較為合適。(2)Zelen精確概率法不值得提倡，其原因在于不僅檢驗(yàn)效能低，而且需要的運(yùn)算量、運(yùn)算時(shí)間也較大。(3)基于超幾何分布求解得到公共OR值的檢驗(yàn)方法，AC方法檢驗(yàn)效能尚可，但是SD方法檢驗(yàn)效能較差。(4)廣泛應(yīng)用于Meta分析的WLS法，當(dāng)樣本量以

12、及層數(shù)較大時(shí)值得采用。(5)當(dāng)檢驗(yàn)水準(zhǔn)設(shè)定為0.10時(shí)，這11種檢驗(yàn)方法的檢驗(yàn)效能間的比較，較0.05時(shí)并未發(fā)生較大的改變。
　　 結(jié)論:
　　 (1)當(dāng)樣本量較小時(shí)，尚且沒有合適的檢驗(yàn)方法值得采用。
　　 (2)SD方法不提倡使用。
　　 (3)Zelen精確概率法過于保守，不推薦使用。
　　 (4)近似非條件Score法以及近似非條件似然比檢驗(yàn)法最為合適，其次為M-H方法。
　　 (5

13、)M-H方法由于計(jì)算的便捷性，在樣本量進(jìn)一步增大，層數(shù)進(jìn)一步增加時(shí)，可以視為最佳的選擇。Breslow-Day-Tarone法比Breslow-Day法更為保守，所以相對(duì)來(lái)說，本文不建議對(duì)Breslow-Day法進(jìn)行校正。
　　 (6)AU方法與AC方法相比，綜合I型錯(cuò)誤以及檢驗(yàn)效能來(lái)說，AU方法優(yōu)于AC方法。
　　 (7)WLS方法相對(duì)于M-H以及AU方法來(lái)說偏于保守，但是當(dāng)樣本量很大時(shí)，由于其計(jì)算的便捷性，因此，被廣