病例对照研究数据分析
1,不匹配不分層
第一,計(jì)算兩組暴露率有無統(tǒng)計(jì)學(xué)差異,卡方計(jì)算使用四個(gè)表專用公式,該校正校正
$X^{2}=frac{left ( ad-bc ight )^{2}n}{left (a+b ight )left ( c+d ight )left ( a+c ight )left ( b+d ight )}$
X2=7.70,p<0.01,兩組暴露率有統(tǒng)計(jì)學(xué)差異。
第二, 計(jì)算OR值,ad/bc=2.20
第三,OR的CI(confidence interval)
①Woolf對數(shù)轉(zhuǎn)換法
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 ( 1pm 1.96sqrt{X^{2}} ight )}$ = (1.26,3.84)
2,不匹配分層
第一,計(jì)算各層OR
OR1=(21*59)/(17*26)=2.80
OR2=(18*95)/(7*88)=2.78
兩層的OR均較不分層時(shí)OR值(2.20)大
第二,排除暴露因素與疾病的干擾
排除暴露的影響,OR=2.10,卡方=7.27,說明年齡與MI的發(fā)生有關(guān)聯(lián)
排除疾病的干擾,OR=3.91,卡方=8.98,說明年齡與是否口服避孕藥也有聯(lián)系
另外年齡也不是口服避孕藥與心肌梗死的中間環(huán)節(jié)。所以認(rèn)為年齡是研究OC與MI關(guān)系時(shí)的混雜因素。可以用分層分析的方法來控制年齡的混雜作用。
用woolf齊性檢驗(yàn)檢驗(yàn)OR1與OR2是否同質(zhì)(homogeneous)
第三,計(jì)算總的OR
若同質(zhì),計(jì)算總的OR,Mantel-Haenszel提出,ORMH
第四,計(jì)算總的卡方值
X2MH
第五,計(jì)算總的OR的可信區(qū)間
Mietten卡方值法
3,分級暴露
如能獲得暴露因素不同暴露水平的資料,可用來分子暴露和疾病的劑量反應(yīng)關(guān)系
計(jì)算了各級OR和總的OR
第一,各級OR,逐級變大,呈現(xiàn)明顯的劑量-反應(yīng)關(guān)系。
第二,計(jì)算趨勢卡方
Cochran-Armitage趨勢檢驗(yàn)或者上文說到的MH卡方檢驗(yàn)
4,匹配資料
介紹1:1配對
卡方用McNemar公式,OR=b/c,95%CI用Mietten法。
歸因分值類似AR,歸因危險(xiǎn)度,≈(OR-1)/ OR
Valar morghulis
總結(jié)
以上是生活随笔為你收集整理的病例对照研究数据分析的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: volatile的适用场景
- 下一篇: Python 读写文件中w与wt, r与