病例对照研究数据分析


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