方差分析操作课件

SPSS

在医学研究数据管理与分析中的应用(二)

主要内容描述性统计分析统计推断t检验、方差分析卡方检验一.描述性统计分析资料类型计量资料计数资料等级资料(2)计数资料以定性方式表达每个观察单位的某项观察指标，如血型、性别等。特点：①各观察单位间或者相同，或者存在质的差别；②有质的差别者之间无连续性。

(3)等级资料以等级表达每个观察单位的某项观察指标，如疗效分级、心功能分级等。特点：①各观察单位间或者相同，或者存在质的差别；②各等级间只有顺序，而无数值大小，故等级之间不可度量。资料类型的判断病例号年龄(岁)性别身高(cm)血型心电图尿WBC职业RBC1012/L135女1.65A正常－教师4.67244男1.74B正常－工人5.21326男1.80O正常＋职员4.10425女1.61AB正常＋农民3.92541男1.71A异常++工人3.49645女1.58B正常++工人5.48750女1.60O异常++干部6.78828男1.76AB正常+++干部7.10931女1.62O正常＋军人5.24计量资料的统计描述描述集中趋势的指标：众数、中位数、算术均数和几何均数。描述离散趋势的指标：极差、四分位间距、方差、标准差和变异系数。描述数据分布形态的指标：偏度和峰度数据特征描述的图示法：直方图、箱图正态分布用算术均数+标准差描述，偏态分布用中位数+四分位间距描述等级资料的统计描述频数表列联表表4-1100包药的重量（g）0.890.860.880.920.980.890.910.950.850.920.890.970.860.920.870.900.930.910.880.910.860.990.850.890.821.030.930.810.960.920.950.880.900.840.870.980.880.850.860.910.900.930.950.920.950.860.870.920.870.940.950.820.840.800.940.860.920.860.870.930.970.910.880.920.890.890.870.930.910.980.880.900.920.870.880.950.940.890.780.840.880.870.940.900.960.980.890.920.901.060.870.910.870.840.891.000.940.900.870.92AnalyzeDescriptiveStatisticsFrequencies集中趋势离散趋势分布百分位数直方图输出结果：统计量、直方图例4-2

从某医院随机抽取100名护士，调查其学历，数据见表4-2，试进行频数分布分析，列出频数表。输出结果Descriptives——描述变量集中趋势和离散趋势

AnalyzeDescriptiveStatisticsDescriptivesExplore——探索性分析统计图的绘制茎叶图直方图箱图正态性检验方差齐性检验AnalyzeDescriptiveStatisticsExplore判断有无离群点或极端值提示数据转换方法了解组间差异的特征Explore——探索性分析例4-3

为考察中药葛根对心脏功能的影响，配制每100ml含葛根1g，1.5g，3g，5g的药液，用来测定27只大鼠离体心脏在药液中7－8分钟时间内心脏冠脉血流量，数据见表4-2。试对不同剂量的葛根的心脏冠脉血流量作探索性分析。直方图箱式图正态性方差齐性输出结果：Descriptives表：4组统计量TestofNormality表：正态性结果输出结果：TestofHomogeneityofVariance

方差齐性检验Crosstabs——列联表资料分析生成列联表R*C表的统计分析AnalyzeDescriptiveStatisticsCrosstabsCrosstabs——列联表资料分析例4-4

肿瘤转移病人全量放疗后用两种药物（试验药、对照药）配合治疗，观察肿瘤消除情况，结果见表4-3，试生成药物和消除情况各水平组合的频数分布表，并计算两种药物的全消率。Crosstabs——列联表资料分析输出结果：Crosstabulation输出结果：Chi-SquareTests计数资料统计描述图示法饼图、直条图例16-4

某医院护士的学历构成的数据文件例4-5.sav，绘制护士的学历构成饼图与直条图。Graphs输出结果——饼图输出结果——直条图主要内容——DescriptiveStatisticsFrequenciesDescriptivesExploreCrosstabsAnalyze

