谁会计算ridit分析中的R值？

凤鸾

请教：谁会计算ridit分析中的R值？？
我计算了好久，都是错的～～:Q
先谢谢～～

fendou

Ridit分析程序在SPSS13.0中的实现

文字Ridit分析程序在SPSS13.0中的实现
王孝权 （ 安徽省五河县临北乡石家村卫生室 邮编 233316 ） e-mail;wxq1638@sina.com
因为Ridit分析从实验设计到实施、推断的整个过程比较简单，用途广泛，目前在各种医学杂志出现的频次有增多倾向。但一些著名的统计软件至今没有此分析程序模块，如SPSS13.0及SAS9.0。由于依赖手工计算，在多组，多分类的情况下，带来大量的计算麻烦，乃至准确度难以保证。鉴于此，笔者利用SPSS13.0中的函数与数据管理系统编辑了Ridit分析程序，运行后获得成功，现将其报告。
实       例
  采用《中国卫生统计》2006年第三期280页的文章为依据，进行程序编辑，原始数据整理见表1 （注：本表省略了单元数据的“.00”，实际使用时不可省略可，以免程序不能执行）。
表1 2004年某医院不同科室疗效     
疗效级别             年 份 （1995～2003）               
  j       x1       x2       x3       x4     x5       x6       x7       x8       x 9     
1         183     176   143     155     127       102     132     129       116
2         533     521     757   860     999     1110   1133   1175     1297
3       3428   3689   3855   4850   5562     5842   6389   8069     8610
4       8387   8038   8630   9631   9569     9817 10143   11087   11299
                Ridit程序的实现
  首先把数据按表1形式录入SPSS13.0（或SPSS15.0结果相同），根据例1编辑程序如下：
WEIGHT   By   x9   .
RANK VARIABLES = j / PROPROT   .
COMPUTE y1 = x1 * Pj   .
COMPUTE y2 = x2 * Pj   .
COMPUTE y3 = x3 * Pj   .
COMPUTE y4 = x4 * Pj   .
COMPUTE y5 = x5 * Pj   .
COMPUTE y6 = x6 * Pj   .
COMPUTE y7 = x7 * Pj   .
COMPUTE y8 = x8 * Pj   .
COMPUTE y9 = x9 * Pj   .
FLIP VARIABLES =y1 y2 y3 y4 y5 y6 y7 y8 y9 x1 x2 x3 x4 x5 x6 x7 x8 x9/NEWNAME = j   .
COMPUTE n = SUM(K_1.00,K_2.00,K_3.00,K_4.00)   .
COMPUTE R = LAG(n,9) / n   .
COMPUTE CI.95X = R - 1 / SQRT(3 * n) .
COMPUTE CI.95S = R + 1 / SQRT(3 * n) .
EXECUTE .
程序运行完成后将结果于原文比较见表2；
        原 文   结 果                         程 序 结 果     
年度   R         CI.95X   CI.95S           R         CI.95X   CI.95S
95   .5669     .5617     .5671         .5669       .5617     .5720
96   .5567     .5515     .5619         .5567       .5515     .5619
97   .5539     .5489     .5589         .5532       .5482     .5582
98   .5428     .5382     .5474         .5428       .5382     .5474
99   .5266     .5221     .5311         .5266       .5221     .5312
00   .5229     .5185     .5273         .5229       .5185     .5274
01   .5175     .5132     .5218         .5175       .5132     .5218
02   .5061     .5021     .5101         .5061       .5021     .5102
03   .5000     .4960     .5040         .5000       .4960     .5040
由表2的比较可见，程序运行成功。并检测出在97年有错误（该组的总和原作者多加了50），导致 及其置信区间求解异常，其余置信区间上的微小差异与计算工具、四舍五入有关。
                  讨     论
Ridit分析法是将等级分组的资料转换为一组连续性计量资料的非参数统计方法，在疗效评价，医院综合管理中有较多应用。本次从变量中挑选x9(2003年)作为标准组，也有用各变量频次的和作为标准组，此时，首命令语句是;COMPUTE xi = x1 + x2 + … + xi-1 ,其余不变。不论有多少组数（变量）编辑的Ridit程序，仅表现参与运行变量个数的变化；函数n对应括弧内ki的个数表明每组分类个数。对Ridit法95%置信区间的求解有两种方式，本次采用 R*1 ± 1/SQRT（3n）*2 ［1］，若编程中的置信区间使用R*1 ± 1.96s*3   ［2］后发现，两者结果略有不同。目前，已经有人在SAS［3］、SPSS［4］中编辑了Ridit分析程序。在SAS［3］编程中的置信区间使用后一种方式，但对于同一例题而言，其命令语句的量是本程序的7倍，即使在本次SPSS 的Ridit程序把两种置信区间先后编辑一起出来，语句量仍然不足SAS［3］的1/3 。仔细考察后可以发现，SPSS［4］没有解决置信区间问题，在仅有R值 的情况下，无法做出推断，致使结果失去意义。值得注意的是，有人已经在http://www.acrpnet.cn/viewthread.php?tid=159 上发表了一篇Ridit分析的文章（SAS程序），但是在将其原文复制到SAS9.0运行后，没有得出结果，日志提示了种种错误，且其仍然不能摆脱命令语句海量的弊病。

参考文献
1.   朱本耀.运用Ridit分析法评价及分析医院医疗质量.中国卫生统计.2006，23
（3）：280.
2.   胡良平. 医学统计实用手册.第一版.北京：人民卫生出版社，2004，126.
3.   胡良平. 医学统计实用手册.第一版.北京：人民卫生出版社，2004，208 – 210.
4.   武松，叶冬青.多组单向有序资料的Ridit分析在SPSS中的实现.中国卫生统计.2006，23（6）554.
注意：因为与word文档的不兼容性，有些符号被改变，现解释如下;
*1=平均r值
*2=3倍例数的平方根
*3=总体标准误
可把上述程序复制到SPSS中运行一下，以增进了解。ok!

凤鸾

终于搞定了.

:victory:

凤鸾

我把Ｒidit的具体算法，和在spss中的操作过程过程写出来了。

有兴趣的话，可以来看看偶的博客。
http://blog.businessanalysis.cn/space.php