医用sas统计分析（三）

《医用sas统计分析（三）》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、医用SAS统计分析第三讲9/9/20211一、分类变量的统计推断—卡方检验在SAS系统中，对分类变量资料的基本统计分析方法主要通过FREQ过程实现的。FREQ过程的主要功能有：1.产生一维或多维频数表；2.计算各种表中格子的理论频数、构成比和各种率；3.对分类变量资料作相应的假设检验。9/9/20212FREQ过程的语句及说明主要格式：PROCFREQ[选择项]；TABLES表达式/[选择项]；WEIGHT<变量名>；说明：PROC过程[选择项]data=数据集；规定PROCFREQ语句使用的数据集。formchar(1,2,7)=‘

2、-+’;规定用来

3、构造列联表单元的轮廓线和分隔线的字符（只有三个字符）。1为垂线，2为水平线，7为水平与垂直的交叉线。9/9/20213TABLES语句PROCFREQ过程中可有多条TABLES语句，TABLES语句后可接多个表格表达式，每个表达式可包含任何数量的变量，从而得到所需的表格。如果TABLES语句缺省，则FREQ过程对数据集中的所有变量都给出相应的一维频数表。不规定任何选项时，若需某变量的一维频数，FREQ给出该变量每一水平的频数(freqency)、累积频数(cumulativefreqency)、频数的百分比(percent)和累积百分比(cumulat

4、ivepercent)；若需二维频数表，FREQ产生交叉分组列表，即包括各格的频数、总频数的格百分数、行频数的格百分数和列频数的格百分数。9/9/20214TABLES语句的表达式表达式是要求FREQ过程分析处理的一维或多维表的清单。一维表有一个变量名表示，二维表由星号“*”联接两个变量名表示，如a*b表示变量a与b的二维表。三维表的形式为a*b*c。TABLES语句的[选择项]1.普通选项out＝数据集建立一个包含变量值和频数计数的输出数据集。如果TABLES语句中不止一个表达式，数据集的内容相应于TABLES语句中最后一个表达式的表格。9/9/20

5、2152.统计分析主要选项chisq对每层作χ2检验，包括Pearsonχ2、似然比χ2和Mantel-Haenszelχ2。此外还给出与χ2检验有关的关联指标包括Phi系数、列联系数和Cramer’sV。对于2×2表，给出Fisher精确概率。agree进行配对χ2检验(McNemar’s检验);一致性检验的Kappa值。exact对大于2×2的列联表计算Fisher精确概率。同时也给出CHISQ选项的全部统计量。cmh给出Cochran-Mantel-Haenszel统计量，主要用于行×列表的统计分析。对于分层2×2表，cmh过程给出总体相对危险度

6、估计及其可信区间，还给出各层关联度指标是否齐性的Breslow检验。9/9/20216measures对每层的二维表计算一系列关联指标及相应的标准误，包括Pearson和Spearman相关系数，以及Gamma和Kendall系数等。对于2×2表，还给出常用的危险度指标及其可信区间。all给出chisq，measures，cmh所请求的全部统计量。Alpha=给出α检验水准。缺省为0.05.expected给出期望频数。9/9/202173.禁止输出选项nofreq不给出列联表中的格频数nopct不给出列联表中的格百分数norow不给出列联表中各格的行

7、百分数nocol不给出列联表中各格的列百分数nocum不给出频数表的累积频数和累积百分数noprint不给出表格，但给出CHISQ、CMH等语句所指定的统计量。WEIGHT语句指明该变量为频数。只能使用一个WEIGHT语句，且该语句作用于所有的表。9/9/20218四格表资料的卡方检验书中例10-1datachisq1;dor=1to2;doc=1to2;inputf@@;output;end;end;cards;63173168;procfreqFormchar(1,2,7)=‘

8、-+’;weightf;tablesr*c/chisq;run;9/9

9、/20219卡方检验结果TABLEOFRBYCRCFrequency

10、二维表每个格子的频数Percent

11、每个格子的频数在总频数中的百分比。RowPct

12、行百分数，每格子频数占该行合计频数的百分比。ColPct

13、*1

14、2

15、Total---------+--------+--------+1

16、63

17、17

18、80

19、35.20

20、9.50

21、44.69

22、78.75

23、21.25

24、*列百分数，每格子频数占

25、67.02

26、20.00

27、该列合计频数的百分比。---------+--------+--------+2

28、31

29、68

30、99

31、17.32

32、37.99

33、55.31

34、31

35、.31

36、68.69

37、

38、32.98

39、80.00

40、---------+--------+-----