spss方差分析

1、方差分析王友信一教樓307Tel: 83911497E-mail: 實(shí)習(xí)目的 掌握單因素和雙因素方差分析的操作方法； 正確解釋輸出結(jié)果。 一、單因素方差分析一、單因素方差分析 例例1 1 將18只大鼠隨機(jī)分為三組，用二氧化硅（SiO2）50mg染塵，分別于染塵后1、3、6個(gè)月處死，稱量全肺濕重（見表4.2），試說明染塵后 1、3、6個(gè)月三個(gè)時(shí)期的全肺濕重是否有變化？表1 染塵后3個(gè)時(shí)期大鼠全肺濕重(g)觀測(cè)結(jié)果染塵時(shí)間： 1個(gè)月3個(gè)月6個(gè)月 3.34.43.6 3.64.44.4 4.33.45.1 4.14.25.0 4.24.75.5 3.34.24.7全肺濕重 Practice in S

2、PSS for Windows（一）（一） 建立數(shù)據(jù)文件建立數(shù)據(jù)文件 設(shè)變量設(shè)變量 “ “ Group ” Group ” 取值取值1 13 3，分別表示不同月份組，分別表示不同月份組 變量變量 “ “weight ” weight ” 各組全肺濕重測(cè)量值各組全肺濕重測(cè)量值Practice in SPSS for Windows( (二二) ) 統(tǒng)計(jì)分析統(tǒng)計(jì)分析 首先要對(duì)首先要對(duì)3 3組數(shù)據(jù)作正態(tài)性檢驗(yàn)組數(shù)據(jù)作正態(tài)性檢驗(yàn) 單因素方差分析及方差齊性檢驗(yàn)單因素方差分析及方差齊性檢驗(yàn)兩組獨(dú)立樣本兩組獨(dú)立樣本 t t 檢驗(yàn)，建檢驗(yàn)，建立數(shù)據(jù)文件的方法與方差立數(shù)據(jù)文件的方法與方差分析檢驗(yàn)有所異同？分析檢

3、驗(yàn)有所異同？。數(shù)據(jù)文件分析步驟分析步驟. .錄入數(shù)據(jù)，建立兩個(gè)檢驗(yàn)變量錄入數(shù)據(jù)，建立兩個(gè)檢驗(yàn)變量groupgroup和和weightweight，groupgroup取值取值1 1、2 2、3 3，分別代表染塵后分別代表染塵后1,3,61,3,6個(gè)月；個(gè)月；weight weight 為染塵后為染塵后3 3個(gè)時(shí)期大鼠全肺濕重個(gè)時(shí)期大鼠全肺濕重(g)(g) 。正態(tài)性檢驗(yàn)正態(tài)性檢驗(yàn)2. 首先對(duì)資料進(jìn)行正首先對(duì)資料進(jìn)行正態(tài)性檢驗(yàn)，方法同獨(dú)立態(tài)性檢驗(yàn)，方法同獨(dú)立樣本的樣本的T檢驗(yàn)檢驗(yàn): 選擇選擇DataSplit FileCompare groups把把“group”調(diào)入調(diào)入“Groups Based

4、 on”下，下，OK。AnalyzeNonparametric Tests1-Sample K-S將變量將變量weight移入移入Test Variable List框內(nèi)；框內(nèi)；在在Test Distribution中中激活激活“Normal”。單擊。單擊OK按鈕。則得出輸出按鈕。則得出輸出結(jié)果。結(jié)果。3. 3. 再回到再回到DataDataSplit FileSplit FileAnalyze all cases, do Analyze all cases, do notnotcreate groupscreate groupsOKOK。4.4.選擇選擇AnalyzeAnalyzeCompa

