重復(fù)測量的多因素方差分析

1、SPSS學(xué)習(xí)筆記之一一重復(fù)測量的多因素方差分析轉(zhuǎn)1、概述重復(fù)測量數(shù)據(jù)的方差分析是對同一因變量進(jìn)行重復(fù)測量的一種試驗設(shè)計技 術(shù)。在給予一種或多種處理后，分別在不同的時間點上通過重復(fù)測量同一個受試 對象獲得的指標(biāo)的觀察值，或者是通過重復(fù)測量同一個個體的不同部位 （或組織） 獲得的指標(biāo)的觀察值。重復(fù)測量數(shù)據(jù)在科學(xué)研究中十分常見。分析前要對重復(fù)測量數(shù)據(jù)之間是否存在相關(guān)性進(jìn)行球形檢驗。如果該檢驗結(jié)果為P> 0.05,則說明重復(fù)測量數(shù)據(jù)之間不存在相關(guān)性，測量數(shù)據(jù)符合 Huynh-Feldt條件，可以用單因素方差分析的方法來處理；如果檢驗結(jié)果 P< 0.05 ,則說明重復(fù)測量數(shù)據(jù)之間是存在相關(guān)性

2、的，所以不能用單因素方差分析的 方法處理數(shù)據(jù)。在科研實際中的重復(fù)測量設(shè)計資料后者較多，應(yīng)該使用重復(fù)測量設(shè)計的方差分析模型。球形條件不滿足時常有兩種方法可供選擇：（1）采用MANOVA多變量方差 分析方法）；（2）對重復(fù)測量ANOV檢驗結(jié)果中與時間有關(guān)的F值的自由度進(jìn) 行調(diào)整。2、問題新生兒胎糞吸入綜合征（MAS是由于胎兒在子宮內(nèi)或著生產(chǎn)時吸入了混有 胎糞的羊水，從而導(dǎo)致呼吸道和肺泡發(fā)生機械性阻塞，并伴有肺泡表面活性物質(zhì)失活，而且肺組織也會發(fā)生化學(xué)性炎癥， 胎兒出生后出現(xiàn)的以呼吸窘迫為主， 同 時伴有其他臟器受損現(xiàn)象的一組綜合征11。血管內(nèi)皮生長因子（vascular endothelial g

3、rowth factor , VEGF是一種有絲分裂原，它特異作用于血管內(nèi)皮 細(xì)胞時，能夠調(diào)節(jié)血管內(nèi)皮細(xì)胞的增殖和遷移，從而使血管通透性增加。而本實驗旨在通過觀察分析給予外源性肺表面活性物質(zhì)治療前后胎糞吸入綜合征患兒 血清中VEGF勺含量變化，評價藥物治療的效果。將收治的診斷胎糞吸入綜合癥的新生兒共 42名。將患兒隨機分為肺表面活 性物質(zhì)治療組（PS組）和常規(guī)治療組（對照組），每組各 21例。PS組和對照組 兩組所有患兒均給予除用藥外的其他相應(yīng)的對癥治療。PS組患兒給予牛肺表面活性劑PS 70mg/kg治療。采集PS組及對照組患兒0小時，治療后24小時和72 小時靜脈血2ml,離心并提取上清液

4、后保存?zhèn)溆貌⒂涗浹逯蠽EGFF勺含量變化情況。結(jié)果如下：表3.：血清中tEGF的原始數(shù)據(jù)空翦 組I 03 1 77200W0,5SW1 3141.5160-987L030,S3d0.7460.55d 1.12 1 IdI 211 0 9S3 0 P34 0 82S 1 164 0.5P4 1 0: 0 SS424，J 的1 703 1 548 I 23 I 135 1 422 1 12? 1 345 1 的8 1 235 1 33"3m X31.2141.2341.313 1.444 1.211.3541 24g 1.B4 0,7672 WL.703 L54S 1,23 1.435

