注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

公卫人

博学而笃志 切问而近思

 
 
 

日志

 
 

卡方分析 趋势分析 SAS学习  

2013-07-25 10:56:23|  分类: 统计分析 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
很多同学在做趋势卡方检验时感觉很不方便,现讨论如下:
资料:(赵耐青-医学统计学 例15.3)
表 病程与依沙酰胺疗效的关系
病程 疗效 合计
治愈 好转 无效
<1月 24 10 3 37
1~3月 10 4 2 16
3月~5年 35 26 10 71
>5年 10 8 4 22
合计 79 48 19 146
文章原文:资料为双向有序,宜采用趋势卡方检验考查病程与依沙酰胺疗效的相关关第(Nonzero correlation)计算得出卡方值为3.302
1.SPSS做法:其Linear-by-linear Association结果用于线性相关检验,两变量均为等级资料,且均从小到大排列时有意义,其他情况可忽略。
在SPSS内用选Chisq选项其给出如下结果
Chi-Square Tests
Value df Asymp. Sig. (2-sided)
Pearson Chi-Square 3.880Angel 6 .693
Likelihood Ratio 3.944 6 .684
Linear-by-Linear Association 3.302 1 .069
N of Valid Cases 146
a 3 cells (25.0%) have expected count less than 5. The minimum expected count is 2.08.
即卡方值为3.302,P=0.069>0.05,不拒绝0相关
2.SAS做法:
proc freq;
weight count;
tables bingc*liaox/cmh;
run;
在SAS里面非零相关及行平均得分及一般相关结果:
Cochran-Mantel-Haenszel Statistics (Based on Table Scores)
Statistic Alternative Hypothesis DF Value Prob
1 Nonzero Correlation 1 3.3020 0.0692
2 Row Mean Scores Differ 3 3.3621 0.3391
3 General Association 6 3.8537 0.6965
其非零相关结果与SPSS一致。
在SAS里面有个trend命令(在CMH位置),其也可以做趋势检验,其结果给出的是Z值,其中Z值的平方与卡方值相等,而这其中只是两个变量之间的趋势,多个变量我没有运算出来。
但是做两个变量间卡方趋势时,其非零相关值与趋势之间有一点点小差异:比如我做非零相关75.7034,而趋势值是75.79,也就是Z^2,Z值是-8.7057(倪宗瓒-医学统计学例6.8),其非零相关是与SPSS结果一致的
3.STATA做法:
没有找到相就的命令:我用ALL这个命令只给出了
Pearson chi2Devil = 3.8803 Pr = 0.693
likelihood-ratio chi2Devil = 3.9440 Pr = 0.684
Cramér's V = 0.1153
gamma = 0.2205 ASE = 0.115
Kendall's tau-b = 0.1371 ASE = 0.072
这几组的统计量,应该是要编程了。
另外,关于相关与回归是否用到与方差分析方面的原理,希望能有人交流,相就不足之处指正和补充
  评论这张
 
阅读(1788)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017