《衛(wèi)生統(tǒng)計學(xué)》2觀察性研究的設(shè)計與實施

1、郭艷郭艷 中山大學(xué)公共衛(wèi)生學(xué)院中山大學(xué)公共衛(wèi)生學(xué)院 2021-7-221 l 不對研究對象施加任何干預(yù)措施的情況下，不對研究對象施加任何干預(yù)措施的情況下， 通過觀察或訪問的方法，客觀地記錄被研通過觀察或訪問的方法，客觀地記錄被研 究事物的狀況。究事物的狀況。 也稱非實驗研究也稱非實驗研究 (non- experimental study) l 描述疾病或治療模式的現(xiàn)狀描述疾病或治療模式的現(xiàn)狀 評估某種暴露對健康的影響評估某種暴露對健康的影響 監(jiān)測某種轉(zhuǎn)歸（如藥物不良反應(yīng)）監(jiān)測某種轉(zhuǎn)歸（如藥物不良反應(yīng)） 探索病因和預(yù)防措施探索病因和預(yù)防措施 2021-7-222 l描述性研究（描述性研究（des

2、criptive study） 最常見：橫斷面研究最常見：橫斷面研究 l分析性研究（分析性研究（analytical study） 病例病例-對照研究（對照研究（case-control study） - 回顧性回顧性 隊列研究（隊列研究（cohort study） - 前瞻性前瞻性 2021-7-223 2021-7-224 （一）單純隨機抽樣調(diào)查資料的分析（一）單純隨機抽樣調(diào)查資料的分析 例例16-1 從某中學(xué)從某中學(xué)5400名學(xué)生中抽取名學(xué)生中抽取270人作為樣本人作為樣本 進行近視眼情況調(diào)查，以估計該中學(xué)學(xué)生的平均視進行近視眼情況調(diào)查，以估計該中學(xué)學(xué)生的平均視 力及患近視眼學(xué)生的比例。

