病例對照研究數據分析


1,不匹配不分層

 

 第一,計算兩組暴露率有無統計學差異,卡方計算使用四個表專用公式,該校正校正

$X^{2}=\frac{\left ( ad-bc \right )^{2}n}{\left (a+b  \right )\left ( c+d \right )\left ( a+c \right )\left ( b+d \right )}$

X2=7.70,p<0.01,兩組暴露率有統計學差異。

 

第二, 計算OR值,ad/bc=2.20

 

第三,OR的CI(confidence interval)

Woolf對數轉換法

Var(lnOR)=1/a+1/b+1/c+1/d=0.0826

lnOR 95%CI=ln OR ±1.95* sqr(Var(lnOR))=ln2.2 ± 1.96 * 0.2874 = (0.2252, 1.3218)

OR 95% CI=(1.25,3.75)

Mietten卡方值法

$OR95perCI=OR^{\left ( 1\pm 1.96\sqrt{X^{2}} \right )}$ = (1.26,3.84)

 

2,不匹配分層

 

 

 第一,計算各層OR

OR1=(21*59)/(17*26)=2.80

OR2=(18*95)/(7*88)=2.78

兩層的OR均較不分層時OR值(2.20)大

 

第二,排除暴露因素與疾病的干擾

排除暴露的影響,OR=2.10,卡方=7.27,說明年齡與MI的發生有關聯

 

 

 排除疾病的干擾,OR=3.91,卡方=8.98,說明年齡與是否口服避孕葯也有聯系

 

 

 另外年齡也不是口服避孕葯與心肌梗死的中間環節。所以認為年齡是研究OC與MI關系時的混雜因素。可以用分層分析的方法來控制年齡的混雜作用。

用woolf齊性檢驗檢驗OR1與OR2是否同質(homogeneous)

 

第三,計算總的OR

若同質,計算總的OR,Mantel-Haenszel提出,ORMH

 

第四,計算總的卡方值

X2MH

 

第五,計算總的OR的可信區間

Mietten卡方值法

 

3,分級暴露

如能獲得暴露因素不同暴露水平的資料,可用來分子暴露和疾病的劑量反應關系

計算了各級OR和總的OR

 

第一,各級OR,逐級變大,呈現明顯的劑量-反應關系。

第二,計算趨勢卡方

Cochran-Armitage趨勢檢驗或者上文說到的MH卡方檢驗

4,匹配資料

介紹1:1配對

卡方用McNemar公式,OR=b/c,95%CI用Mietten法。

歸因分值類似AR,歸因危險度,≈(OR-1)/ OR

 

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM