1、第十二章 重復(fù)測量設(shè)計(jì)資料的方差分析Repeated measures ANOVA,暨南大學(xué)醫(yī)學(xué)院醫(yī)學(xué)統(tǒng)計(jì)學(xué)教研室 林漢生,重復(fù)測量的概念,本章的重復(fù)測量實(shí)驗(yàn)： 對同一觀察對象（病人，受試者、動物、標(biāo)本，等）的同一觀察指標(biāo)在不同時間點(diǎn)上進(jìn)行的多次測量。 同一受試者的多次測量之間可能存在某種相關(guān)性，用通常的統(tǒng)計(jì)方法不能揭示內(nèi)在的特點(diǎn)，有時甚至?xí)贸鲥e誤的結(jié)論。,第一節(jié) 重復(fù)測量資料的數(shù)據(jù)特征,前后測量設(shè)計(jì) 設(shè)立對照的前后測量設(shè)計(jì) 重復(fù)測量設(shè)計(jì) 設(shè)立對照的重復(fù)測量設(shè)計(jì),一、前后測量設(shè)計(jì),前后測量設(shè)計(jì)的統(tǒng)計(jì)方法,前提條件：在假定測量時間對觀測結(jié)果沒有影響時，才能推斷處理是否有效。 統(tǒng)計(jì)方法：用配對

2、t 檢驗(yàn)或配對秩和檢驗(yàn) 高血壓患者治療前后的舒張壓平均下降了16mmHg。 經(jīng)配對t檢驗(yàn)，t =16.18， P0.000 能說明治療有效嗎？ 住院休息，環(huán)境和情緒的改變？考慮了嗎？,二、設(shè)立對照的前后測量設(shè)計(jì),高血壓患者治療前后的舒張壓（mmHg）,確定療效的前后測量設(shè)計(jì)必須增加平行對照,平行對照：處理組和對照組的受試對象來自同一人群，試驗(yàn)同時進(jìn)行。 對兩組治療前后的血壓差值進(jìn)行兩樣本t檢驗(yàn)，信息損失較大。如編號1（130-11416），編號8（138-12216）。,三、重復(fù)測量設(shè)計(jì),當(dāng)前后測量設(shè)計(jì)的重復(fù)測量次數(shù)m3時，稱重復(fù)測量設(shè)計(jì)或重復(fù)測量數(shù)據(jù)。 重復(fù)測量設(shè)計(jì) 設(shè)立對照的重復(fù)測量設(shè)計(jì)

3、未設(shè)立對照的重復(fù)測量設(shè)計(jì),1. 設(shè)立對照的重復(fù)測量設(shè)計(jì),將手術(shù)要求基本相同的15名患者隨機(jī)分3組，在手術(shù)過程中分別采用A，B，C三種麻醉誘導(dǎo)方法，在T0（誘導(dǎo)前）、T1、T2、T3，T4 五個時相測量患者的收縮壓，數(shù)據(jù)記錄見表。,2. 未設(shè)立對照的重復(fù)測量數(shù)據(jù),越是相鄰的點(diǎn)，相關(guān)系數(shù)越大還是越??？,越是相鄰的時間點(diǎn)，數(shù)據(jù)之間的相關(guān)性越大,各列數(shù)據(jù)處理之間存在不同程度的相關(guān)性，而且越是相鄰的時間點(diǎn)，數(shù)據(jù)之間的相關(guān)性越大。 理論上講，當(dāng)重復(fù)測量數(shù)據(jù)之間相關(guān)，用隨機(jī)區(qū)組設(shè)計(jì)方差分析考察各時間點(diǎn)之間的差別，會增加I類錯誤的概率。 如果只重復(fù)2次，不涉及多重比較，故不涉及該問題。,隨機(jī)區(qū)組設(shè)計(jì)數(shù)據(jù),先將

4、全部受試對象按某種或某些特征分為若干個區(qū)組（block），使每個區(qū)組內(nèi)的觀察對象隨機(jī)地接受研究因素某一水平的處理。,重復(fù)測量數(shù)據(jù)與隨機(jī)區(qū)組設(shè)計(jì)的區(qū)別,第二節(jié) 重復(fù)測量數(shù)據(jù)的兩因素兩水平分析,高血壓患者治療前后的舒張壓（mmHg）,可以分析：,時間主效應(yīng)：治療前、后的舒張壓均數(shù)是否不同？ 組間主效應(yīng)：處理組與對照組的舒張壓均數(shù)是否不同？ 交互作用：組間的舒張壓均數(shù)差異是否隨時間改變（或各組舒張壓的下降幅度是否相同）？,AnalyzeGeneral Linear ModelRepeated Measures,將“治療前”,“治療后”選入將“分組”選入,分別 單擊Plots 、Options,將分組

