spss方差分析

1、spssspss 方差分方差分析析（ANOVAANOVA）河北醫(yī)科大學(xué)公共衛(wèi)生學(xué)院河北醫(yī)科大學(xué)公共衛(wèi)生學(xué)院衛(wèi)生統(tǒng)計(jì)學(xué)教研組衛(wèi)生統(tǒng)計(jì)學(xué)教研組2022年7月4日星期一18時(shí)08分45秒例例5-1 5-1 在大腸濕熱征模型研究中，將在大腸濕熱征模型研究中，將3232只大鼠隨機(jī)只大鼠隨機(jī)分為分為4 4組，每組組，每組8 8只，分別為正常對(duì)照組只，分別為正常對(duì)照組( (甲組甲組) )、大、大腸桿菌模型組腸桿菌模型組( (乙組乙組) )、輪狀病毒模型組、輪狀病毒模型組( (丙組丙組) )和大和大腸桿菌加輪狀病毒混合模型組腸桿菌加輪狀病毒混合模型組( (丁組丁組) )。測(cè)得大鼠血。測(cè)得大鼠血清白細(xì)胞介素清白

2、細(xì)胞介素-2(IL-2)-2(IL-2)水平見(jiàn)表水平見(jiàn)表5-15-1，試比較不同，試比較不同大鼠模型的血清大鼠模型的血清IL-2IL-2水平是否有差別。水平是否有差別。 （一）完全隨機(jī)設(shè)計(jì)的ANOVA 2022年7月4日星期一18時(shí)08分45秒表5-1 不同大鼠模型血清IL-2水平(ng/ml)的比較 2022年7月4日星期一18時(shí)08分46秒方法一方法一: : IL_2IL_2: :血清白細(xì)胞介素血清白細(xì)胞介素-2-2水平；水平； group:group:分組，分組，1=1=甲組，甲組，2=2=乙組，乙組，3=3=丙組，丙組，4=4=丁組。丁組。 均數(shù)比較均數(shù)比較單因素方差分析單因素方差分析

3、2022年7月4日星期一18時(shí)08分46秒因變量因變量?jī)蓛杀容^兩兩比較選項(xiàng)（方差齊性檢驗(yàn)）選項(xiàng)（方差齊性檢驗(yàn)）因素因素2022年7月4日星期一18時(shí)08分46秒假定方差齊假定方差齊假定方假定方差不齊差不齊2022年7月4日星期一18時(shí)08分46秒方差齊性檢驗(yàn)方差齊性檢驗(yàn)2022年7月4日星期一18時(shí)08分46秒方差齊性檢驗(yàn)方差齊性檢驗(yàn)P P值值F F值值P P值值SSSSMSMS組間組間組內(nèi)組內(nèi)總總F F值值2022年7月4日星期一18時(shí)08分46秒P P值值多重比較多重比較2022年7月4日星期一18時(shí)08分46秒P P值值2022年7月4日星期一18時(shí)08分46秒SNKSNK相同列，無(wú)差別

