醫(yī)學(xué)統(tǒng)計(jì)學(xué)13

1、醫(yī)學(xué)統(tǒng)計(jì)學(xué)（13）——復(fù) 習(xí),季聰華2012.12.26,Meta分析,什么是Meta分析？,是對具備特定條件的、同課題的諸多研究結(jié)果進(jìn)行定量合并的一類統(tǒng)計(jì)方法。優(yōu)點(diǎn)：Meta分析整合大量的研究報(bào)告，增加了樣本量，增加了結(jié)論的統(tǒng)計(jì)功效，解決分歧意見，比綜述的結(jié)論更客觀全面, 并可引出新見解。,Meta分析研究步驟：,遵循一個結(jié)構(gòu)式程序: – 提出研究問題 – 綜合搜索資料 – 嚴(yán)格評價(jià)資料 – 提取資料數(shù)據(jù)

2、– 分析并形成結(jié)果 （Meta分析方法）,步驟 1：提出問題,選題：是最重要最困難的事情需要考慮以下幾個方面： 重要性 爭議性 創(chuàng)新性 可行性,步驟 1：提出問題,P population/patients 人群/病人群 關(guān)心的是哪一類人群或病人I intervention/exposure 干預(yù)/暴露

3、 進(jìn)行什么治療、處于何種暴露C comparison or control 比較物或?qū)φ?對照物、處理策略、試驗(yàn)或暴露是什么？O outcome 結(jié)果 結(jié)局如何？,PICO原則,步驟 2：全面查找資料,1 制定檢索詞、檢索策略2 嚴(yán)格納入和排除標(biāo)準(zhǔn)3 全面檢索資料： Pubmed、Cochrane、EBMbase、OVID、Google、國內(nèi)全文數(shù)據(jù)庫、報(bào)告、會議、摘要

4、等所有可用資源4 文獻(xiàn)管理軟件：Note express、EndNote,步驟 3：資料評價(jià),對檢索資料進(jìn)行評價(jià)、評分，分5級。,Jadad評分量表1 隨機(jī)分組序列的產(chǎn)生方法2分：通過計(jì)算機(jī)產(chǎn)生的隨機(jī)序列或隨機(jī)數(shù)表產(chǎn)生的序列；1分：試驗(yàn)提到隨機(jī)分配，但產(chǎn)生隨機(jī)序列的方法未予交待；0分：半隨機(jī)或準(zhǔn)隨機(jī)試驗(yàn)，指采用交替分配病例的方法，如入院順序、出生日期單雙數(shù)；2 雙盲法2分：描述了實(shí)施雙盲的具體方法并且被認(rèn)為是恰當(dāng)?shù)?，如采?/p>

5、完全一致的安慰劑等；1分：試驗(yàn)僅提及采用雙盲法；0分：試驗(yàn)提及采用雙盲，但方法不恰當(dāng)，如比較片劑與注射劑而未提及使用雙偽法；3 退出與失防1分：對退出與失防的病例數(shù)和退出理由進(jìn)行了詳細(xì)的描述；0分：沒有提到退出與失防；,評價(jià)治療或篩查的證據(jù)質(zhì)量,預(yù)防、診斷、預(yù)后、治療和危害研究等領(lǐng)域的研究評價(jià),,,,,,步驟 4：抽取資料,,,,步驟 5：數(shù)據(jù)分析整理,計(jì)量資料：均數(shù)之差 例1 氟對女童 I I 掌骨皮質(zhì)厚度的影響,標(biāo)

6、準(zhǔn)化均數(shù)之差: di =( 1i - 2i) /Sci , i = 1, 2, 3 ,k權(quán)重系數(shù): wi = n 1i + n 2i,,步驟 6：結(jié)果展示,步驟 7：討論與總結(jié),評估該系統(tǒng)評價(jià)的 局限性 實(shí)用性 論證強(qiáng)度 對醫(yī)學(xué)研究的意義,Meta分析軟件：,RevMan 操作簡易 系統(tǒng)綜述和meta分析Stata 功能強(qiáng)大 發(fā)表偏倚定量及圖片編輯Met