5、 1.422 L823 1.345 1.41: 1.43? 1336 LW1343I 3541334 1 413 I 464 1.271 1356 L298 13341.423治 療 組山陽O SS 1 051.2261.5580,854 L03 0 778 1 03 0.8560 7S0.S68 1.030 64 L03 0J2 1 272 1 00C1 OS 0.8621.05 0 8S424小時I 203 1 068 1 23 1335 1322 1423 0 921 0J87 1 23: L236 I 432 1 143L2140J341.213 1.1641.21 0.9561.29

6、S 1.154 0.76二1時1.203 1.368 1.423 !335 I 322 I IB I 463 I 3S7 t 235 I 256 L452 1 1431.2141.3341 213 1.3641.210.9S6 UM 1.1340J ?63、統(tǒng)計分析建立數(shù)據(jù)文件變量視圖：名稱類型寬度T標(biāo)簽值二|gmup數(shù)值(N)80分組1,對照組ti數(shù)值(N)83。小時無t2數(shù)值(N)1324小時無13數(shù)值(N)8372小時無數(shù)據(jù)視圖:groupUt2t3i11.0301.70317031I-1.7721.5481.548I1.0941.2301.230_r1.5961.135143511.3

7、141.4221.42211.5161.1231.8231.9871.3451.345L 11.0初1.0981.4121.8361.2351.435菜單選擇:轉(zhuǎn)換（D 分析的）直銷網(wǎng)）圖昨售）并用程序畫口也0 耐助報告描述線計表（D比較均值圖）一眼段性模型0）義俄性模型溫含銀型Q9相關(guān)回歸首先進(jìn)入如下對話框，在“被試內(nèi)因子名稱”中輸入“ time”，“級別數(shù)” 輸入3,因為每個患者重復(fù)測量了 3次。后點擊“添加”按鈕。此時下方“定義”按鈕變?yōu)榭捎?，點擊進(jìn)入下列對話 框：目重復(fù)度量區(qū)）步分蛆jgroup 6 0小時tif- 9 2到、時【 / 72小時R3群體內(nèi)I部變量）模型睡一 對比例）繪制

8、兩兩比軾也）一保存- - - -1選項血m二 -重置晅）取消 幫助將“group”選入“因子列表”框，t1-t3分別選入“全體內(nèi)變量（time） 框內(nèi)，如下圖所示：點擊右上角“模型”按鈕，進(jìn)入以下對話框，選擇“設(shè)定”，將“ time”選 入“全體內(nèi)模型”框，“ group”選入“群體間模型”框，“構(gòu)建項”選擇“主 效應(yīng)”。下方的平方和選“類型III ”，這是對于平衡數(shù)據(jù)。如果兩組樣本量不 等，則選擇“類型IV”。點擊“繼續(xù)”返回，點擊“繪制”按鈕。進(jìn)入下面對話框：將“ time”選 入“水平軸”，group選入“單圖”，然后點擊“添加”按鈕，下面框中會顯示 “time*group ”。點擊“繼

9、續(xù)”返回，點擊“兩兩比較”按鈕，將 group選入右側(cè)“兩兩比 較檢驗”框中，選中復(fù)選框“ LSU 。點擊“繼續(xù)”返回，點擊“選項”按鈕，進(jìn)入下面對話框：將 time選入 “顯示均值框”，選中“比較主效應(yīng)”復(fù)選框，選中下方“描述統(tǒng)計”復(fù)選框。下方顯著性水平設(shè)為0.05。點擊“繼續(xù)”返回，點擊“確定”輸出結(jié)果4、結(jié)果解讀：描述性統(tǒng)計量分組均值標(biāo)準(zhǔn)偏差N口小時對照狷.07107,35672021治療組.97924,20085721忠計,97545.235950422 4小時對照組L28495,16373321治療蛆1.15876,14153721總計1.22136J 64096427 2小時對照蛆

