16统计量计算

频数过程
频数过程（FREQ）用于计算各种形式的频数及 一些检验统计量。
频数过程句法
PROC FREQ options; OUTPUT <OUT= SAS-data-set><output-statistic-list>; TABLES requests / options; WEIGHT variable; EXACT statistic-keywords; BY variable-list;
例16.13 创建包含卡方统计量的数据集。
options nodate pageno=1 pagesize=60; proc freq data=ResDat.color order=data; weight count; tables eyes*hair /chisq expected cellchi2 norow nocol; output out=chisqdat pchi lrchi n nmiss; title 'Chi-Square Tests for 3 by 5 Table of Eye and Hair Color'; run;
表格没有轮廓线和分隔线的FORMCHAR=选项： FORMCHAR(1,2,7)= ' ' ; /* 三个空格 */
例16.6 按格式化值的顺序排列。 proc format; value $sfmt 'M' = 'male ' 'F ' ='female'; proc freq data=ResDat.class order=formatted; table sex; format sex $sfmt.; run;
语句说明：
BY EXACT OUTPUT TABLES TEST WEIGHT 对BY变量定义的观测组分别计算其相应的频数或相等交叉制表 对特定统计量作精确检验 产生包含特定统计量的数据集 产生多变量交叉表并对关联度进行度量和检验 要求对关联度和一致性度量进行近似检验 规定一个变量，其值为每一观测的权数
PROC CORR <option(s)>; BY <DESCENDING> variable-1<...<DESCENDING> variable-n> <NOTSORTED>; FREQ frequency-variable; PARTIAL variable(s); VAR variable(s); WEIGHT weight-variable; WITH variable(s);
其它语句
VAR语句 VAR variable-list; 列出要计算相关系数的变量。 WITH语句 WITH variable-list; 该语句和VAR语句联合使用计算变量间特殊组合的相关系数。用 VAR语句列出的变量放在相关阵的上方，而用WITH语句列出的 变量放在相关阵左边。 PARTIAL语句 PARTIAL variable-list; 计算Pearson偏相关，Spearman偏秩序相关，或Kendall偏tau-b。该 语句给出偏相关变量的名子。
例中，产生区间[0, 1]上均匀分 布的随机数1000个，分别将区间 [0, 1]均分成3和4个小区间，并把 产生的随机数按所属区间转换为整 数。然后对这些整数作频数分析。

例16.10 对One-Way频率表作卡方检验。
proc sort data=ResDat.color; by region; run; proc freq data=ResDat.color order=data; weight count; tables hair/nocum testp=(30 12 30 25 3); by region; title 'Hair Color of European Children'; run;
PROC FREQ语句
PROC FREQ options;
选项说明：
Data= Compress Formchar= Noprint Order= Page 规定输入数据集 在下一个单向频数表不适合页面的空间时强迫在当前页输出 规定用来构造列联表单元的轮廓线和分隔线的字符 规定不输出任何描述统计量 规定输出频数表时分类变量的排序方式 规定每页只输出一张表，否则按每页行数允许的空间输出多张表
proc print data=chisqdat noobs; title 'Chi-Square Statistics for Eye and Hair Color'; title2 'Output Data Set from the FREQ Procedure'; run;
第16章 统计量计算
清华大学经管学院 朱世武 Zhushw@ Resdat样本数据： SAS论坛：
本章将介绍的统计量计算过程包括：
相关过程； 频数过程； 均值过程； 单变量过程。
相关过程
相关过程（CORR）用于计算变量间的相关系数。 相关过程句法
BY语句 BY variable-list; 对BY变量定义的观测组分别计算其相应的简单统计量。当使用BY 语句时，要求输入数据集已按BY变量排序的次序排列，除非指定 NOTSORTED。
应用举例
例16.1 计算Pearson相关系数及其它关联测度。 proc corr data=ResDat.fitness pearson spearman hoeffding; var weight oxygen runtime; title 'Measures of Association for'; title2 'a Physical Fitness Study'; run;
例中，计算数据集ResDat.CLASS中变量SEX的分布，并以格 式化值的顺序排列。
FREQ 过程 累积 累积 Sex 频数 百分比 频数 百分比 --------------------------------------------------female 9 47.37 9 47.37 male 10 63 19 100.00
例16.3 计算两个数据集中相同变量之间的相关系数。 data a; /*数据集准备 */ merge ResDat.Idx000001(keep=date oppr clpr) ResDat.szcz(keep=date oppr clpr rename=(oppr=oppr_sz clpr=clpr_sz) ); by date; run; proc corr data=a nomiss cov; var oppr_sz clpr_sz; with oppr clpr; title2 '长方形的COV和CORR阵'; run; proc corr data=a cov csscp outp=oup; title2 '从含有缺失值的数据集中计算CSSCP和COV'; run; 例中，对上证指数ResDat.Idx000001和深证成指ResDat.SZCZ中的变 量开盘价和收盘价作相应的计算。
WEIGHT语句
WEIGHT variable; 该语句规定一个WEIGHT变量，它的值表示相应 观测的权数。该变量的值应大于零。若这个值<0 或缺失，假定该值为0。
BY语句
BY variable-list; 对BY变量定义的观测组分别计算其相应的简单统计 量。当使用BY语句时，要求输入数据集已按BY变 量排序的次序排列，除非指定NOTSORTED。

应用举例
例16.8 随机数频数分析。 data a; do I=1 to 1000; X=int(uniform(8888)*3)+1; Y=int(uniform(8888)*4)+1; output; end; proc freq data=a(drop=i); title '没有TABLES语句'; run; title; proc freq; tables x x*y/chisq; run; proc freq; tables x*y/list; run;
语句说明：
BY FREQ PARTIAL VAR WEIGHT WITH 分别对每一BY组计算相关系数 规定一个数值变量, 其值为每一观测值出现的频数 给出Pearson, Spearman或Kendall偏相关系数的变量名 给出要计算相关系数矩阵的变量和顺序 计算加权的乘积矩相关系数时给出权数变量名字 计算变量组合之间的相关系数
OUTPUT语句
OUTPUT <OUT= SAS-data-set><output-statistic-list>; 该语句创建一个由PROC FREQ过程输出统计量的SAS数据 集。OUTPUT创建的数集可以包括由TABLES语句规定的任 意统计量。 PROC FREQ过程每一次只允许使用一个OUTPUT语句。当 规定多个TABLES语句时，OUTPUT语句创建的数据集内容 对应于最后那个TABLES语句，当一个TABLES语句中规定多 个表时，OUTPUT创建的数据集内容对应于最后那个表。 选项说明： OUT=规定输出数据集； output-statistic-list规定输出数据集中所包含的统计量。
PROC CORR语句
PROC CORR <option(s)>; PROC CORR语句选项<option(s)>说明由下页表给出。
选项说明：
ALPHA COV DATA= NOPRINT OUTP= OUTS= PEARSON 输出Cronbach系数 输出协方差 输入数据集名 禁止打印输出 规定创建存放Pearson相关系数的数据集 规定创建存放Spearman相关系数的数据集 输出Pearson相关系数
ORDER=选项及说明：