4、相同列，無(wú)差別不同列，有差別不同列，有差別2022年7月4日星期一18時(shí)08分46秒方法二方法二廣義線(xiàn)性模型廣義線(xiàn)性模型單因變量單因變量2022年7月4日星期一18時(shí)08分46秒兩兩比較兩兩比較選項(xiàng)選項(xiàng)2022年7月4日星期一18時(shí)08分46秒假定方差齊假定方差齊假定方假定方差不齊差不齊2022年7月4日星期一18時(shí)08分46秒方差齊性檢驗(yàn)方差齊性檢驗(yàn)2022年7月4日星期一18時(shí)08分46秒 2022年7月4日星期一18時(shí)08分46秒(二）、隨機(jī)區(qū)組設(shè)計(jì)的兩因素ANOVA 例例5 5- -6 6 為研究不同卡環(huán)對(duì)牙齒的固定效果，以為研究不同卡環(huán)對(duì)牙齒的固定效果，以1010顆取自新鮮尸體的牙齒

5、為觀察單位。每顆牙齒分顆取自新鮮尸體的牙齒為觀察單位。每顆牙齒分別采用別采用3 3種普通卡環(huán)、種普通卡環(huán)、RPIRPI卡環(huán)和卡環(huán)和Y Y型卡環(huán)進(jìn)行固定，型卡環(huán)進(jìn)行固定，測(cè)試卡環(huán)抗拉強(qiáng)度，結(jié)果見(jiàn)表測(cè)試卡環(huán)抗拉強(qiáng)度，結(jié)果見(jiàn)表5-5-8 8。試分析。試分析3 3種卡種卡環(huán)的固定效果有無(wú)差異？環(huán)的固定效果有無(wú)差異？ 2022年7月4日星期一18時(shí)08分46秒表5-8 不同大拉環(huán)的抗拉強(qiáng)度(牛頓)2022年7月4日星期一18時(shí)08分47秒block:block:區(qū)組變量；卡環(huán)類(lèi)別區(qū)組變量；卡環(huán)類(lèi)別: :分組變量，分組變量，1=1=普通卡環(huán)，普通卡環(huán)，2=2=RPIRPI卡環(huán)，卡環(huán)，3=3=Y Y型卡環(huán)；

6、型卡環(huán)； 抗拉強(qiáng)度抗拉強(qiáng)度: :（牛頓）。（牛頓）。2022年7月4日星期一18時(shí)08分47秒選擇模型選擇模型兩兩比較兩兩比較2022年7月4日星期一18時(shí)08分47秒主效應(yīng)主效應(yīng)全因素模型全因素模型自定義模型（不包括交互效應(yīng)）自定義模型（不包括交互效應(yīng)）2022年7月4日星期一18時(shí)08分47秒假定方差齊假定方差齊假定方假定方差不齊差不齊2022年7月4日星期一18時(shí)08分47秒?yún)^(qū)組區(qū)組處理處理誤差誤差. .SSSS. .MSMS. .F F值值. .P P值值. .2022年7月4日星期一18時(shí)08分47秒P P值值多重比較多重比較2022年7月4日星期一18時(shí)08分47秒Homogene

7、ous SubsetsSNKSNK抗拉強(qiáng)度抗拉強(qiáng)度卡環(huán)類(lèi)別NSubset12Student-Newman-Keulsa,bY型卡環(huán)107.1300普通卡環(huán)107.4700RPI卡環(huán)108.8100Sig.1791.000Means for groups in homogeneous subsets are displayed. Based on observed means. The error term is Mean Square(Error) = .296.a. Uses Harmonic Mean Sample Size = 10.000.b. Alpha = 0.05.相同列，無(wú)差別

8、相同列，無(wú)差別不同列，有差別不同列，有差別SNKSNK2022年7月4日星期一18時(shí)08分47秒（三）、析因設(shè)計(jì)資料的ANOVA 例例5-8 5-8 為研究藥物治療附加磁場(chǎng)對(duì)體內(nèi)磁性物質(zhì)為研究藥物治療附加磁場(chǎng)對(duì)體內(nèi)磁性物質(zhì)分布的影響，安排兩個(gè)藥物組：實(shí)驗(yàn)組：分布的影響，安排兩個(gè)藥物組：實(shí)驗(yàn)組：“絲裂絲裂霉素霉素+ +高分子物質(zhì)高分子物質(zhì)+ +磁性物質(zhì)磁性物質(zhì)+ +磁場(chǎng)磁場(chǎng)”，對(duì)照組：，對(duì)照組：“絲裂霉素絲裂霉素+ +高分子物質(zhì)高分子物質(zhì)+ +磁性物質(zhì)磁性物質(zhì)”，每組分別，每組分別于給藥后于給藥后1515分鐘和分鐘和6060分鐘處死實(shí)驗(yàn)動(dòng)物小鼠，檢分鐘處死實(shí)驗(yàn)動(dòng)物小鼠，檢測(cè)測(cè)2424只小鼠肝臟組

9、織磁性物質(zhì)鐵的濃度，分析附只小鼠肝臟組織磁性物質(zhì)鐵的濃度，分析附加磁場(chǎng)、藥物作用時(shí)間分別對(duì)小鼠肝臟組織鐵含加磁場(chǎng)、藥物作用時(shí)間分別對(duì)小鼠肝臟組織鐵含量是否有影響？附加磁場(chǎng)與藥物作用時(shí)間是否存量是否有影響？附加磁場(chǎng)與藥物作用時(shí)間是否存在交互效應(yīng)？在交互效應(yīng)？2022年7月4日星期一18時(shí)08分47秒 實(shí)驗(yàn)組（A1） 對(duì)照組（A2） 15min(B1) 60min(B2) 15min(B1) 60min(B2) 0.5540.5500.5780.7060.6860.6511.0151.0051.0711.1061.1551.1450.3370.2760.3130.3870.4310.3620.50