7、a-dicsR SAS,Meta分析涉及的幾個統(tǒng)計(jì)處理,統(tǒng)計(jì)處理流程：a．異質(zhì)性檢驗(yàn)（齊性檢驗(yàn)）。 b．統(tǒng)計(jì)合并效應(yīng)量并進(jìn)行統(tǒng)計(jì)推斷。 c．采用“倒漏斗圖”了解潛在的發(fā)表偏倚。,異質(zhì)性檢驗(yàn)（齊性檢驗(yàn)）,異質(zhì)性檢驗(yàn)基本思想： 假如研究資料間的真實(shí)效應(yīng)量一致，那么實(shí)際效應(yīng)量間的差異可認(rèn)為是由誤差引起的。常用Q檢驗(yàn)方法檢驗(yàn)，P值50%,說明存在比較明顯的異質(zhì)性。,異質(zhì)性檢驗(yàn)（齊性檢驗(yàn)）,識別異質(zhì)性：觀察幾個獨(dú)立研究的差異

8、情況，是否合并？ Q檢驗(yàn) p值<0.1表示有異質(zhì)性,,亞組分析 Subgroup Analysis,根據(jù)性別、年齡、種族等進(jìn)行分組分析,,,,合并效應(yīng),RevMen軟件計(jì)算結(jié)果：,,,總體效應(yīng)檢驗(yàn)，p<0.05有意義,,發(fā)表偏倚 Publication Bias,陽性的研究結(jié)果發(fā)表的機(jī)會更多，過分強(qiáng)調(diào)陽性結(jié)果可能造成誤導(dǎo)。漏斗圖描述發(fā)表偏倚一項(xiàng)有爭議的研究得出來的是一個倒置漏斗圖：如果有發(fā)表偏倚，則圖形不對稱

9、：,,發(fā)表偏倚 Publication Bias,繪制漏斗圖分析 P<0.05 表示偏倚存在,,Stata 輸出：,Stata,Revman,R軟件,Egger's test------------------------------------------------------------------------------ Std_Eff | Coef.

10、Std. Err. t P>|t| [95% Conf. Interval]-------------+---------------------------------------------------------------- slope | .4270146 .286694 1.49 0.197 -.3099558 1

11、.163985 bias | .0547229 .947982 0.06 0.956 -2.382143 2.491589,如何提高,檢索,系統(tǒng)綜述/Meta分析,,英語,醫(yī)學(xué),統(tǒng)計(jì),,,,做一個系統(tǒng)綜述研究常需要多人協(xié)作,更多了解……,丁香園循證醫(yī)學(xué)專版：http://ebm.dxy.cn/bbs/post/page?bid=96&sty=0&age

12、=0,醫(yī)學(xué)統(tǒng)計(jì)學(xué),一、緒論二、統(tǒng)計(jì)描述三、概率分布四、參數(shù)估計(jì)與假設(shè)檢驗(yàn)五、t檢驗(yàn)六、方差分析七、雙變量相關(guān)與回歸分析八、x2檢驗(yàn),九、基于秩次的非參數(shù)檢驗(yàn)十、實(shí)驗(yàn)設(shè)計(jì)十一、調(diào)查設(shè)計(jì)十二、協(xié)方差分析十三、多重線性回歸分析十四、logistic回歸分析十五、生存分析,,,《中華醫(yī)學(xué)雜志》對來稿統(tǒng)計(jì)學(xué)處理的有關(guān)要求,醫(yī)學(xué)統(tǒng)計(jì)學(xué),一、緒論二、統(tǒng)計(jì)描述三、概率分布四、參數(shù)估計(jì)與假設(shè)檢驗(yàn)五、t檢驗(yàn)六、方差分析

13、七、雙變量相關(guān)與回歸分析八、x2檢驗(yàn),九、基于秩次的非參數(shù)檢驗(yàn)十、實(shí)驗(yàn)設(shè)計(jì)十一、調(diào)查設(shè)計(jì)十二、協(xié)方差分析十三、多重線性回歸分析十四、logistic回歸分析十五、生存分析,統(tǒng)計(jì)工作的步驟及內(nèi)容,統(tǒng)計(jì)描述統(tǒng)計(jì)推斷,參數(shù)估計(jì)假設(shè)檢驗(yàn),統(tǒng)計(jì)指標(biāo)統(tǒng)計(jì)圖表,,,,研究設(shè)計(jì)收集資料整理資料分析資料,,統(tǒng)計(jì)描述的常用指標(biāo),例數(shù)（n）均數(shù)±標(biāo)準(zhǔn)差（x±s）最小值，最大值（min，max）中位數(shù)（M，

