方差分析的SAS操作.docx

方差分析

常用于方差分析的主要过程有ANOVA和GLM（广义线性模型），对于平衡数据资料（各水平下等重复，数据没有丢失），一般用ANOVA过程，对于非平衡数据，应采用GLM过程.

1、 ANOVA过程格式及使用说明过程格式：

PROCANOVA[选项]；

CLASS 处理因素；

MODEL因变量=效应表[/选择项]；

MEANS 效应表[/选择项]；过程说明：

PROCANOVA 语句的选项主要有：

DATA=数据集名 指明要分析的SAS数据集，缺省时SAS将使用最近建立的数据集.

OUTSTAT=输出数据集 指定分析计算结果输出的数据集名.

CLASS

语句指明分类变量，是ANOVA过程的必需语句，并且

必须出现在MODEL语句之前.分类变量可以为数值型或字符型，分类变量的个数表示方差分析的因素个数.

MODEL语句定义分析所用的效应模型，即方差分析的因变量和效应变量.在方差分析过程中，关键在于定义线性数学模型，常用的模型定义语句有：

MODELy=a 单因素一元方差分析

MODELy=ab双因素无交互作用一元方差分析

MODELy=aba*b双因素有交互作用一元方差分析

MEANS(格式：因素/选择项)语句用来计算该语句所列的每个效应所对应的因变量均值，其选项用于设定多重比较的方法以及方差齐性检验。

可以选择的检验方法有：

T/LSD法：对means语句中出现的所有因素的各水平进行两两T检验，当每一水平的观测数相等时，T检验变成Fisher的最小显著差检验。

BON法：对MEANS语句中出现的所有因素的各水平均值之差进行Bonferroni的T检验。

TUKEY法：对MEANS语句中出现的所有因素的各水平均值进行TUKEY的学生化极差检验。

DUNCAN法：对MEANS语句中出现的所有因素的各水平均值进行DUNCAN的极差检验。

REGWF法：对MEANS语句中出现的所有因素的各水平均值进行多重极差检验。

HOVTEST：进行方差齐性检验。

2、GLM即广义线性模型（GeneralLinerModel）过程，它使用最小二乘法对数据拟合广义线性模型.该过程功能强大，可用于多种不同的统计分析中.GLM过程用于方差分析时，主要语句和使用格式与上述ANOVA过程类似，详见例3.

一、单因素方差分析1．单因素等重复方差分析(ANOVA过程)应用举例:

编写程序如下：

Data exam; /*建立数据集*/

Do trt=1to3;Do I=1to5;Inputx@@;Output;

End;End;Cards;

4046384244

2634302832

3940434850

;

Procsort;Bytrt;Run;

Procunivariatenormal;Varx;

/*3个水平（trt）分别为1、2、3*/

/*每个水平下5次重复*/

Bytrt;Run;

Procanova; /*调用方差分析过程*/

Classtrt; /*定义处理为分类变量*/

Modelx=trt; /*定义效应模型*/

Meanstrt/thovtestbon; /*要求计算每一水平下的均值，进行方差齐性检验，多重比较T检验和BON检验*/

Run;

Title方差分析;Run;

补充：‘试验错误率MEER’与‘比较错误率CER’。例如:

有5种施肥方法，均值做两两比较有10种组合，如要控制10种比较的总错判率，就称为“试验错误率”。

如果想对10种比较中的每一种都控制比较错判机会，就称为“比较错误率”。

例2 在4种不同的肥料处理（k1，k2，k3，k4）下测得土

壤的含氮量（nit），每个处理下重复5次，分析各处理间土壤含氮量的均值是否有显著差异（数据见程序中）.

DmLOG;CLEAR;OUTPUT;CLEAR;;

Datanew2;

Inputstr$nit@@; /*str表示处理，是字符型变量*/Cards; /*每一个数据前的符号ki(i=1,2,3,4)表明

该数据属于哪一个处理*/

k119.4k132.6k127k132k133

k217.7k224.8k227.9k225.2k224.3

k317k319.4k39.1k311.9k315.8

k420k421k420.5k418.8k418.6

;

Procprint;

Ti