DescriptiveStatistics二、统计推断研究设计配对设计：自身配对：同一受试对象处理前后或不同部位测定值的比较；同一受试对象(或样品)分别接受两种不同的处理。两两配对：将条件近似的观察对象两两配成对子，对子中的两个个体分别给予不同的处理。成组设计：将实验对象完全随机地分配到实验组与对照组或几个对比组中随机化随机抽样随机分组随机排序1.设定种子数TransformRandomnumbergenerate2.利用随机函数生成随机数TransformComputeVariable(选择函数RandomnumbersRV.Uniform(0,1))1、随机抽样例：某医院有护士100人，今欲了解其近视率，从中抽取20%即20人作调查。先将全院护士从1到100编号，然后用计算机产生100个0~1之间的随机数，按随机数排序，前20位随机数字对应用编号的护士组成本次调查的样本。2、随机分组例：欲将符合试验要求的实验动物20只，随机分配为A，B组。先将动物1到20编号，取种子数20030514，用计算机产生20个0~1之间的随机数，前10个对应于A组，后10个对应于B组。按随机数的大小进行排序，随机数的秩次就是动物的编号，再按动物的编号从小到大排列，对应的组别即为分配方案。方法的选择计量资料计数资料等级资料多个样本均数的比较2个样本均数的比较方差分析秩和检验卡方检验方差齐方差不齐t检验t′检验配对t检验成组设计配对设计成组设计配对设计配对卡方检验成组设计配对设计配对秩和检验计量资料的假设检验一、t检验1、样本均数与总体均数的比较2、完全随机设计（成组设计）的两样本均数的比较3、配对设计的两样本均数比较二、方差分析1、完全随机设计（成组设计）的单因素方差分析2、随机区组设计（配伍设计）的两因素方差分析样本均数与已知总体均数的比较例：通过大量调查,已知某地正常男婴体重为3.26Kg,某医生随机抽取20名难产男婴，测得出生体重如下：问，该地难产男婴体重是否不同于本地正常男婴？3.53.53.23.53.33.03.33.23.42.73.43.63.52.83.42.93.53.54.04.0

SPSS分析命令：Analyze→CompareMeans→One-SampleTTest…样本均数与已知总体均数的比较总体均数样本均数与已知总体均数的比较结论：不能认为难产男婴的出生体重与正常男婴不同。均数标准差标准误P

值自由度t值完全随机设计（成组设计）两个样本均数的比较例：某医生测得12名正常人和13名病毒性肝炎患者血清转铁蛋白含量（g/L）结果如下：问肝炎患者和正常人血清转铁蛋白含量有无差异？病毒性肝炎患者：2.342.472.222.312.362.382.152.572.192.252.282.312.42正常人：2.612.712.732.642.682.812.762.552.912.852.712.64

SPSS分析命令：Analyze→CompareMeans→Independent-SamplesTTest…完全随机设计（成组设计）两个样本均数的比较分析变量分组变量完全随机设计（成组设计）两个样本均数的比较统计描述方差齐性检验结论：病毒性肝炎患者与正常人血清转铁蛋白含量有差别。配对设计的两个样本均数比较例：为比较某新药与常规药降血脂的效果，将性别相同、血清总胆固醇水平相近的高血脂患者配成对子，每对中随机抽取一个人服用新药，另一个人服用常规药。服用一段时间后，测得血清总胆固醇含量（mmol/L）如下：问新药与常规药降血清总胆固醇效果是否相同？新药：6.576.466.276.896.217.617.607.046.687.42常规药：6.006.835.977.286.306.647.387.006.037.22

SPSS分析命令：Analyze→CompareMeans→Paired-SamplesTTest…配对设计的两个样本均数比较配对设计的两个样本均数比较结论：不能认为新药与常规药降低血清总胆固醇的效果不同。配对差值方差分析AnalysisofVariance,ANOVA方差分析的目的：推断多个总体均数是否相等方差分析的使用条件：各处理组样本来自正态总体各样本是相互独立的随机样本各处理组的总体方差相等，即方差齐性方差分析的结果解释：方差分析的F检验，当P≤0.05，可以认为各组总体均数不等或不全相等，但并不以为着任意两组总体均数都有差别。需要进一步作两两比较才能知道哪些组间有差别。多个样本均数间的两两比较当样本组数大于2时，不宜再用前述的t检验分别作两两比较，否则会增大犯第一类错误的概率。方差分析与t检验的联系

t检验可以看作时方差分析的特例:方差分析的注意事项完全随机设计（成组设计）的单因素方差分析1个因素，k个水平例4：某社区随机抽取30名糖尿病患者、糖耐量减低（IGT）者和正常人进行载脂蛋白(mg/dl)测定，结果如下：问三种人载脂蛋白水平有无差别？糖尿病患者：85.70105.20109.5096.00115.2095.30110.00100.00125.60111.00106.50IGT异常者：96.00124.50105.1076.4095.30110.0095.2099.00120.00正常人：144.00117.00110.00109.00103.00123.00127.00121.00159.00115.00

