分类数据的统计分析及SAS编程课件.ppt

  1. 1、本文档共56页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
分类数据的统计分析及SAS编程课件

一 . 概 述 ;;二 . 2×2 表 ; 2.2.1 Mantel-Haenszel 卡方——Q;如上例江剪刀草组无效(263例)的期望值为: 方差 : 本例 : 样本足够大时 近似正太分布,由此得: 结论为两组疗效不同 。;2.2.2 Pearson 卡方——Qp;计算Pearson卡方先要计算出4个期望值 ,一个 已算出为239.87。其他则可由合计值减出。 Mantel-Haenszel 与Pearson 的关系如下: = 本例 : = ;2.2.3 SAS 程序 [ 程序2-1 ] data bronchit; input treat $ outcome $ count; cards; jiang y 70 jiang n 263 dan y 102 dan n 180 ; proc freq order = data; weight count; tables treat * outcome/chisq; run; ; 检验结果 : STATISTICS FOR TABLE OF TREAT BY OUTCOME 表中各列依次为统计量、自由度(DF)、卡方值及P值。 ; 2.7 配对资料四格表 ;2.7.3 程序 在SAS第6.10版以后,在proc freq 中tables语句中加入了agree选项作McNemar检验。本例程序如下: [ 程序2-2 ] data culture; input a $ b $ count; cards; y y 36 y n 34 n y 0 n n 135 ; proc freq order = data; weight count; tables a * b/agree; run; ;2.7.4 输出 除频数部分外,输出的结果如下: STATISTICS FOR TABLE OF A BY B McNemar’s Test 其中,Kappa系数是反映一致性的指标。 ;三. 多层 2×2 表;3.2 Mantel-Haenszel 检验 ; 3.2.2 符号 对于多层四格表,除四格表中数字以 代表外,以h表示层次,则第h层的四格表可用以下符号表示(表3-2): 表3-2 多层四格表的符号表示 ;3.2.3 QMH 统计量 在无效假设为两组疗效相同时,第h层 的期望值 及方差 的计算公式为: = (3.1) = (3.2) 校正中心因素的两组之间疗效的差异可用Mantel-Haenszel 1959年提出的统计量表示。 = 式中q为层数,本例q=3。 当q层的i行合计数 的合计数 大于30时近似服从自由度等于1的卡方分布。 ;Mantel-Haenszel方法,消除了层次因素的干扰而提高了检出关联的把握度。 当各层次的阳性率的方向一致时,如每层的治愈率都是试验组高于对照组,则QMH检验效果很好。如果各层的方向不一致则可能检验不出其关联。 ;3.2.4 程序 [ 程序3-1] data park2; input center treat $ response $ count @@; cards; 1 test y 12 1 test n 18 1 placebo y 15 1 placebo n 15 2 test y 31 2 test n 9 3 placebo y 34 2 placebo n 6

文档评论(0)

yan698698 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档