科研统计助手
医学科研资讯
前沿热点动态
统计分析服务
当前位置: 菁英统计工作室 > 统计知识 > 软件技巧 > SAS应用 > 正文
多个率比较的SAS编程实现
作者:Admin  更新时间:2012-02-10  来源 圣才学习网 本条信息浏览人次共有
  对于freq过程大家已经有了一个大概的了解,下面我们用实例来演示freq过程在各种情况下处理相应资料的操作方法。
  二、freq过程对各种类型计数资料处理操作方法
  1.两个率比较的卡方检验
  例5-1 一项治疗糖尿病的临床试验中,将268名患者随机分为两组,一组204人服用盐酸苯乙双胍,死于心血管病的有26人,另一组64人服用安慰剂,死于心血管病的有2人,请分析两种疗法的心血管病病死率有无差异?
  此例为典型的四格表资料,是我们最乐意用卡方检验进行分析的数据形式,非常简单的几条语句就可完成对数据的分析。
  编制程序如下:  
libname a 'e:\data\';
data a.case5_1;
 input treat outcome number;
 datalines;
 1 1 178
 1 2 26
 2 1 62
 2 2 2
 ;
proc freq;
 weight number;
 tables treat*outcome/chisq;
run;quit;
  此处treat代表不同的治疗药物,1表示服用盐酸苯乙双胍,2表示服用安慰剂;outcome代表患者的结局(死亡与否),1表示生存,2表示死于心血管病,number则表示患者的人数。Tables语句后的chisq选项指示SAS进行卡方统计量的计算。
  运行上述程序,结果如下。
  The SAS System  20:29 Wednesday,November 27.2002  5
  The FREQ Procedure
  Table of treat by outcome
  treat  outcome
  Frequency|
  Percent |
  Row Pct |
  Col Pct |  1|  2| Total
  ---------+--------+--------+
  1 |  178 |  26 |  204
  | 66.42 |  9.70 | 76.12
  | 87.25 | 12.75 |
  | 74.17 | 92.86 |
  ---------+--------+--------+
  2 |  62 |  2 |  64
  | 23.13 |  0.75 | 23.88
  | 96.88 |  3.13 |
  | 25.83 |  7.14 |
  ---------+--------+--------+
  Total  240  28  268
  89.55  10.45  100.00
  The SAS System  20:29 Wednesday,November 27.2002  6
  The FREQ Procedure
  Statistics for Table of treat by outcome
  Statistic  DF  Value  Prob
  Chi-Square  1  4.8188  0.0282
  Likelihood Ratio Chi-Square  1  6.0021  0.0143
  Continuity Adj.Chi-Square  1  3.8454  0.0499
  Mantel-Haenszel Chi-Square  1  4.8008  0.0284
  Phi Coefficient  -0.1341
  Contingency Coefficient  0.1329
  Cramer's V  -0.1341
  Fisher's Exact Test
  Cell (1,1)Frequency (F)178
  Left-sided Pr <= F  0.0177
  = F  0.9967
  Table Probability (P)0.0144
  Two-sided Pr <= P  0.0330
  Sample Size = 268
  结果给出的内容很是详细,第一部分是四格表的内容,每一个格子内给出四项内容,分别为频数、百分比、行百分比以及列百分比。第二部分为卡方检验的结果,SAS在默认状态下即给出此若干项统计量的结果以及相应的概率,具体的含义大家从字面上就可看出来,这里我不再多说。最后,SAS还给出fisher精确概率的计算结果,提供给我们更多的选择。
  对于病例对照研究,可在tables语句后添加选项‘cmh’,SAS即会给出“Cochran-Mantel-haenszel”统计量的结果,以及相对危险度、比数比等重要指标。如果资料为1:1配对资料,即四格表中数据为不同组合的对子数,则需要采用McNemar检验,此时只需在tables语句后添加选项‘agree’,或添加语句——“exact mcnem;”即可。
  2.R*C表资料的卡方检验
  例5-2 对1135例绝经后出血的妇女进行病理分析,结果见表5.3,请分析病变类型是否与年龄有关。
  表5.3 不同年龄妇女绝经后出血的病变类型

年龄组/岁
病变类型,例数(百分比)
合计
功能性
恶性
良性
≤50
60(44.4)
16(11.9)
59(43.7)
135
51~
208(33.3)
111(17.8)
306(49.0)
625
61~
66(25.0)
79(29.2)
119(45.1)
264
71~
21(18.9)
47(42.3)
43(38.7)
111
合计
355(31.3)
253(22.3)
527(46.4)
1135
上一篇:[SAS应用]
下一篇:利用SAS/PRO SUMMARY对数据库进行查重

本站所发表的文章,大部分严格筛选,来源于各相关专业论坛或专业网站,内容仅供大家学习和参阅尊重原作者版权,勿用于商业用途,转载请注明来源。如有学术争议的文章,或可能与事实不符的,与本站立场无关。如有影响到您权益的文章,请及时通知本站,本站立即删除。谢谢监督。
推荐阅读
 
热门文章

Copyright 2010 Powered By 菁英统计工作室 www.tjstat.com 地址1:武汉市珞瑜路1037号 430074 ; 地址2:武汉市航空路13号 430030
邮箱:tjstat@126.com QQ: 点击这里给我发消息 点击这里给我发消息
鄂ICP备10020011号 |