科研统计助手
医学科研资讯
前沿热点动态
统计分析服务
当前位置: 菁英统计工作室 > 统计知识 > 软件技巧 > SAS应用 > 正文
多个率比较的SAS编程实现
作者:Admin  更新时间:2012-02-10  来源 圣才学习网 本条信息浏览人次共有

  此例为4*3表资料,对于年龄和病变类型两个属性间的独立性检验,完全可以使用和前述四格表资料分析方法相同的pearson卡方检验,程序的步骤和内容也是基本相同的。
  对此例资料进行分析的程序如下:
libname a 'e:\data\';
data a.case5_2;
 do age=1 to 4;
  do pathtype=1 to 3;
  input count@@;
  output;
  end;
 end;
 datalines;
 60 16 59 208 111 306
 66 79 119 21 47 43
 ;
run;
proc freq;
 tables age*pathtype/chisq;
 weight count;
run;
quit;
  大家可以看到,和上例的程序相比,此例并无多少差异,SAS对于交叉表形式的计数资料,进行卡方检验的程序是完全相同的。其不同的地方是根据我们对资料分析需要的不同而有所区别,比如我们对于线性趋势的检验一般只针对R*C表资料。
  提交以上程序,结果如下。
  The SAS System  20:54 Thursday,November 28.2002 26
  The FREQ Procedure
  Table of age by pathtype
  age  pathtype
  Frequency|
  Percent |
  Row Pct |
  Col Pct |  1|  2|  3| Total
  ---------+--------+--------+--------+
  1 |  60 |  16 |  59 |  135
  |  5.29 |  1.41 |  5.20 | 11.89
  | 44.44 | 11.85 | 43.70 |
  | 16.90 |  6.32 | 11.20 |
  ---------+--------+--------+--------+
  2 |  208 |  111 |  306 |  625
  | 18.33 |  9.78 | 26.96 | 55.07
  | 33.28 | 17.76 | 48.96 |
  | 58.59 | 43.87 | 58.06 |
  ---------+--------+--------+--------+
  3 |  66 |  79 |  119 |  264
  |  5.81 |  6.96 | 10.48 | 23.26
  | 25.00 | 29.92 | 45.08 |
  | 18.59 | 31.23 | 22.58 |
  ---------+--------+--------+--------+
  4 |  21 |  47 |  43 |  111
  |  1.85 |  4.14 |  3.79 |  9.78
  | 18.92 | 42.34 | 38.74 |
  |  5.92 | 18.58 |  8.16 |
  ---------+--------+--------+--------+
  Total  355  253  527  1135
  31.28  22.29  46.43  100.00
  The SAS System  20:54 Thursday,November 28.2002 28
  The FREQ Procedure
  Statistics for Table of age by pathtype
  Statistic  DF  Value  Prob
  Chi-Square  6  58.9049  <.0001
  Likelihood Ratio Chi-Square  6  56.1468  <.0001
  Mantel-Haenszel Chi-Square  1  3.6219  0.0570
  Phi Coefficient  0.2278
  Contingency Coefficient  0.2221
  Cramer's V  0.1611
  Sample Size = 1135
  结果也和例5-1基本相同,有所不同的是此处没有经连续性校正的卡方值(这里不需要校正),也没有自动给出fisher精确概率值。
  在得出两因素不相互独立的结论后,可进一步对因素的各水平进行两两之间的差异检验,这就涉及到R*C表的拆分问题,关于此内容统计学教科书上多数都有介绍,感兴趣的可以自己研究一下。
 
上一篇:[SAS应用]
下一篇:利用SAS/PRO SUMMARY对数据库进行查重

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

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