3、力及患近視眼學(xué)生的比例。 1均數(shù)的區(qū)間估計均數(shù)的區(qū)間估計 2概率的區(qū)間估計概率的區(qū)間估計 2021-7-225 表表 16-1 單純隨機抽樣標準誤及區(qū)間估計計算公式單純隨機抽樣標準誤及區(qū)間估計計算公式 總體總體 類型類型 總體均數(shù)總體均數(shù) 標準誤標準誤 區(qū)間估計區(qū)間估計 總體概率總體概率 標準誤標準誤 區(qū)間估計區(qū)間估計 無限無限 總體總體 n S S X X n StX 1, 2/ 1 1 n PP SP )( P SZP 2/ 有限有限 總體總體 n S N nN S X 1 X n StX 1, 2/ 1 1 1 n PP N nN SP )( P SZP 2/ 注：注： 有限總體校正因子

4、有限總體校正因子 1 N nN n 較大時，可用較大時，可用 2/ Z代替代替 1, 2/n t 2021-7-226 等比例分配等比例分配 (proportional allocation) 大層多抽，小層少抽；大層多抽，小層少抽； 各層抽取的比例與該層在總體中所占比例相同各層抽取的比例與該層在總體中所占比例相同 N N n n ii 非等比例分配非等比例分配 (non- proportional allocation) 各層抽取相同數(shù)目的個體各層抽取相同數(shù)目的個體 最優(yōu)分配分層隨機抽樣（最優(yōu)分配分層隨機抽樣（optimum allocation） 大層多抽，小層少抽；大層多抽，小層少抽；

5、變異大多抽，變異小少抽變異大多抽，變異小少抽 分層抽樣資料分析方法：加權(quán)平均！分層抽樣資料分析方法：加權(quán)平均！ 2021-7-22 7 2021-7-228 l成組設(shè)計成組設(shè)計：比較兩組的暴露史：比較兩組的暴露史 病例組病例組 對照組對照組 分層成組設(shè)計分層成組設(shè)計： 病例組病例組 層層1 層層2 層層k 對照組對照組 層層1 層層2 層層k l匹配設(shè)計匹配設(shè)計 病例病例: 對照對照 = 1 : 1 ; 病例病例: 對照對照 = 1 : m l病例與對照在非處理因素方面具有可比性病例與對照在非處理因素方面具有可比性 2021-7-229 2021-7-2210 例例16-7 食管癌與飲酒的病例

6、食管癌與飲酒的病例-對照研究對照研究 每日飲酒量每日飲酒量(g/天天) 組別組別 80+ 0-79 合計合計 病例組病例組 96 (a) 104 (b) 200 (a+b) 對照組對照組 109 (c) 666 (d) 775 (c+d) 合計合計 205 (a+c) 770 (b+d) 975 (1) 優(yōu)勢比的估計優(yōu)勢比的估計 645 109104 66696 . bc ad OR 暴露，暴露，2水平水平 2021-7-2211 （2）優(yōu)勢比的假設(shè)檢驗優(yōu)勢比的假設(shè)檢驗 0 H：OR1， 1 H：OR 1 檢驗統(tǒng)計量 14.110 770205775200 )10910466696)(1975

7、( )()()( )(1( 2 2 2 dcdbcaba bcadN 注：注：檢驗統(tǒng)計量分子上檢驗統(tǒng)計量分子上:（N-1)， 而以前所學(xué)而以前所學(xué)Pearson 2 統(tǒng)計量的分子上：統(tǒng)計量的分子上：N N 較大時，兩者等價較大時，兩者等價 2021-7-2212 Woolf 法法 dcba ORVar 1111 )(ln lnOR： )(ln.ln ORVarOR961 OR： )(ln.ORVar eOR 961 例例 16-7 03070 770 1 205 1 775 1 200 1 .ln ）（ORVar OR： )95. 7 ,00. 4(64. 5 0307. 096. 1 )(l

8、n96. 1 e eOR ORVar 2021-7-2213 Miettinen 法法 OR： 2 96. 1 1 OR 例例 16-7 64. 5,14.110 2 OR OR： )79. 7 ,08. 4(64. 5 14.110 96. 1 1 96. 1 1 2 OR 2021-7-2214 例例16-8 食管癌與飲酒關(guān)系病例食管癌與飲酒關(guān)系病例-對照研究對照研究 年齡組年齡組（歲）（歲） 2544 4554 5564 65+ 合計合計 病病 對對 病病 對對 病病 對對 病病 對對 病病 對對 每日飲酒量每日飲酒量 (g/天天) 例例 照照 例例 照照 例例 照照 例例 照照 例例

9、照照 80+ 5 35 25 29 42 27 24 18 96 109 079 5 270 21 138 34 139 44 119 104 666 hOR 7.71 5.67 6.36 3.61 5.64 hhh a dn 4.2857 16.1972 24.1240 13.9317 58.5386 hhh b ch 0.5556 2.8592 3.7934 3.8634 11.0716 Th 1.2698 11.6620 21.6694 13.9317 48.5329 Vh 1.0768 6.8576 10.6701 7.4392 26.0437 （1）公共優(yōu)勢比公共優(yōu)勢比 1 1 58

10、.5386 5.29 11.0716 H hh h h H hh h h a d n OR b c n 對象分層對象分層 2021-7-2215 （2）對公共優(yōu)勢比對公共優(yōu)勢比 OR 的的假設(shè)檢驗假設(shè)檢驗 0 H：OR1， 1 H：OR 1 檢驗檢驗統(tǒng)計量統(tǒng)計量 H h h H h h H h h MH V Ta 1 2 112 )( Th：對應(yīng)于：對應(yīng)于 ah的的理論頻數(shù)理論頻數(shù) Vh：對應(yīng)于：對應(yīng)于 ah的方差的方差 1010 3 hhhh h hh n n m m V nn 例例 16-8 22 2 ()9648.5329 86.51 26.0437 hh MH h aT V （） 20

11、21-7-2216 （3）公共優(yōu)勢比公共優(yōu)勢比 OR 的的 95%置信區(qū)間置信區(qū)間 Miettinen 法估計：法估計： 2 1.96 1 MH OR 例例 16-8 OR： 2 1.96 1 MH OR = 51.86 96. 1 1 29. 5= (3.72 , 7.51) 2021-7-2217 例例 16-9 食管癌與飲酒關(guān)系食管癌與飲酒關(guān)系病例病例-對照對照研究研究 每日飲酒量每日飲酒量(g/天天) 0-39 40-79 80-119 120+ 合計合計 Xk 0 1 2 3 病例數(shù)病例數(shù) k a(Tk) 29(85.13) 75(72.82) 51(28.31) 45(13.74)

12、 200( 1 n) 對照數(shù)對照數(shù) k b 386 280 87 22 775(n0) 合計人數(shù)合計人數(shù) k m 415 355 138 67 975(n) 優(yōu)勢比優(yōu)勢比 OR 1.0 3.57 7.80 27.23 H0:不存在劑量不存在劑量-反應(yīng)線性關(guān)系反應(yīng)線性關(guān)系； H1:存在劑量存在劑量-反應(yīng)線性關(guān)系反應(yīng)線性關(guān)系 檢驗統(tǒng)計量檢驗統(tǒng)計量 23 2 0 22 10 00 () () () K kkk k KK kkkk kk XaTnn n n nm Xm X =152.99 df=1 暴露分暴露分 多個水平多個水平 2021-7-2218 例例 16-10 子宮內(nèi)膜癌與服用雌激素關(guān)系子宮

13、內(nèi)膜癌與服用雌激素關(guān)系的的病例病例-對照研對照研 究究 正常正常 (對照對照) 服用雌激素服用雌激素 未未服服 合計合計 服用雌激素服用雌激素 27 對對 29 對對 56 對對 子宮內(nèi)膜癌子宮內(nèi)膜癌 (病例病例) 未未服服 3 對對 4 對對 7 對對 合合 計計 30 對對 33 對對 63 對對 2021-7-2219 例例 16-11 子宮內(nèi)膜癌與服用雌激素關(guān)系的子宮內(nèi)膜癌與服用雌激素關(guān)系的 1:4 病例病例-對照研究對照研究 正常正常 (對照對照) 服用服用雌激素雌激素 未服未服 合計合計 1 2 3 4 0 服用雌激素服用雌激素 17 組組 16 組組 15 組組 5 組組 3 組

14、組 56 組組 子宮內(nèi)膜癌子宮內(nèi)膜癌 (病例病例) 未服未服 4 組組 1 組組 1 組組 1 組組 0 組組 7 組組 63 組組 病例病例：63 名（其中名（其中 56 名服用名服用雌激素雌激素） ；） ； 對照：對照：634 = 252 名名 匹配組：匹配組：63 個個 暴露分暴露分2個水平個水平; 1:4小組中服用小組中服用 人數(shù)不同人數(shù)不同 2021-7-2220 2021-7-2221 暴露暴露 Exposure 暴露人群暴露人群 Cohort 危險因素危險因素 Risk factor 暴露結(jié)局暴露結(jié)局 Outcome 對照人群對照人群 Cohort 危險因素危險因素 Risk f

15、actor 對照結(jié)局對照結(jié)局 Outcome 2021-7-2222 例例16-15 為研究鼻咽癌與為研究鼻咽癌與EB病毒關(guān)系，某市進行病毒關(guān)系，某市進行 了了EB病毒抽樣調(diào)查，共檢查病毒抽樣調(diào)查，共檢查25-70歲居民歲居民41471人，人， 其中其中38550人人EB病毒抗體陰性，病毒抗體陰性，2921人陽性。以人陽性。以 EB病毒抗體陰性和陽性者形成兩個組，進行了為病毒抗體陰性和陽性者形成兩個組，進行了為 期期12年的隊列研究。年的隊列研究。 隨訪人群隨訪人群 發(fā)病數(shù)發(fā)病數(shù) 暴露人年數(shù)暴露人年數(shù) 人年發(fā)病率（人年發(fā)病率（/10 萬）萬） 抗體陽性抗體陽性 57 20477 278.36

16、抗體陰性抗體陰性 54 437579 12.34 合計合計 111 458056 24.23 56.22 34.12 36.278 RR 人年資料人年資料 暴露組暴露組 2021-7-2223 隨訪人群隨訪人群 發(fā)病人數(shù)發(fā)病人數(shù) 觀察人年數(shù)觀察人年數(shù) 暴露組暴露組 1 d 1 T 非暴露組非暴露組 0 d 0 T 合計合計 D T 0 H：RR=1， 1 H：RR1 Mantel-Haenszel 2 統(tǒng)計量統(tǒng)計量 2 211 10 () MH d TT D T DT 0 H成立時，近似地服從自由度為成立時，近似地服從自由度為 1 的的 2 分布分布 人年人年 資料：資料： 2 211 10

17、() MH d TT D T DT = 2 (5745805620477 111) 20477 111 437579 =571.25 拒絕零假設(shè)，可以認為，拒絕零假設(shè)，可以認為，RR 1 2021-7-2224 RR 的（的（1-）置信區(qū)間：）置信區(qū)間： 2 2/ 1 MH Z RR AR 的（的（1-）置信區(qū)間：）置信區(qū)間： 2 2/ 1 MH Z AR 相對危險度相對危險度RR=風險之比風險之比 歸因危險度歸因危險度AR=風險之差風險之差 RR 的的 95%置信區(qū)間：置信區(qū)間： 47571 961 1 1 5622 2 2 . . . / MH Z RR AR 的的 95%置信區(qū)間：置信區(qū)

18、間： ) 47.571 96. 1 1(02.266 1 2 2/ MH Z AR 2021-7-2225 隨訪人群隨訪人群 發(fā)病人數(shù)發(fā)病人數(shù) 未發(fā)病人數(shù)未發(fā)病人數(shù) 總?cè)藬?shù)總?cè)藬?shù) 暴露組暴露組 a c a+c 非暴露組非暴露組 b d b+d 頻率資料：頻率資料： 2 2 (1)() ()()()() MH Nadbc ab ac cd bd 2021-7-2226 例例16-16 按年齡分層的人時資料按年齡分層的人時資料 H h h hh H h h hh T Tb T Ta RR 1 1 1 0 = 28 46149/64399785 19898/68712 22 18250/643999

