spss方差分析

1、方差分析王友信一教楼307Tel: 83911497E-mail: 实习目的 掌握单因素和双因素方差分析的操作方法； 正确解释输出结果。 一、单因素方差分析一、单因素方差分析 例例1 1 将18只大鼠随机分为三组，用二氧化硅（SiO2）50mg染尘，分别于染尘后1、3、6个月处死，称量全肺湿重（见表4.2），试说明染尘后 1、3、6个月三个时期的全肺湿重是否有变化？表1 染尘后3个时期大鼠全肺湿重(g)观测结果染尘时间： 1个月3个月6个月 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（一）（一） 建立数据文件建立数据文件 设变量设变量 “ “ Group ” Group ” 取值取值1 13 3，分别表示不同月份组，分别表示不同月份组 变量变量 “ “weight ” weight ” 各组全肺湿重测量值各组全肺湿重测量值Practice in SPSS for Windows( (二二) ) 统计分析统计分析 首先要对首先要对3 3组数据作正态性检验组数据作正态性检验 单因素方差分析及方差齐性检验单因素方差分析及方差齐性检验两组独立样本两组独立样本 t t 检验，建检验，建立数据文件的方法与方差立数据文件的方法与方差分析检验有所异同？分析检

3、验有所异同？。数据文件分析步骤分析步骤. .录入数据，建立两个检验变量录入数据，建立两个检验变量groupgroup和和weightweight，groupgroup取值取值1 1、2 2、3 3，分别代表染尘后分别代表染尘后1,3,61,3,6个月；个月；weight weight 为染尘后为染尘后3 3个时期大鼠全肺湿重个时期大鼠全肺湿重(g)(g) 。正态性检验正态性检验2. 首先对资料进行正首先对资料进行正态性检验，方法同独立态性检验，方法同独立样本的样本的T检验检验: 选择选择DataSplit FileCompare groups把把“group”调入调入“Groups Based

4、 on”下，下，OK。AnalyzeNonparametric Tests1-Sample K-S将变量将变量weight移入移入Test Variable List框内；框内；在在Test Distribution中中激活激活“Normal”。单击。单击OK按钮。则得出输出按钮。则得出输出结果。结果。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的的 主对主对话框。话框。5.5.将因变量将因变量weightweight移入移入Dependent List Dependent List 框内。框内。GroupGroup变量变量进入进入 FactorFactor框中。框中。6. 6. 单击单击Option Option 按钮，点击按钮，点击“Descriptive”Descriptive”，即做基本，即做基本统计描述；点击统计描述；点击“Homogeneity of variance test”Homogeneity of variance te

6、st”， 即做方差齐性检验；点击即做方差齐性检验；点击ContinueContinue，回到主对话框，点，回到主对话框，点OKOK。输出结果输出结果7. 7. 如需要进行如需要进行各组均数的多各组均数的多重比较，在主重比较，在主对话中点击对话中点击Post HocPost Hoc按钮。按钮。进入进入Post Hoc Post Hoc Multiple Multiple ComparisonsComparisons对对话框。如果方话框。如果方差齐，常用差齐，常用LSDLSD、SNKSNK等方法对各等方法对各均值进行比较，均值进行比较，选中后，点击选中后，点击ContinueContinue，回，

7、回到主对话框，到主对话框，点击点击OKOK。输出结果输出结果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 结果结果输出结果输出结果SNKSNK结果结果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二、随机区组设计的方差分析二、随机区组设计的方差分析窝别窝别 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对镉诱导大鼠睾丸损伤的保护作用，将同一窝别的对镉诱导大鼠睾丸损伤的保护作用，将同一窝别的3 3只大鼠随机地只大鼠随机地分到三组，进行不同处理（对照组、氯化镉组、分到三组，进行不同处理（对照组、氯化镉组、RglRgl+ +氯化镉组）氯化镉组）, , 共观察了共观察了1010个窝别大鼠的睾丸个窝别

11、大鼠的睾丸MTMT含量（含量（g/gg/g），数据如下表所示，试问不同处理对大），数据如下表所示，试问不同处理对大鼠鼠MTMT含量有无影响？含量有无影响？ （一）建立数据文件（一）建立数据文件 取变量取变量 group group 表示表示 T1 、 T2 、 T3三个水平组，三个水平组， block block 表示表示10 10 个不同窝别，个不同窝别， mtmt 表示表示 MT MT 含量值。含量值。 Practice in SPSS for Windows（二）统计分析（二）统计分析1. 1. 对对 T T1 1 、 T T2 2 、 T T3 3 三个不同处理的数据，作正态性检验三个

12、不同处理的数据，作正态性检验 2. 作随机区组设计的方差分析作随机区组设计的方差分析数据文件录入数据正态性检验2. 2. 对资料进行正态性检验，方法同前。对资料进行正态性检验，方法同前。正态性检验结果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主对话框。主对话框。4. 4. 定义定义mtmt为因变量，移入为因变量，移入Dependent Variable Dependent Variable 框内。定义框内。定义blockblock和

15、和groupgroup为固定变量，进入为固定变量，进入Fixed FactorsFixed Factors框中。框中。5.5.单击单击ModelModel按钮按钮. .选择选择CustomCustom项，在项，在Build TermsBuild Terms栏内的参数框中选择栏内的参数框中选择Main Main effecteffect项，从项，从Factor & CovariatesFactor & Covariates框中分别将框中分别将blockblock和和groupgroup选入选入ModelModel框中。框中。6.6.再击再击ContinueContinue，回到主

16、对话框，单击，回到主对话框，单击OkOk按钮，执行多元按钮，执行多元方差分析过程。方差分析过程。方差分析结果方差分析结果第一栏为变异来源：第一栏为变异来源：BlockBlock不同窝别为配伍组，不同窝别为配伍组， groupgroup不同不同营养素为处理组，营养素为处理组，ErrorError为误差，为误差，Corrected TotalCorrected Total为总变为总变异。后面各栏为各个便以相应的离均差平方和、自由度、异。后面各栏为各个便以相应的离均差平方和、自由度、均方、均方、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.如需进行个组均数间的多重比较，则在主对话框中单如需进行个组均数间的多重比较，则在主对话框中单击击Post HocPost Hoc按钮。进入按钮。进入Post Hoc Multiple Post Hoc Multiple ComparisonsComparisons对话框。方法同单因素方差分析。对话框。方法同单因素方差分析。三、三、 2 22 2 析因设计的方差分析析因设计的方差分析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小鼠被随机分为小鼠被随机分为4 4组组, ,接受不同的处接受不同的处理：理： 、 、 、 ，观测小鼠，观测小鼠BALFBALF中中IL-4IL-4（pg/mlpg/ml）水平，观）水平，观测结果见下表：测结果见下表：11BA21BA12BA22BA （一）建立数据文件（一）建立数据文件 本例要建立三个变量，本例要建立三个变量， 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（二）统计分析（二）统计分析数据文件1在评价某药物耐受性及安全性的I期临床试验中，对符合纳入标准的40名健康自愿者随机分为4组，每组10名，各组注射剂量分别为0.5U、1U、2U、3U，观察48小时后部分凝血活酶时间（s）。试比较任意两两剂量间的部分凝血活酶时间有无差别？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小时部分凝血活酶时间（s）2为探讨小剂量地塞米松对急性肺损伤动物模型肺脏的