生物统计常用sas程序

《生物统计常用sas程序》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、生物统计学常用的SAS程序62目录第二章几种常用的SAS程序§2.1利用SAS软件描述样本数据2.1.1用MEANS语句描述数据2.1.2频数分布表的编制2.1.3频数分布图的绘制§2.2统计假设检验的SAS程序2.2.1单个样本的t检验2.2.2配对数据t检验2.2.3成组数据t检验§2.3χ2检验的SAS程序§2.4单因素方差分析的SAS程序§2.5常用实验设计方差分析的SAS程序2.5.1三因素交叉分组实验的方差分析2.5.2随机化完全区组实验的方差分析2.5.3拉丁方与正交拉丁方实验的方差分析2.5.4两因素随机化

2、区组实验的方差分析2.5.5裂区实验设计的方差分析2.5.6套设计的方差分析2.5.7正交设计的方差分析§2.6相关与回归分析的SAS程序2.6.1一元线性回归分析2.6.2一元线性回归的图形描述2.6.3一元非线性回归2.6.4相关系数和偏相关系数的计算§2.7多元回归分析的SAS程序2.7.1多元回归方程计算2.7.2逐步回归分析62第二章几种常用的SAS程序§2.1利用SAS软件描述样本数据用SAS可以对样本数据进行全面描述，得出样本的各种特征数以及频数分布图。在阅读以下内容之前请先阅读第一章“SAS软件基本操作”。

3、2.1.1用MEANS语句描述数据例2.1计算课本上习题1.2的平均数和标准差。解在进行分析之前，最好先建立一个外部数据文件。如果不建立外部数据文件，在作业流中输入数据也可以。一般来说，在数据量比较大，有可能重复使用时，最好建立一个数据文件。在这里我们创建一个称为2-1data.dat的外部数据文件，存储在A盘中。用MEANS过程描述数据的最基本的程序如下：optionslinesize=76;dataabc;infile‘a:2-1data.dat’;inputx;run;procmeans;run;提交SAS运行后，

4、所得结果如下：表2－1：例2.1输出的结果TheSASSystem AnalysisVariable:XNMeanStdDevMinimumMaximum25063.27600003.013994155.000000070.000000062打印的结果中包括：数据个数，平均数，标准差，最小值和最大值。这是MEANS语句在缺省时得到的基本结果。若需要对数据作更详细的描述，则要指明所需的统计量（关于统计量的概念见课本2.2.4）。用MEANS过程所计算的统计量，在这里也一并列出。下面还会反复提到“变量”这一名词，关于变量的概念

5、见课本2.2.1。MEANS过程所计算的统计量是用关键词表示，这些关键词及其含义如下：N：输入的观测值（观测值的概念见课本2.2.1）个数NMISS：每个变量所含缺失值的个数MEAN：变量的平均数STD：变量的标准差MIN：变量的最小值MAX：变量的最大值RANGE：变量的极差SUM：变量所有值的和VAR：变量的方差USS：每一变量原始数据的平方和（未校正平方和）CSS：每一变量的离均差平方和（校正平方和）CV：变异系数STDERR：每一变量的标准误差（平均数的标准差）T：在H0:μ=0时的t值（见课本5.1.4）PRT：

6、在H0:μ=0的假设下，统计量t大于t临界值绝对值的概率（见课本4.1.1）SKEWNESS：偏斜度KURTOSIS：峭度CLM：置信区间的上限和下限（见课本6.2.1）LCLM：置信区间的下限UCLM：置信区间的上限另外，在PROCMEANS语句中还有12个选项，其中几个主要选项如下：DATA＝（SAS数据集）：指出SAS数据集的名称，若省略，则使用最近产生的数据集MAXDEC＝（数字）：指出所输出的结果中，小数部分的最大位数（0－8），缺省时为8位FW＝（域宽）：指出打印的结果中每个统计量的域宽，缺省时为12VARDE

7、F＝（DF/N）：VARDEF＝DF为缺省值，表示计算方差时，使用n-1作分母，VARDEF＝N表示计算方差时，使用观测值个数n作分母ALPHA＝（α值）：指出在计算置信区间时，选用的显著水平62例2.2计算课本上习题1.2的离均差平方和，方差，偏斜度和平均数的0.95置信区间。解仍然使用2-1data.dat外部数据文件。MEANS过程如下：optionslinesize=76;dataabc;infile‘a:2-1data.dat’;inputx;run;procmeansdata=abcmaxdec=4fw=8a

8、lpha=0.05cssvarskewnessclm;run;运行的结果如下：表2－2：例2.2输出的结果TheSASSystem AnalysisVariable:X CSSVarianceSkewnessKurtosisLower95.0%CLMUpper95.0%CLM2261.969.0842-