5、re MeansCompare MeansOne-Way ANOVAOne-Way ANOVA的的 主對(duì)主對(duì)話框。話框。5.5.將因變量將因變量weightweight移入移入Dependent List Dependent List 框內(nèi)?？騼?nèi)。GroupGroup變量變量進(jìn)入進(jìn)入 FactorFactor框中?？蛑小?. 6. 單擊單擊Option Option 按鈕，點(diǎn)擊按鈕，點(diǎn)擊“Descriptive”Descriptive”，即做基本，即做基本統(tǒng)計(jì)描述；點(diǎn)擊統(tǒng)計(jì)描述；點(diǎn)擊“Homogeneity of variance test”Homogeneity of variance te

6、st”， 即做方差齊性檢驗(yàn)；點(diǎn)擊即做方差齊性檢驗(yàn)；點(diǎn)擊ContinueContinue，回到主對(duì)話框，點(diǎn)，回到主對(duì)話框，點(diǎn)OKOK。輸出結(jié)果輸出結(jié)果7. 7. 如需要進(jìn)行如需要進(jìn)行各組均數(shù)的多各組均數(shù)的多重比較，在主重比較，在主對(duì)話中點(diǎn)擊對(duì)話中點(diǎn)擊Post HocPost Hoc按鈕。按鈕。進(jìn)入進(jìn)入Post Hoc Post Hoc Multiple Multiple ComparisonsComparisons對(duì)對(duì)話框。如果方話框。如果方差齊，常用差齊，常用LSDLSD、SNKSNK等方法對(duì)各等方法對(duì)各均值進(jìn)行比較，均值進(jìn)行比較，選中后，點(diǎn)擊選中后，點(diǎn)擊ContinueContinue，回，

7、回到主對(duì)話框，到主對(duì)話框，點(diǎn)擊點(diǎn)擊OKOK。輸出結(jié)果輸出結(jié)果Post Hoc TestsMultiple ComparisonsDependent Variable: WEIGHT-.2500.29926.417-.8879.3879-.7500*.29926.024-1.3879-.1121.2500.29926.417-.3879.8879-.5000.29926.115-1.1379.1379.7500*.29926.024.11211.3879.5000.29926.115-.13791.1379(J) GROUP3月6月1月6月1月3月(I) GROUP1月3月6月LSDMeanDi

8、fference(I-J)Std. ErrorSig.Lower BoundUpper Bound95% Confidence IntervalThe mean difference is significant at the .05 level.*. LSD LSD 結(jié)果結(jié)果輸出結(jié)果輸出結(jié)果SNKSNK結(jié)果結(jié)果WEIGHT63.966764.216764.7167.059GROUP1月3月6月Sig.Student-Newman-KeulsaN1Subsetfor alpha= .05Means for groups in homogeneous subsets are displayed.

9、Uses Harmonic Mean Sample Size = 6.000.a. Homogeneous Subsets二、隨機(jī)區(qū)組設(shè)計(jì)的方差分析二、隨機(jī)區(qū)組設(shè)計(jì)的方差分析窩別窩別 T1T1 T2T2 T3 T3 1 140.640.678.378.3116.3116.32 244.844.886.086.0124.6124.63 336.736.772.172.1149.0149.04 449.949.995.495.4128.8128.85 559.859.899.299.2134.1134.16 654.554.595.995.9133.0133.07 738.438.476.476.

10、4115.6115.68 841.641.679.979.9117.0117.09 946.846.886.586.5128.4128.4101044.744.785.385.3124.3124.3表表14-20 14-20 三組大鼠三組大鼠MTMT含量值（含量值（g/gg/g） 為探討為探討RglRgl對(duì)鎘誘導(dǎo)大鼠睪丸損傷的保護(hù)作用，將同一窩別的對(duì)鎘誘導(dǎo)大鼠睪丸損傷的保護(hù)作用，將同一窩別的3 3只大鼠隨機(jī)地只大鼠隨機(jī)地分到三組，進(jìn)行不同處理（對(duì)照組、氯化鎘組、分到三組，進(jìn)行不同處理（對(duì)照組、氯化鎘組、RglRgl+ +氯化鎘組）氯化鎘組）, , 共觀察了共觀察了1010個(gè)窩別大鼠的睪丸個(gè)窩別