14、median）四分位間距、百分位數(shù)頻數(shù)、百分比,,,,,,,,min,median,p75,max,p25,t檢驗(yàn),資料由什么設(shè)計(jì)而來：完全隨機(jī)設(shè)計(jì)、配對設(shè)計(jì)統(tǒng)計(jì)學(xué)描述（例數(shù)、均數(shù)、標(biāo)準(zhǔn)差）統(tǒng)計(jì)學(xué)檢驗(yàn)方法（單組、兩組、配對）統(tǒng)計(jì)學(xué)檢驗(yàn)結(jié)果（根據(jù)方差是否齊，選擇不同結(jié)果），t值，p值專業(yè)意義,【例1】為了驗(yàn)證腎上腺素有無降低呼吸道阻力的作用，用豚鼠12只進(jìn)行支氣管灌流法實(shí)驗(yàn)，在注入定量腎上腺素前后，測得豚鼠支氣管功能反映在儀器

15、上的灌流速度(每分鐘灌流滴數(shù))，結(jié)果見下表 。問：腎上腺素能否降低豚鼠呼吸道阻力?,,結(jié)果解讀,結(jié)果解讀1,N：用藥前、用藥后樣本例數(shù)=12均值：用藥前、用藥后均值分別為44.83、52.83標(biāo)準(zhǔn)差：用藥前、用藥后標(biāo)準(zhǔn)差分別為9.815、6.952,結(jié)果解讀2,t：統(tǒng)計(jì)量t=-2.653Sig（雙側(cè)）：p值=0.022均值差值：兩個均數(shù)的差值=-8.000差值的95%CI：-14.636~-1.364,【例2】為觀察中成藥青

16、黛明礬片對急性黃疸肝炎的退黃效果，以單用輸液保肝的患者作為對照進(jìn)行了完全隨機(jī)設(shè)計(jì)觀察，受試對象為黃疸指數(shù)在 30~50 之間的成年患者，觀測結(jié)果為退黃天數(shù)，數(shù)據(jù)見下表。試比較中藥組與對照組退黃天數(shù)有無差別?,結(jié)果解讀,結(jié)果解讀1,N：組1、組2的樣本例數(shù)分別為8、9；均值：組1、組2的均值分別為11.75、21.44；標(biāo)準(zhǔn)差：組1、組2的標(biāo)準(zhǔn)差分別為5.392、6.405；,結(jié)果解讀2,方差齊性檢驗(yàn)F=0.002，p=0.968

17、，p>0.05方差相等，t檢驗(yàn)結(jié)果選第一行。,,,結(jié)果解讀3,t：統(tǒng)計(jì)量t=-3.351Sig（雙側(cè)）：p值=0.004，p<0.05均值差值：兩個均數(shù)的差值=-9.694差值的95%CI：-15.861~-3.528,,,方差分析,資料由什么設(shè)計(jì)而來統(tǒng)計(jì)學(xué)描述（例數(shù)、均數(shù)、標(biāo)準(zhǔn)差）統(tǒng)計(jì)學(xué)檢驗(yàn)方法,兩兩比較統(tǒng)計(jì)學(xué)檢驗(yàn)結(jié)果，F(xiàn)值，p值專業(yè)意義,【例3】研究單味中藥對小白鼠細(xì)胞免疫機(jī)能的影響，把 40 只小白鼠隨

18、機(jī)分成四組，每組 10 只，雌雄各半，用藥15天后，進(jìn)行E -玫瑰花結(jié)形成率( E-SFC )測定，結(jié)果見下表，試比較三種單味中藥的免疫作用是否一樣?,結(jié)果解讀1,描述了4組分別的和總的均值、標(biāo)準(zhǔn)差、標(biāo)準(zhǔn)誤、95%CI、最小、最大值。,結(jié)果解讀2,進(jìn)行了方差齊性檢驗(yàn)，p=0.048，p<0.05，表示方差不齊。,結(jié)果解讀3,總的組間差異比較，F(xiàn)=83.546，p=0.000，p<0.01，表明組間差異有統(tǒng)計(jì)學(xué)意義。,結(jié)果解讀