5、和time選入，單擊Add，單擊Continue,選擇“Descriptive statistics”單擊Continue,單擊OK,結(jié)果解釋,SPSS結(jié)果解釋,這是多變量方差分析結(jié)果，不是本章討論內(nèi)容,“球?qū)ΨQ”檢驗(yàn)結(jié)果：因只重復(fù)兩次測量，不能進(jìn)行“球?qū)ΨQ”檢驗(yàn)。,“球?qū)ΨQ”假設(shè),球?qū)ΨQ假設(shè)：各時間點(diǎn)數(shù)據(jù)的兩兩相關(guān)系數(shù)為常數(shù)。 兩個時間點(diǎn)的重復(fù)測量，不進(jìn)行“球?qū)ΨQ”檢驗(yàn)。,治療前后舒張壓均數(shù)有差別( F =55.01, P=0.000)； 組間的舒張壓均數(shù)差異隨時間改變( F =18.78, P =0.000)，處理組舒張壓均數(shù)的下降幅度大于對照組。,處理組與對照組舒張壓均數(shù)的差別無統(tǒng)計(jì)學(xué)意

6、義（F = 1.57, P = 0.226）。,第三節(jié) 重復(fù)測量數(shù)據(jù)的兩因素多水平分析,例12-3：將手術(shù)要求基本相同的15名患者隨機(jī)分3組，在手術(shù)過程中分別采用A，B，C三種麻醉誘導(dǎo)方法，在T0（誘導(dǎo)前）、T1、T2、T3，T4 五個時相測量患者的收縮壓，數(shù)據(jù)記錄見表。試進(jìn)行方差分析。,數(shù)據(jù)文件,Analyze General Linear Model Repeated Measures,結(jié)果解釋,“球?qū)ΨQ”假設(shè)成立（ P=0.178 ），可選擇Sphericity Assumed法,不同麻醉誘導(dǎo)時相間的收縮壓均數(shù)存在差別（F = 106.56, P = 0.000）； 不同誘導(dǎo)方法患者收縮

7、壓均數(shù)的差別隨時相改變（F = 19.10, P = 0.000），其中A組不同誘導(dǎo)時相的收縮壓較為穩(wěn)定(見圖)。,A組不同誘導(dǎo)時相的收縮壓較為穩(wěn)定,統(tǒng)計(jì)圖的文字說明最好用英文，不僅可以在SPSS中進(jìn)行各種編輯，而且著名的醫(yī)學(xué)期刊一般都要求統(tǒng)計(jì)圖用全英文說明。,不同麻醉誘導(dǎo)方法患者收縮壓均數(shù)的差別有統(tǒng)計(jì)學(xué)意義（F = 5.78, P =0.017）。,每個患者的收縮壓變化：單擊Plots,單擊Add，返回,（1）哪一組的收縮壓變化較為穩(wěn)定？ （2）為什么要比較個體的變化情況,第四節(jié) 重復(fù)測量設(shè)計(jì)（無對照組）,與隨機(jī)區(qū)組設(shè)計(jì)的方差分析比較,數(shù)據(jù)文件,Analyze General Linear

8、Model Repeated Measures,沒有分組變量！單擊OK,“球?qū)ΨQ” 假設(shè)不成立（ P=0.008 ），可選擇Epsilon ()對自由度進(jìn)行校正?？蛇x擇0.528，或0.657，或0.333。折中的方法是選擇Greenhouse-Geisser 0.528。,“球?qū)ΨQ”假設(shè),球?qū)ΨQ假設(shè)：各時間點(diǎn)數(shù)據(jù)的兩兩相關(guān)系數(shù)為常數(shù)。 如不滿足球?qū)ΨQ假設(shè)，應(yīng)采用校正方法。方法是校正后的自由度變小，F(xiàn)界值變大。如未校正前F0.05,(3.21)=3.07，當(dāng)校正系數(shù)為0.536時，校正后F界值為F0.05,(1.6,11.3) F0.05,(1,11 ) 4.84。從而減小I型錯誤。,F=79.14 盡管兩者的P值在小數(shù)點(diǎn)3位時一樣：P0.000，但精確的P值不一樣。 Sphericity Assumed: P=1.303629031394 e-011 Greenhouse-Geisser: P=5.36754297376 e-007 e-007 的含義是 10-7,第五節(jié) 重復(fù)測量數(shù)據(jù)統(tǒng)計(jì)分析常見的誤用情況,重復(fù)進(jìn)行各時點(diǎn)的 t 檢驗(yàn) 差值缺乏效度 忽略個體曲線變化特征,作業(yè)：P199-200,簡答題：1，3 計(jì)算分析題：2，5,謝謝大家！,感謝您的觀看！,工作