11、大鼠的睪丸MTMT含量（含量（g/gg/g），數(shù)據(jù)如下表所示，試問不同處理對(duì)大），數(shù)據(jù)如下表所示，試問不同處理對(duì)大鼠鼠MTMT含量有無影響？含量有無影響？ （一）建立數(shù)據(jù)文件（一）建立數(shù)據(jù)文件 取變量取變量 group group 表示表示 T1 、 T2 、 T3三個(gè)水平組，三個(gè)水平組， block block 表示表示10 10 個(gè)不同窩別，個(gè)不同窩別， mtmt 表示表示 MT MT 含量值。含量值。 Practice in SPSS for Windows（二）統(tǒng)計(jì)分析（二）統(tǒng)計(jì)分析1. 1. 對(duì)對(duì) T T1 1 、 T T2 2 、 T T3 3 三個(gè)不同處理的數(shù)據(jù)，作正態(tài)性檢驗(yàn)三個(gè)

12、不同處理的數(shù)據(jù)，作正態(tài)性檢驗(yàn) 2. 作隨機(jī)區(qū)組設(shè)計(jì)的方差分析作隨機(jī)區(qū)組設(shè)計(jì)的方差分析數(shù)據(jù)文件錄入數(shù)據(jù)正態(tài)性檢驗(yàn)2. 2. 對(duì)資料進(jìn)行正態(tài)性檢驗(yàn)，方法同前。對(duì)資料進(jìn)行正態(tài)性檢驗(yàn)，方法同前。正態(tài)性檢驗(yàn)結(jié)果One-Sample Kolmogorov-Smirnov Test1045.78007.25684.154.154-.105.486.9721085.50009.07206.162.156-.162.514.95510127.110010.18162.146.146-.129.462.983NMeanStd. DeviationNormal Parametersa,bAbsolutePositi

13、veNegativeMost ExtremeDifferencesKolmogorov-Smirnov ZAsymp. Sig. (2-tailed)NMeanStd. DeviationNormal Parametersa,bAbsolutePositiveNegativeMost ExtremeDifferencesKolmogorov-Smirnov ZAsymp. Sig. (2-tailed)NMeanStd. DeviationNormal Parametersa,bAbsolutePositiveNegativeMost ExtremeDifferencesKolmogorov-

14、Smirnov ZAsymp. Sig. (2-tailed)GROUP123MTTest distribution is Normal.a. Calculated from data.b. 3.3.選擇選擇AnalyzeAnalyzeGeneral Linear ModelGeneral Linear ModelUnivariateUnivariate 展開展開UnivariateUnivariate主對(duì)話框。主對(duì)話框。4. 4. 定義定義mtmt為因變量，移入為因變量，移入Dependent Variable Dependent Variable 框內(nèi)。定義框內(nèi)。定義blockblock和

15、和groupgroup為固定變量，進(jìn)入為固定變量，進(jìn)入Fixed FactorsFixed Factors框中?？蛑?。5.5.單擊單擊ModelModel按鈕按鈕. .選擇選擇CustomCustom項(xiàng)，在項(xiàng)，在Build TermsBuild Terms欄內(nèi)的參數(shù)框中選擇欄內(nèi)的參數(shù)框中選擇Main Main effecteffect項(xiàng)，從項(xiàng)，從Factor & CovariatesFactor & Covariates框中分別將框中分別將blockblock和和groupgroup選入選入ModelModel框中?？蛑小?.6.再擊再擊ContinueContinue，回到主

16、對(duì)話框，單擊，回到主對(duì)話框，單擊OkOk按鈕，執(zhí)行多元按鈕，執(zhí)行多元方差分析過程。方差分析過程。方差分析結(jié)果方差分析結(jié)果第一欄為變異來源：第一欄為變異來源：BlockBlock不同窩別為配伍組，不同窩別為配伍組， groupgroup不同不同營(yíng)養(yǎng)素為處理組，營(yíng)養(yǎng)素為處理組，ErrorError為誤差，為誤差，Corrected TotalCorrected Total為總變?yōu)榭傋儺?。后面各欄為各個(gè)便以相應(yīng)的離均差平方和、自由度、異。后面各欄為各個(gè)便以相應(yīng)的離均差平方和、自由度、均方、均方、F F值、概率。值、概率。Tests of Between-Subjects EffectsDepende