19、5 48814/68712 =3.36 按年齡分層按年齡分層 50 50 未按年齡分層未按年齡分層 危險因素危險因素 危險因素危險因素 危險因素危險因素 指標指標 有有(1) 無無(0) 小計小計 有有 無無 小計小計 有有 無無 小計小計 發(fā)病人數(shù)發(fā)病人數(shù) 28(a) 22(b) 50(m) 785 95 880 813 122 935 觀察人年數(shù)觀察人年數(shù) 18250(T1) 46149(T0) 64399(T) 48814 19898 68712 67064 66047 133111 RR 3.22 3.38 6.56 （1）RR 的估計的估計 每層有一個每層有一個 四格表四格表 202

20、1-7-2227 H0：RR=1 H1：RR1 Mantel-Haenszel 2 統(tǒng)計量統(tǒng)計量 按年齡分層按年齡分層 50 50 未按年齡分層未按年齡分層 危險因素危險因素 危險因素危險因素 危險因素危險因素 指標指標 有有(1) 無無(0) 小計小計 有有 無無 小計小計 有有 無無 小計小計 發(fā)病人數(shù)發(fā)病人數(shù) 28(a) 22(b) 50(m1) 785 95 880 813 122 935 觀察人年數(shù)觀察人年數(shù) 18250(T1) 46149(T0) 64399(T) 48814 19898 68712 67064 66047 133111 RR 3.22 3.38 6.56 （2）R