19、4：組間均值的兩兩比較,,,結(jié)果解讀5,S-N-K檢驗(yàn)結(jié)果將無統(tǒng)計(jì)學(xué)意義的比較組列在同一列中。即樣本均數(shù)顯示在同一列時(shí)，表示兩組總體均數(shù)差別無統(tǒng)計(jì)學(xué)意義。,【例4】按2X2析因設(shè)計(jì)進(jìn)行臨床研究，把 20名胃癌患者隨機(jī)分到各實(shí)驗(yàn)組，治療 3 個月后，應(yīng)用胃癌患者QOL的測評量表評價(jià)方法，以總分評定A、B兩種中藥配合抗癌藥物治療胃癌患者的生存質(zhì)量，結(jié)果如下表。試分析A、B 兩藥的治療效果，以及 A、B 兩藥間是否存在交互作用?,結(jié)果解讀1,

20、A藥有抗胃癌作用，B藥有抗胃癌作用，兩藥同時(shí)用抗胃癌作用更好（協(xié)同作用）。,第七步：完成，解讀結(jié)果,,,,簡單線性回歸,直線趨勢：散點(diǎn)圖回歸方程是否有意義：方差分析結(jié)果回歸系數(shù)的意義寫出回歸方程相關(guān)系數(shù)及其意義,【例5】11名糖尿病患者的血清總膽固醇含量( minol/L)與空腹血糖(mmol/L)的測量值如下，試進(jìn)行直線回歸分析。,第3步：作散點(diǎn)圖（4）,查看散點(diǎn)圖結(jié)果：有直線趨勢，異常點(diǎn)不明顯。,第4步：線性回歸（5）,結(jié)果

21、解讀1：r=0.852，p=0.001決定系數(shù)R2=0.726，表示Y的變異有72.6%與回歸有關(guān)。,第4步：線性回歸（5）,結(jié)果解讀2：直線回歸的方差分析F=23.789，p=0.001,第4步：線性回歸（5）,結(jié)果解讀3：直線回歸方程的截距a=3.433，回歸系數(shù)b=1.620；T檢驗(yàn)t=4.877，p=0.001。直線回歸方程可以寫成血糖=3.433+1.620*總膽固醇含量,x2檢驗(yàn),應(yīng)用條件：不同性質(zhì)的資料，不同的方法

22、統(tǒng)計(jì)學(xué)描述：例數(shù)、百分比x2值、p值的選擇專業(yè)意義,【例6】某醫(yī)院把慢性支氣管炎患者376名，隨機(jī)分為2組，分別用中西醫(yī)結(jié)合法和西醫(yī)法治療，結(jié)果見表。問兩種療法治療慢性支氣管炎病人的治愈率是否有差別？,結(jié)果解讀（1）,結(jié)果解讀：中西醫(yī)組的治愈率為98.2%，西醫(yī)組的治愈率為74.0%。,,,結(jié)果解讀（2）,結(jié)果解讀：x2=56.772，p=0.000兩種療法治療慢性支氣管炎病人的治愈率的差別有統(tǒng)計(jì)學(xué)意義,四格表x2檢驗(yàn)結(jié)果的選擇,

23、1）當(dāng)n≥40，所有理論值≥ 5時(shí)，用pearson卡方檢驗(yàn)；或者確切概率法。2）當(dāng)n≥40，但有理論頻數(shù)1≤理論值＜5時(shí)，用連續(xù)校正的卡方檢驗(yàn)；或者確切概率法。3） n<40或有理論值<1，或P≈α?xí)r，用確切概率法。,非參數(shù)檢驗(yàn),應(yīng)用條件統(tǒng)計(jì)量和p值的選擇,【例7】20名正常人和32名鉛作業(yè)工人尿鉛定性檢查結(jié)果如表。問鉛作業(yè)工人尿鉛是否高于正常人?,結(jié)果解讀,結(jié)果解讀：Z=4.503,P=0.000,,,常用的實(shí)驗(yàn)設(shè)