方法一：Analyze→CompareMeans→One-WayANOVA…完全随机设计（成组设计）的单因素方差分析1个因素，k个水平分组变量分析变量指定输出描述统计量指定进行方差齐性检验完全随机设计（成组设计）的单因素方差分析1个因素，k个水平方差齐性检验组间变异组内变异总变异SSMS结论：可以认为三种人血清载脂蛋白水平有差别。两两比较假定方差齐时假定方差不齐时两两比较多重比较结论：糖尿病患者与正常人、IGT异常者与正常人的载脂蛋白水平有差异；糖尿病患者与IGT异常者的载脂蛋白水平无差异。两两比较结论：糖尿病患者与正常人、IGT异常者与正常人的载脂蛋白水平有差异；糖尿病患者与IGT异常者的载脂蛋白水平无差异。SNK同一列，无差别不同列，有差别方法二：Analyze→GeneralLinearModel→Univariate…分析变量分组变量组间变异组内变异总变异除了方差分析表不同以外，方法一和方法二的其他输出结果是一致的：

SPSS分析命令：Analyze→GeneralLinearModel→Univariate…随机区组设计（配伍组设计）资料的方差分析（1个研究因素a个水平，1个配伍组因素b个水平）例：为探索丹参对肢体缺血再灌注损伤的影响，将30只纯种新西兰实验用大白兔，按窝别相同、体重相近划分为10个区组。每个区组3只大白兔随机采用A、B、C三种处理方案，即在松止血带前分别给予丹参2ml/kg、丹参1ml/kg、生理盐水2ml/kg，在松止血带前及松后1小时分别测定血中白蛋白含量（g/L），算出白蛋白减少量。问三种方案的处理效果是否不同？分组变量分析变量分析所有主效应或交互作用分析所有主效应和交互作用随机区组设计资料的方差分析组间区组组内结论：认为10个区组的总体均数相同。但是，A、B、C三种方案的处理效果不全相同，即三个总体均数中至少有两个不同。至于三个总体均数中哪些不同，同样需要进行多个均数间的两两比较。随机区组设计资料的方差分析结论：认为A、B种方案的处理效果没有差异，但均与C种处理方案的效果有差异，松止血带前给予丹参能使血中白蛋白含量降低。计数资料的假设检验和非参数检验一、分类变量的假设检验

1、完全随机设计（成组设计）的两样本率（构成比）比较2、配对设计的两样本率（构成比）比较3、完全随机设计的多个样本率（构成比）比较二、非参数检验

1、完全随机设计的两样本比较2、配对设计的两样本比较3、完全随机设计的多个样本比较4、随机区组设计多个样本比较完全随机设计两个样本率的比较两个样本率比较方法选择条件：n≥40，T≥5，用2；n≥40，但1≤T<5，用校正2。n<40，或T<1，用确切概率。完全随机设计两个样本率的比较例5-1：某中药在改变剂型前曾在临床观察152例，治愈129例，未治愈23例；改变剂型后又在临床观察130例，治愈101例，未治愈29例，能否得出新剂型疗效与旧剂型不同的结论？疗效合计治愈未治愈新剂型12923152旧剂型10129130合计23052282完全随机设计两个样本率的比较完全随机设计两个样本率的比较SPSS分析命令：Analyze→DescriptiveStatistics→Crosstabs…行变量列变量例5-2：某矿石粉厂当生产一种矿石粉时，在数天内即有部分工人患职业性皮炎，本生产季节开始，随机抽取15名车间工人穿上新防护服，其中1名患皮炎，其余28名工人仍穿旧防护服，其中10名患皮炎。生产一段时间后，检查两组工人皮炎患病率，问两组工人皮炎患病率有无差异？皮炎患病情况合计患病未患病新防护服11415旧防护服101828合计113243例5-3：某医师A药治疗9例病人，治愈7人；用B药治疗10例病人，治愈1人，问两药疗效是否有差别？疗效合计治愈未治愈A药729B药1910合计81119配对设计两个样本率的比较例6-1：分别用反向血凝法和酶标法对200名献血员进行HBsAg检测，结果如下，问两种检验方法检出率有无差异？反向血凝法合计+-酶标法+-3018482150152合计32168200方法一：Crosstabs方法一：Crosstabs方法一：Crosstabs方法一：Crosstabs方法二：非参数检验NonparametricTest方法二：非参数检验NonparametricTest方法二：非参数检验NonparametricTest精确概率法完全随机设计的多个样本率的比较例：某省从水氟含量不同的地区随机抽取10-12岁儿童，进行第一恒齿患病率的调查，问3个地区第一恒齿患病率是否不同？地区患齿数未患齿数合计高氟区452671干预区311546低氟区43548合计11946165SPSS分析命令：Analyze→DescriptiveStatistics→Crosstabs…完全随机设计的多个样本率的比较练习秩和检验主要适用于以下情况：数值变量资料：严重偏离正态分布或者分布未知数值变量资料：完全随机设计中方差不齐有序分类资料（等级资料）注：但数据满足参数检验的条件时