21、R 的檢驗的檢驗 64.156 ) 68712 1989848814880 64399 461491825050 ( ) 68712 48814880 64399 1825050 ()78528( 22 2 1 2 011 2 1 11 1 2 H h h hhh H h h hh H h h MH T TTm T Tm a 理論頻數(shù)理論頻數(shù) 理論頻數(shù)理論頻數(shù) 實際頻數(shù)實際頻數(shù) 28 （3）RR 的的95%置信區(qū)間置信區(qū)間 （Miettinen法）法） 64.156 96. 1 1 1 36. 3 2 2/ MH Z RR =（2.74，4.10） 注：頻數(shù)資料的分層分析注：頻數(shù)資料的分層分析

22、 第第 h 層四格表一般形式層四格表一般形式 危險因素危險因素 發(fā)病與否發(fā)病與否 有有 無無 小計小計 發(fā)病發(fā)病 ah bh m1h 未發(fā)病未發(fā)病 ch dh m0h 小計小計 n1h n0h nh H h h hh H h h hh n nb n na RR 1 1 1 0 H h hh hhhh H h h hhhh MH nn mmnn n cbda 1 2 0101 1 2 2 1)( 2021-7-2229 例例 16-17 某冶煉廠工人的砷暴露年份與呼吸系統(tǒng)腫瘤死亡關(guān)系的資料某冶煉廠工人的砷暴露年份與呼吸系統(tǒng)腫瘤死亡關(guān)系的資料 砷暴露年份砷暴露年份 指標指標 0 1 1 5 5 1

23、5 15 合計合計 呼吸系統(tǒng)腫瘤的死亡數(shù)呼吸系統(tǒng)腫瘤的死亡數(shù) 51 17 13 34 115 觀察人年數(shù)觀察人年數(shù) 56105 7696 7223 11616 82640 線性趨勢線性趨勢檢驗（檢驗（trend 2 test） 0 H ：事件的發(fā)生率與劑量不相關(guān)，事件的發(fā)生率與劑量不相關(guān)， 1 H ：事件的發(fā)生率與劑量相關(guān)事件的發(fā)生率與劑量相關(guān) 暴露不是有或無，暴露不是有或無， 而是有幾個劑量而是有幾個劑量 2021-7-2230 暴露因素水平暴露因素水平 E0 E1 EK 取取 值值 X0 X1 XK 合計合計 出現(xiàn)結(jié)局事件的人數(shù)出現(xiàn)結(jié)局事件的人數(shù) d0 d1 dK D 觀察人年數(shù)觀察人年數(shù)

24、 T0 T1 TK T 線性趨勢檢驗（線性趨勢檢驗（trend 2 test） 0 H ：事件的發(fā)生率與劑量不相關(guān)，事件的發(fā)生率與劑量不相關(guān)， 1 H ：事件的發(fā)生率與劑量相關(guān)事件的發(fā)生率與劑量相關(guān) 22 200 22 00 (1)() ()() KK kk kk ki KK kkkk kk X T D TTX d T D TD TX TX T 2 )1( 本例本例， 50.30 )5699014113282640(82525115 )5699082640/115145(8263982640 2 22 2 P0.01，拒絕，拒絕 H0，可以認為砷暴露年限與呼吸系統(tǒng)腫瘤死亡率存在，可以認為砷暴露

25、年限與呼吸系統(tǒng)腫瘤死亡率存在 劑量劑量-反應(yīng)關(guān)系。反應(yīng)關(guān)系。 31 l觀察性研究觀察性研究(observational study)： 不對研究對象施加任何干預(yù)措施的情況下，通過不對研究對象施加任何干預(yù)措施的情況下，通過 觀察或訪問的方法，客觀地記錄被研究事物的狀況。觀察或訪問的方法，客觀地記錄被研究事物的狀況。 也稱非實驗研究（也稱非實驗研究（non- experimental study） 2. 最常見的觀察性研究：最常見的觀察性研究： 描述性：橫斷面研究描述性：橫斷面研究 分析性：病例分析性：病例-對照研究對照研究 隊列研究隊列研究 2021-7-2232 3. 橫斷面研究，常通過抽樣調(diào)

26、查，估計當前總體中橫斷面研究，常通過抽樣調(diào)查，估計當前總體中 某個量的平均值，或某個事件出現(xiàn)的概率；某個量的平均值，或某個事件出現(xiàn)的概率； l 常用抽樣方法包括：單純隨機抽樣、分層隨機抽樣常用抽樣方法包括：單純隨機抽樣、分層隨機抽樣 和多階段整群抽樣；和多階段整群抽樣； l 對于有限總體，估計總體參數(shù)時，必須考慮校正對于有限總體，估計總體參數(shù)時，必須考慮校正 因子因子 ； l 分層抽樣時，參數(shù)估計必須加權(quán)平均分層抽樣時，參數(shù)估計必須加權(quán)平均 1 N nN 2021-7-2233 4. 病例病例-對照研究有兩種設(shè)計：對照研究有兩種設(shè)計： 組與組病例組與組病例-對照對照 和個體匹配病例和個體匹配病

27、例-對照；對照； l 組與組病例組與組病例-對照，若暴露只分為有或無，不考慮對照，若暴露只分為有或無，不考慮 其他影響因素時，數(shù)據(jù)歸結(jié)為單個四個表；其他影響因素時，數(shù)據(jù)歸結(jié)為單個四個表； 可以計算：優(yōu)勢比的點估計、區(qū)間估計和假設(shè)檢驗可以計算：優(yōu)勢比的點估計、區(qū)間估計和假設(shè)檢驗 l 組與組病例組與組病例-對照，若暴露只分為有或無，按照某個對照，若暴露只分為有或無，按照某個 重要因素分層時，可歸結(jié)為多個四個表；可以計算：重要因素分層時，可歸結(jié)為多個四個表；可以計算： 公共優(yōu)勢比的點估計、區(qū)間估計和假設(shè)檢驗公共優(yōu)勢比的點估計、區(qū)間估計和假設(shè)檢驗 l 組與組病例組與組病例-對照，暴露分為對照，暴露分為C個水平時，可歸結(jié)為個水平時，可歸結(jié)為 一個一個2C 列聯(lián)表；可以分析暴露的劑量