24、計(jì)類型,完全隨機(jī)設(shè)計(jì) 配對設(shè)計(jì) 隨機(jī)區(qū)組設(shè)計(jì)拉丁方設(shè)計(jì)交叉設(shè)計(jì) 析因設(shè)計(jì)　 正交設(shè)計(jì)　重復(fù)測量設(shè)計(jì),研究時(shí)有無設(shè)計(jì)干預(yù)因素？,實(shí)驗(yàn)性研究,觀察性研究,是否隨機(jī),有無對照組,分析性研究,描述性研究,時(shí)間方向,,,,,,,,,,,,隊(duì)列研究,非隨機(jī)對照試驗(yàn),,,隨機(jī)對照試驗(yàn),病例對照研究,橫斷面研究,,,,有,無,有,無,否,是,暴露→結(jié)局,結(jié)局→暴露,,處理因素、受試對象、實(shí)驗(yàn)(試驗(yàn))效應(yīng)是實(shí)驗(yàn)設(shè)計(jì)的三個基本

25、要素。能說出具體案例中的三個基本要素,多元線性回歸,回歸方程的建立：逐步回歸方法回歸方程的檢驗(yàn)：方差分析偏回歸系數(shù)的意義寫出回歸方程,【例8】觀察30名兒童血中血紅蛋白、鈣、鎂、鐵的含量。試研究血紅蛋白與微量元素的關(guān)系？,結(jié)果解讀（1）,結(jié)果解讀：逐步回歸指標(biāo)進(jìn)入的過程,結(jié)果解讀（2）,結(jié)果解讀：復(fù)相關(guān)系數(shù)和決定系數(shù),結(jié)果解讀（3）,結(jié)果解讀：回歸模型的方差分析,結(jié)果解讀（4）,結(jié)果解讀：回歸方程：血紅蛋白=1.129+0.0

26、31*血鐵-0.043*血鈣,,生存分析,資料特點(diǎn)生存函數(shù)曲線統(tǒng)計(jì)描述：中位生存時(shí)間多因素分析：COX回歸RR值及其95%CI專業(yè)意義,【例9】50例某惡性腫瘤的生存時(shí)間及有關(guān)預(yù)后因素資料如表。X1為患者年齡、X2為病理類型（0=低分化，1=高分化），X3為有無淋巴結(jié)浸潤（0=無，1=有），X4為緩解出院后的鞏固治療（0=無，1=有），T為患者的生存時(shí)間，S為患者的隨訪結(jié)局（0=截尾，1=死亡）。試進(jìn)行生存分析。,結(jié)果解讀（1

27、）,結(jié)果解讀：有意義的因素，相對危險(xiǎn)度RR及其95%CI。,,,X3為有無淋巴結(jié)浸潤（0=無，1=有），RR>1,危險(xiǎn)因素X4為緩解出院后的鞏固治療（0=無，1=有），RR<1,保護(hù)因素,結(jié)果解讀（2）,結(jié)果解讀（3）,2011年試卷,配對設(shè)計(jì)的隨機(jī)分組方法方法：先將 60個受試對象按某種或某些特征排列并依次編號1-60，特征相近的 2個受試對象配成1對，這樣分成30對，再將每對的2例隨機(jī)分配到 A、B兩種處理中。,配對

28、t檢驗(yàn)（-5.37104~-2.76229）,有序分類資料基于秩次的非參數(shù)檢驗(yàn),P=0.612，兩組的不良反應(yīng)發(fā)生率差異無統(tǒng)計(jì)學(xué)意義（p>0.05）,方差不齊，采用隨機(jī)區(qū)組設(shè)計(jì)的非參數(shù)檢驗(yàn),S-N-K法皮質(zhì)醇3優(yōu)于2、4，優(yōu)于1,3、4低于12與1差異沒有統(tǒng)計(jì)學(xué)意義,析因設(shè)計(jì)60,液氮-196，保存5天活力最大,時(shí)間*溫度時(shí)間與溫度有交互作用,,應(yīng)變量Y為連續(xù)變量自變量X可為連續(xù)、有序分類或無序分類變量,1.線性趨勢