10、30.6120.5930.6040.6400.5603.7250.6212.33636.4971.0837.05562.1060.3510.75433.5120.5852.067215.840.49512.2134iTix21kijix 小鼠肝臟組織的鐵濃度檢測(cè)結(jié)果2022年7月4日星期一18時(shí)08分47秒a因素:藥物，1=實(shí)驗(yàn)組，2=對(duì)照組；b因素:給藥后時(shí)間，1=15min，2=60min； y:磁性物質(zhì)鐵的濃度（mg/g）2022年7月4日星期一18時(shí)08分47秒選擇模型選擇模型兩兩比較兩兩比較輪廓圖輪廓圖選項(xiàng)選項(xiàng)（均數(shù)及方差齊性檢驗(yàn)）（均數(shù)及方差齊性檢驗(yàn)）固定因素固定因素2022年7月

11、4日星期一18時(shí)08分47秒全模型全模型（同時(shí)分析主效應(yīng)和交互效應(yīng)）（同時(shí)分析主效應(yīng)和交互效應(yīng)）2022年7月4日星期一18時(shí)08分47秒橫軸橫軸區(qū)分線(xiàn)區(qū)分線(xiàn)2022年7月4日星期一18時(shí)08分47秒2022年7月4日星期一18時(shí)08分47秒2022年7月4日星期一18時(shí)08分47秒2022年7月4日星期一18時(shí)08分47秒2022年7月4日星期一18時(shí)08分47秒2022年7月4日星期一18時(shí)08分47秒（四）、重復(fù)測(cè)量資料的ANOVA例例5.9 5.9 臨床上為指導(dǎo)腦梗塞患者的治療和預(yù)后，某臨床上為指導(dǎo)腦梗塞患者的治療和預(yù)后，某研究人員對(duì)不同類(lèi)型腦梗塞患者酸性磷脂研究人員對(duì)不同類(lèi)型腦梗塞患

12、者酸性磷脂(AP)(AP)在不在不同時(shí)間點(diǎn)的變化，進(jìn)行了如下觀察：隨機(jī)選取三種同時(shí)間點(diǎn)的變化，進(jìn)行了如下觀察：隨機(jī)選取三種不同類(lèi)型的腦梗塞不同類(lèi)型的腦梗塞( (短暫性腦缺血發(fā)作、腦血栓形成短暫性腦缺血發(fā)作、腦血栓形成、腔隙性腦梗塞、腔隙性腦梗塞) )患者各患者各8 8例，于腦梗塞發(fā)生的第例，于腦梗塞發(fā)生的第2424小時(shí)、小時(shí)、4848小時(shí)、小時(shí)、7272小時(shí)、小時(shí)、7 7天分別采血，測(cè)量血中天分別采血，測(cè)量血中APAP值，結(jié)果見(jiàn)下表。值，結(jié)果見(jiàn)下表。2022年7月4日星期一18時(shí)08分47秒患者類(lèi)型病例編號(hào)測(cè)量時(shí)間( ) 24小時(shí)48小時(shí)72小時(shí)7天短暫性腦缺血發(fā)作1 4.8422.4419.

13、384.04 2 7.2615.8421.826.07 312.1014.52 9.685.28 411.3717.1611.792.34 532.3038.0026.408.21 6 5.88 7.92 6.393.59 7 4.7518.48 9.904.10 8 5.28 7.49 6.492.40腦血栓形成9 5.8313.5315.4311.88 10 6.3411.9726.4023.23 1123.3933.2629.7514.78 12 9.7318.5812.07 8.74 1310.5619.9618.48 7.92 1423.7626.0926.4010.98 1522.

14、6729.9733.5619.69 1617.1622.4425.86 3.17腔隙性腦梗塞1740.1312.3210.8919.71 1859.3851.7413.3417.95 1927.6424.1032.2156.33 2010.19 8.82 5.7210.03 2131.4829.0427.5312.17 2225.90 7.2034.60 7.10 2342.8031.00 8.4020.10 2423.9616.3920.3415.87不同類(lèi)型腦梗不同類(lèi)型腦梗塞塞患者患者APAP值值( (mol/L)mol/L)2022年7月4日星期一18時(shí)08分47秒leixing:1=leixing:1=短暫性腦缺血發(fā)作短暫性腦缺血發(fā)作 2= 2=腦血栓形成腦血栓形成 3= 3=腔隙性腦梗塞腔隙性腦梗塞t0=t0=第第2424小時(shí)小時(shí) t1=t1=第第4848小時(shí)小時(shí)t2=t2=第第7272小時(shí)小時(shí) t3=t3=第第7 7天天2022年7月4日星期一18時(shí)08分47秒20