第六章分类资料的统计描述与简单推断.docVIP

第六章 分类资料的统计描述与简单推断 （医学统计之星） 最后一次更新时间： 在SAS/STAT模块中FREQ、TABULATE和SUMMARY等过程可用于分类资料的统计描述，其中FREQ过程兼具统计描述和统计推断的功能，对分类变量计算频数分布，产生从一维到n维的频数表和列联表；对于二维表，可进行?2检验，对于三维表，可作Mentel-Hanszel分层分析。FREQ过程是SAS用于分析分类资料的一个常用过程。本章将向大家介绍FREQ过程的用法。 §6.1 引 例 例6.1 某医生用国产呋喃硝胺治疗十二指肠溃疡，以甲氰咪胍作对照组，问两种方法治疗效果有无差别（医统第二版P37 例3.10）？ 处 理 愈 合 未愈合 合计 呋喃硝胺 54 8 62 甲氰咪胍 44 20 64 合 计 98 28 126 解：在建立行? 列表的数据集时一般需要三个变量――行变量、列变量和指示每个格子中频数的变量。本例首先建立这样一个数据集，然后调用FREQ过程输出行? 列表，同时利用CHISQ选项进行? 2检验。 ① 设定数据库环境： LIBNAME A C:\USER ； DATA A.YTLI8_1; DATA A.YTLI8_1; 　INPUT ROW COLUMN NUMBER ; DO ROW=1 TO 2； CARDS； DO COLUMN=1 TO 2； 1　1　54 　　　INPUT NUMBER @@; 　1　2　? 8 OUTPUT； 2　1　44 END； 2　2　20 END； ； CARDS； RUN; 　54　8　44　20 ； RUN; FREQ过程，进行? 2检验： PROC FREQ DATA=A.YTLI8_1; 　TABLE ROW*COLUMN / CHISQ; 　WEIGHT NUMBER; RUN; §6.2 FREQ过程 6.2.1 语法格式 PROC FREQ [选项]； TABLES 请求式/[选项]； ? 　WEIGHT 变量名； ? 指定频数指定变量 BY 变量名列； 如：PROC FREQ； TABLES A； 产生变量A的一维频数表。 PROC FREQ； TABLES A*B； 产生A、B两个变量的列联表，*前为行变量，*后为列变量。 在TABLES语句中用*号连接三个或n个变量名，可得到三维或n维列联表，最后一个变量各水平形成表的列，倒数第二个变量各水平形成表的行，其它变量的每一级水平(或水平组合)形成一层，且每一层都形成分离的列联表。 如：PROC FREQ； TABLES A*B*C*D； 列联表的数量为A、B不同水平的组合数，每个表的上部为D的各水平，左侧为C的各水平。 6.2.2 语法说明 【程序步选项】 DATA＝数据集 规定PROC FREQ语句使用的数据集 ORDER＝FREQ 按频数递减顺序排列 ORDER＝DATA 按数据集中出现的顺序排列 ORDER＝INTERNAL 按内部值排列(缺省) ORDER＝FORMATTED 按外部格式值排列 【TABLES语句】 PROC FREQ过程中可有多条TABLES语句，TABLES语句后可接多个表格请求式，每个请求式可包含任何数量的变量，从而得到所需的表格。 如果TABLES语句缺省，则FREQ过程对数据集中的所有变量都给出相应的一维频数表。不规定任何选项时，若需某变量的一维频数，FREQ给出该变量每一水平的频数(freqency)、累积频数(cumulative freqency)、频数的百分比(percent)和累积百分比(cumulative percent)；若需二维频数表，FREQ产生交叉分组列表，即包括各格的频数、总频数的格百分数、行频数的格百分数和列频数的格百分数。 请求式由一个或多个用“*”连接起来的变量名组成。几个变量可放在括号中，如： TABLES A*(B C)；等价于TABLES A*B A*C； TABLES (A-C)*D；等价于TABLES A*D B*D C*D； 下列选项可用于TABLES语句中“/”的后面： 普通选项 OUT＝数据集 建立一个包含变量值和频数计数的输出数据集。如果TABLES语句中不止一个请求式，数据集的内容相应于TABLES语句中最后一个请求。 统计分析选项 CHISQ 对每层作? 2检验，包括Pearson ??2、似然比? 2和Mantel-Haenszel ??2。此外还给出与? 2检验有关的关联指标包括Phi系数、列联系数和Cramer’s V。对于2×2表，给出Fisher精确概率。 AGREE 进行配对? 2检验。 EXACT 对大于2×2的列联表计算Fisher精确概率。同时也给出CHI