关于医学等级资料的分析

发现C-kit在肾盂输尿管连接部梗阻段的表达随积水程度的加重呈下降趋势。

[img]file:///D:\Program Files\Tencent\QQ\Users\675904556\Image\9)@F_{UH)[MQ{CQ~P%Y0~%C.jpg[/img]





用什么方法证明以下结论：随着疾病程度的恶化，c-kit的表达呈下降趋势 。
求统计学方法介绍及sas程序。

1. data a;
   
2. do disease=1 to 3;
   
3. do ckt=0 to 3;
   
4. input num @@;
   
5. output;
   
6. end;
   
7. end;
   
8. cards;
   
9.         4        2        3        1
   
10.         6        2        2        1
    
11.    13        1        0        0
    
12. ;
    
13. run;

复制代码

双向有序且属性不同的RxC资料，根据分析目的的不同，采用的统计方法也不同：
1.比较不同积水程度，C-kit表达强度差别是否具有统计学意义，采用秩和检验
2.分析两变量间是否有相关性，采用Spearman秩相关
根据你的分析目的，应该采用Spearman秩相关

1. proc corr spearman;
   
2.         var disease ckt;
   
3.         freq num;
   
4. run;

复制代码

我觉得可以做秩相关分析
spearman  

要是用 等级资料的秩和检验  只能检验出不同组别的水平是否有差别，但是不能看到两个变量的相关性！
可以考虑两变量的 相关性分析   秩相关

蜗牛和我想的一样哈

the other way to look at the question is in the contingency table. we can look at measures, such as GAMMA tau, etc. or CMH. More than that, to appreciate some extra power for ordinal data is through using some typical scores and thus maximizing the test power. jingju

如果其中一类只是两分类的，直接做趋势检验即可。
现在两类多是有序多分类，直接做不出趋势检验。
/*检验相关性*/
proc freq data=a;
   table disease*ckt/cmh;
   weight num;
run;


/*模拟线性趋势*/
proc reg data=a;
   model disease=ckt;
   freq num;
run;

不过，这么做我也没有把握。

webgu 发表于 2013-4-17 15:17
如果其中一类只是两分类的，直接做趋势检验即可。
现在两类多是有序多分类，直接做不出趋势检验。
/*检验 ...

这种列联表形式，可以从独立性（相关性）角度来做，检验行列指标是否独立；
proc freq data=a;
   table disease*ckt/cmh;
   weight num;
run;

由于行列都属于有序资料，所以结果参考Cochran-Mantel-Haenszel第一个统计量。
P值为0.0069 在alpha=0.05水平下，拒绝原假设，认为行列指标相关。

ziyenano 发表于 2013-4-17 16:37
这种列联表形式，可以从独立性（相关性）角度来做，检验行列指标是否独立；
proc freq data=a;
   tabl ...

但是相关能说明趋势吗？

webgu 发表于 2013-4-17 17:24
但是相关能说明趋势吗？

我也想过这个问题，以前看到的资料，在给出相关的结论后，继而给出了趋势的结论。
我觉得能给出趋势的结论的理由：
1.行列均属于有序资料；
2.既然相关，对于有序指标来说，无非正相关，负相关。
有序结合相关，等同于趋势？