17、nt Variable: MT34355.761a113123.25164.567.000222551.3071222551.3074600.797.0001276.9639141.8852.933.02533078.798216539.399341.919.000870.7021848.372257777.7703035226.46329SourceCorrected ModelInterceptBLOCKGROUPErrorTotalCorrected TotalType III Sumof SquaresdfMean SquareFSig.R Squared = .975 (Adjust

18、ed R Squared = .960)a. 7.7.如需進(jìn)行個(gè)組均數(shù)間的多重比較，則在主對(duì)話框中單如需進(jìn)行個(gè)組均數(shù)間的多重比較，則在主對(duì)話框中單擊擊Post HocPost Hoc按鈕。進(jìn)入按鈕。進(jìn)入Post Hoc Multiple Post Hoc Multiple ComparisonsComparisons對(duì)話框。方法同單因素方差分析。對(duì)話框。方法同單因素方差分析。三、三、 2 22 2 析因設(shè)計(jì)的方差分析析因設(shè)計(jì)的方差分析B1B2BA1A2A37.4337.4335.0435.0432.0232.0237.9437.9439.2139.2139.5839.58 33.3033.30

19、36.8336.8336.8836.8839.7739.7733.6333.6332.6232.62 30.5830.5838.1938.1936.3236.3243.5043.5027.4127.4139.3039.30 36.6036.6046.4346.43 30.5730.5736.2736.2737.7637.7645.1145.1154.1554.1551.6151.6130.6330.6335.9535.9533.9533.9552.4552.4546.7746.7747.6647.6629.6929.6936.0636.0632.4432.4454.9954.9950.6550

20、.6542.7142.7134.6934.6948.6048.60 40 40只只5 5周齡雌性周齡雌性BALB/cBALB/c小鼠被隨機(jī)分為小鼠被隨機(jī)分為4 4組組, ,接受不同的處接受不同的處理：理： 、 、 、 ，觀測(cè)小鼠，觀測(cè)小鼠BALFBALF中中IL-4IL-4（pg/mlpg/ml）水平，觀）水平，觀測(cè)結(jié)果見下表：測(cè)結(jié)果見下表：11BA21BA12BA22BA （一）建立數(shù)據(jù)文件（一）建立數(shù)據(jù)文件 本例要建立三個(gè)變量，本例要建立三個(gè)變量， B BB B1 1 表示用表示用 藥藥 B B2 2 表示用表示用 藥藥 A AA A1 1 表示用表示用 藥藥 A A2 2 表示用表示用

21、藥藥 IL-4 IL-4 表示表示IL-4IL-4（pg/mlpg/ml）水平。）水平。1B2B1A2APractice in SPSS for Windows（二）統(tǒng)計(jì)分析（二）統(tǒng)計(jì)分析數(shù)據(jù)文件1在評(píng)價(jià)某藥物耐受性及安全性的I期臨床試驗(yàn)中，對(duì)符合納入標(biāo)準(zhǔn)的40名健康自愿者隨機(jī)分為4組，每組10名，各組注射劑量分別為0.5U、1U、2U、3U，觀察48小時(shí)后部分凝血活酶時(shí)間（s）。試比較任意兩兩劑量間的部分凝血活酶時(shí)間有無差別？0.5 U1 U2 U3 U36.840.032.933.034.435.537.930.734.336.730.535.335.739.331.132.333.240.134.737.431.136.837.639.134.333.440.233.529.838.338.136.635.438.432.432.031.239.835.633.8表：各劑量組48小時(shí)部分凝血活酶時(shí)間（s）2為探討小劑量地塞米松對(duì)急性肺損傷動(dòng)物模型肺臟的