1、程序： （yp111.sas）编程说明：在数据步中，变量x读取测定值，产生一个差值变量y（y=x-20.7），在过程步中，计算出Y的均值、标准差、标准误差、t值、P值。 结果说明：因t1.0636907，0.05p=0.3125，故用此法测定水中CaCO3的含量的均值与真值间无显著差异。此法可信。 例2：配对比较的t检验。研究食物中维生素E与肝脏中维生素A含量的关系。将大白鼠按性别、体重配对。每对随机分配，一个用正常饲料，一个用缺乏维生素E的饲料。经过一个时期饲养，杀死动物测定肝中维生素A的含量，结果如下表：大白鼠肝脏中维生素A含量（IU/g）配对号12345678正常饲料3550200030
2、0039503800375034503050缺乏E饲料24502400180032003250270025001750 （yp112.sas） 编程说明： 数据步中，把每对数据中的一个作为x，另一个作为y，计算出差数d（d=x-y），在过程步中，计算差数d的均值、标准差、标准误差、t值、P值。因t=4.21 p=0.00400.05，故有非常显著差异，即正常饲料组鼠肝维生素A含量比维生素E缺乏组的含量大。 例3 自身对照比较的t检验。应用克矽平治疗矽肺患者10名，治疗前后血红蛋白的含量如下表，问该药是否会引起血红蛋白的变化?治疗前后血红蛋白的含量（mg%）治疗前11.315.013.512.8
3、10.011.012.013.012.3治疗后14.013.814.711.4 （yp113.sas）因t=1.1989377，O.05F=右边的值，即方差相等假设检验的P值，如果此值大于0.01则在方差相等假设下继续进行，否则只能使用近似T检验。因P=0.0074，说明方差差异显著，即方差不齐。 当方差不齐时，T=2.0000，0.05P0.0644，故无显著差异，即此药物对大白鼠肉瘤无显著影响。第二节 方差分析当试验结果受到多个因素的影响，而且也受到每个因素的各水平的影响，为从数量上反映各因素以至各因素诸水平对试验结果的影响时使用方差分析的方法。 方差分析的基本思想是把全部数据关于总均值的
4、离差平方和分解成几个部分，每一部分表示某因素交互作用所产生的效应，将各部分均方与误差均方相比较，从而确认或否认某些因素或交互作用的重要性。用公式概括为： 总变异组间变异+组内变异其中：组间变异由各因素所引起，组内变异由个体差异所引起的，或者说由误差引起的。 常用的方差分析法有以下4种： （1）完全随机设计资料的方差分析（单因素方差分析） （2）随机区组设计资料的方差分析（二因素方差分析） （3）拉丁方设计资料的方差分析（三因素方差分析）（4）R*C析因设计资料的方差分析（有交互因素的方差分析） SAS系统中，ANOVA过程可以处理以上情形的方差分析，但它要求每个分类因子的组合观察数相等，即数据
5、是均衡的。若不均衡，就要求用GLM过程进行处理。在只考虑组间变异和误差变异时，称为单向方差分析。此时ANOVA会自动处理均衡和非均衡数据。在方差分析中，每次只研究1个指标时，称之为一元方差分析（简称ANOVA），同时考虑多个指标时，称之为多元方差分析（MANOVA）。在这一节里，我们还将讨论协方差分析。9.2.1 均衡数据的方差分析（ANOVA过程） PROC ANOVA 选项 CLASS 变量表； MODEL 因变量表=效应； MEANS 效应选择项； 2使用说明 （1）程序中，CLASS语句和MODEL语句是必需的，而且，CLASS语句必须出现在MODEL语句之前。（2）CLASS语句中的
6、变量是分类变量，可以是数值型，也可以是字符型。（3）MODEL语句指明因变量和自变量（因子变量）效应。效应是分类变量的各种组合，效应可以是主效应、交互效应、嵌套效应和混合效应。对应的效应模型如下： 主效应模型MODEL y=a b c； 模型中，a，b，c是主效应，y是因变量。下同。交互模型MODEL y=a b c a*b a*c b*c a*b*c； 模型中，a*b，a*c,b*c，a*b*c是交互效应。嵌套效应模型MODEL y=a b c（a b）； 模型中，c（a b）是嵌套效应。混合效应模型MODEL y=a b（a） c（a） b*c（a）； （4）MEANS语句是选择语句，计算
7、并输出所列的效应对应的因变量均值，若指明了选择项，则将进行主效应均值间的检验。常用的选择项如下：BON、DUNCAN、LSD、REGWF、REGWO、SNK（Q检验）、SCHEFFE、SIDAK、SMM（GT2）、TUKEY、WALLER。以上选择项在实际应用中，一般选择一种或两种方法即可。ALPHA=p 确定检验的显著性水平。缺省值是0.05。 3举例某劳动卫生研究所研究棉布、府绸、的确凉、尼龙四种衣料吸附十硼氢量。每种衣料各做五次测量，所得数据如下表。试检验各种衣料吸附十硼氢量有没有显著差别?各种衣料间棉花吸附十硼氢量棉布府绸的确凉尼龙2.332.483.064.002.002.345.1
8、32.932.683.004.612.732.662.802.223.60 （yp115.sas）数据步中，用循环控制变量a做分类变量，其水平数是4，分别代表不同的衣料。过程步中，用CLASS语句指明一个因素a，用MODEL语句反映出该因素的效果模型。结果说明：在输出中，找CLASS语句指出的变量的P值。此例中，P0.0003，可得出各衣料组间有非常显著差异。说明各种衣料间吸附十硼氢量是不同的。R-Square（R平方）对单向方差分析时，描述组间变异占总变异的比例，它越接近1，说明变异越归因于组间变异。 （2）随机区组设计资料的方差分析（两因素方差分析）用4种不同方法治疗8名患者，其血浆凝固时
9、间的资料如下表，试分析影响血浆凝固的因素。表：4种方法8名患者血浆凝固时间（分）资料受试者编号（区组）处 理 组8.49.49.812.215.212.914.49.69.111.28.89.98.28.58.610.98.99.09.210.47.98.1（yp116.sas） 数据步中，变量a代表不同治疗方法，其水平数是4，变量b代表区组因素不同患者，其水平数是8。过程步中，用CLASS语句指明两个因素a和b，用MODEL语句指明二因素的效果模型。 总处理间：F14.03，PO.0001，故总体有非常显著差异。 A因素：F=6.62，PO.0025，故认为因素A（治疗方法）对血浆凝固时间影
10、响很大。 B因素：F=17.20，P0.0001，故认为因素B（不同患者）对血浆凝固时间影响很大。 如果想对4种治疗方法进行两两比较，可在MODEL语句后增加一MEANS语句。例如，用SNK法比较，语句为：MEANS ASNK； （yp117.sas）若试验中涉及到3个因素，各因素间不存在交互作用，或交互作用很小可忽略不计，同时各因素的水平数又都相同，此种资料称拉丁方设计资料。例：五种防护服，由五人各在不同的五天中穿着测定脉博数，如下表所示。试比较五种防护服对脉搏数有无不同。五人穿着五种防护服的脉搏数试验日期受 试 者甲乙丙丁戊A 12.98B 116.2C 114.8D 104.0E 100
11、.6B 144.4C 119.2D 113.2E 132.8A 115.2C 143.0D 118.0E 115.8A 123.0B 103.8D 133.4E 110.8A 114.0B 98.0C 110.6E 142.8A 110.6B 105.8C 120.0D 109.8（yp118.sas） 数据步中，变量a代表试验时间因素，变量b代表受试者因素，变量C代表防护服因素。过程步中，用CLASS指明三个因素，MODEL语句反映三个因素的效果模型。 因F=6.80 P0.0011，故总体有非常显著差异，其中只有一个区组B因素（受试者）影响极大（F=16.27 P0.001），而区组A因素
12、（试验时间）和处理因素（防护服）都无显著差异。 （4）R*C析因设计资料的方差分析（有交互因素的方差分析）治疗缺铁性贫血病人12例，分为4组给予不同治疗，一个月后观察红细胞增加数（百万mm），资料如下表。试分析两种药物对红细胞增加的影响。贫血病人治疗后红细胞增加数甲药（A）不用（A0）用（A1）乙药（B）不用（B0）0.80.90.71.31.21.1用（B1）1.02.12.22.0（yp119.sas） 数据步中。变量a代表甲药因素变量b代表乙药因素。过程步中，用CLASS语句指明两个因素，用MODEL语句反映主效应（a b）和交互效应（a*b）模型。因F=98.75， P=0.00010
13、.01，故总体有非常显著差异，其中A因素（甲药）、B因素（乙药）、a*b因素（a和b的交互作用）都对红细胞增加数有非常大的影响。 为进一步分析其影响的强弱，可用MEANS语句： MEANS a b a*b： （yp120.sas）（5）用Dunnet法进行均值间的比较（对照组与多个实验组的比较检验）。Dunnett方法在比较若干“实验水平”与“对照水平”时使用一个临界差距（T）：，其中ms为将方差分析表中交互效应归入残差平方和后的均方，即：n为样本大小, q为来自Dunnett分布的分位数，可根据样本组个数r，自由度df（交互效应df+残差df），以及比较的显著水平，查表确定。Dunnett多
14、重比较仅适用各“实验水平”与“对照水平”之间差异的比较，不适用于各“实验水平”间的相互差异比较。有资料如下表，试比较白血病鼠与正常鼠脾中DNA平均含量（mgg）差别有无显著意义?表；白血病鼠与正常鼠脾中DNA平均含量（mgg）（3）环境影响评价中应用环境标准的原则。脾中DNA平均含量（mgg）A 正常鼠（对照组）3.意愿调查评估法12.32.间接市场评估法13.213.7环境的两个特点：第一节环境影响评价15.4（二） 环境影响经济损益分析的步骤15.816.9B1 自发性白血病鼠10.8既包括天然的自然环境，也包括人工改造后的自然环境。11.612.7（4）环境保护验收。14.8（3）安全现
15、状评价。B2 移植性白血病鼠10.3主体是人类；11.111.7（yp121.sas）数据步中，变量a代表因素。过程步中，用MODEL语句指明一因素的效果模型，用MEANS语句进行a因素水平间的Dunnett法比较。 因F=11.42，P0.0006，故总体或A因素各水平间有非常显著差异。当A因素有意义时，可进一步分析组间的差异。Dunnet法检验时，用“*”指明有显著性，由此可以看出，两组均值（A=2，A=3）与对照组（A=1）相比较，差别有显著意义。若改换对照组，假设是第2组（A=2）。应修改选择项DUNNETT（2）。输出中给出DUNNET t的临界值为2.399、A3与A1均值之差为-
16、1.9，其95置信限的下限和上限分别为-3.597和-O.203。 （6）用SNK法（Q检验）作均值间的两两比较。为研究酵解作用对血糖浓度的影响，从8名健康人中抽取了血液并制备成血滤液。每一个受试者的血滤液又分成4份，然后随机地把各份血滤液分别放置0、45、90、135分钟后测定其中血糖浓度（资料见下表）。试比较放置时间对血糖浓度有无影响。放置时间对血糖浓度的影响受试者编号（区组号）放置时间（分）459013595898394888410610597981021121011039280（yp122.sas）数据步中，变量a代表受试者因素，变量b代表放置时间因素。过程步中，作二因素的方差分析，用
17、CLASS指明二个因素a和b，用MODEL语句指明二因素的效果模型，用MEANS语句对因素b作SNK检验（Q检验）。 因F=43.66,P0.0001，故总体有非常显著差异。 因F=28.74, PO.0001，说明因素A（受试者）对血糖浓度有很大影响。 因F=78.49,P0.000l，说明因素B（放置时间）对血糖浓度有很大影响。 对B因素各组均值作SNK法比较检验:经SNK检验，可以看出：由于O和45分钟的均值间用同一字母A分组，故无显著差别（Means with the same letter are not significantly different），其余都有显著差别。9.2.2
18、. 非平衡数据的方差分析（GLM过程）在SASSTAT中GLM（General Linear Models）过程分析功能最多。它包括： 简单回归（一元回归） 加权回归 多重回归及多元回归 多项式回归 方差分析（尤其对不平衡设计资料更为有效） 偏相关分析 协方差分析 多元方差分析 反应面模型分析 重复测量方差分析这里，我们先介绍有关GLM过程在方差分析中的应用。GLM过程里，采用MODEL语句反映因变量与自变量（即效应）的模型，其形式是：模型说明模型类型MODEL Y=A B C主效应MODEL Y=A B A*B交互效应MODEL Y=A B A（B）嵌套效应MODEL Y1 Y2=A B多元
19、方差分析MODEL Y=A X协方差分析 其中A，B，C是分类变量；X，Y是连续型变量。前面介绍的ANOVA过程只能用于平衡设计资料（每个分类因子的组合观察数相等）的方差分析，当不平衡时，只能用GLM过程。 PROC GLM； MOIDEL 因变量一效应； MEANS 效应选择项；RUN； （1）程序中，CLASS语句和MODEL语句是必需的，而且CLASS语句必须出现在MODEL语句之前。 （2）CLASS语句中的变量是分类变量，可以是数值型，也可以是字符型。 （3）MODEL语句指明因变量和自变量（因子变量）效应，其效应可以是主效应、交互效应、嵌套效应和混合效应。对应的效应模型同ANOVA
20、过程中效应模式相同。 （4）MEANS语句是选择语句，计算所列的效应对应的因变量均值，若指明了选择项，则将进行主效应均值间的检验。 BON、DUNCAN、LSD、REGWF、REGWQ、SNK（Q检验）、SCHEF、SIDAK、SMM（GT2）、TUKEY、WALLER、DUNNETT（1-n检验） 以上选择项在实际应用中，一般选择一种或两种方法即可。ALPHA=p 给出检验的显著性水平。 （1）单因素方差分析健康男子各年龄组淋巴细胞转化率（）如下表，问各组淋巴细胞转化率的均值之间的差异是否显著?健康男子各年龄组淋巴细胞转化率（）1120岁58 61 61 62 63 68 70 70 74
21、784150岁54 57 57 58 60 60 63 64 666175岁43 52 55 56 60（yp123.sas）变量a代表年龄组因子，有3个水平值；变量n代表各组例数。过程步中，用CLASS语句指明因素为分类变量a，用MODEL语句指明单因素的效果模型。因F=9.77, P=0.0010.05，故总体或A因素各水平间有非常显著差异。即说明各年龄组的淋巴细胞转化率有显著差异。 （2）二因素方差方析对下表中的假设数据作二因素方差分析。因 素 BB1B2B3因素AA13.32.61.53.63.11.91.63.25.24.7A24.24.35.32.82.94.43.85.1A33.93.54.92.54.84.65.63.0 （yp124.sas） 数据步中，变量a代表因素a，其水平数是2，变量b代表因素b，其水平数是3，变量n是每组的重复次数。过程步中，用