10、1.41171,13744121治療俎1：25567,1 3660121總計1.3336SJ 5758342這是一個關(guān)于各個時間點的兩組數(shù)據(jù)描述性統(tǒng)計。M刑Kh®的球米度檢物t度量 MEASURE, 1主件內(nèi)力民Mmk附的收遁但卡方dfStQEptiton1Or4*nhou$e- OttSMfHuynh-Feiat下用9014 48920017B3B07500trnv設(shè),即標(biāo)看正女w-團費量的混襄D方若矩際與1-1寺位矩蟀成也倒.可用于®麒著性平均檢在W體內(nèi)/國檢中顯示"BMlb6 iftif 婕 + group 主住內(nèi)設(shè)時t>me這是球形檢驗結(jié)果，p=0.

11、001<0.05 ,所以不滿足球形分布假設(shè)，需要進(jìn)行多 變量方差分析或者自由度調(diào)整，SPSSg下來會給出以上兩種結(jié)果。名變量檢驗b值F哪也誤差dfSigtimePHI即的蜂蹤,67941 297*200039000000Wii心的 Lambda32141 297200039000000Hoteling的蹋豚2 11841.207*200031000000Roy的最大相211641 29732.00039000,000time 1 groupFill川的說雅,0911.946J200039000,156WHKs的 Lambda9091 340J200039D00156Hoteling的跟蹤

12、,100194y2.00039000156Roy的最大相,1001 940*2000I39000J56精船統(tǒng)計量b也計 翻,距+ grojp 主住內(nèi)設(shè)十time這是進(jìn)行多變量方差分析的結(jié)果，給出了 4種統(tǒng)計量，它們的檢驗結(jié)果一致， time的P<0.001,說明各個時間點的數(shù)據(jù)的差異有統(tǒng)計學(xué)意義， time*group的 P>0.05,說明時間和分組無交互作用，說明時間因素（即 0小時、24小時、72 小時）的作用不隨分組（即治療組和對照組）的不同而不同。主體內(nèi)效應(yīng)的程雅度量MEASURE戟山也平方和df均方F$igttme采用的障花度2 52221川142720 1000Gree

13、nhouse-Gelsser26221 5261 84942.720000HuynfvFeldt2S?21 6141.74842.72000029221 ooo282242720,000time * group家用的尊索度,1592,0802 414096G-nhousa&elsser.1591 5261042.414.111HuynhvFeldt1591,614,ogg2 414TTO.1591 000,1392414.128儡羞（time）采用的殍是度264280 1.0332 6"61 064,043HuyntyFelcft2 64 264 571,041284240,

14、000,066所謂“主體內(nèi)”，即是重復(fù)測量的各個時間點。上表是用各個時間點進(jìn)行分 組的方差分析表，給出4種統(tǒng)計量，第一種為滿足球星假設(shè)的情況，后三種對自 由度進(jìn)行了校正，本題目中不滿足球形分布假設(shè)，只能看下面的三種檢驗方法。 結(jié)果解釋同上一個表。這個圖可以直觀地看出測量指標(biāo)隨時間的變化趨勢。治療組與對照組兩組資料隨時間變化的趨勢大致相同，治療組血清中VEGF勺含量較對照組呈下降趨我們還可以給出在每個時間點上兩個分組之間的比較，需要用到多變量方差分析：操作步驟如下：跟之前操作類似，不贅述，看圖就行。I結(jié)果輸出主體間效屈的檢驗源因變量III盟平方和df均方F3ig校正模型必、時OOP1001007933"讖167b11677.13901172小時25B6125613,414,001K距制通39963139 g63476 910OQQ2 4小時627031E270326772E800072>W74,707174 7073919.311,000groupCkNt00110010079332討網(wǎng)1671.1677.13901172-M,2561,25613,414oot(kbB3.35240,06424/m,93740,0237 24時76240,019總計孫阿43,316 14224W63,