请求双向有序且属性不同R*C列联表的线性趋势检验的SAS程序!我的题目如下:谢谢高手指点,
某研究者欲研究年龄与冠状动脉粥样硬化等级之间的关系,将278例尸解资料整理成表7-13,问年龄与冠状动脉粥样硬化等级之间是否存在线性变化趋势?
年龄(岁)X 冠状动脉粥样硬化等级(Y) 合计
— + ++ +++
20- 70 22 4 2 98
30- 27 24 9 3 63
40- 16 23 13 7 59
>=50 9 20 15 14 58
合计 122 89 41 26 278
我查阅了一些文献,程序如下:
data ac;
do y=1 to 4;
do x=1 to 4;
input f @@;
output;
end;
end;
cards;
70 22 4 2
27 24 9 3
16 23 13 7
9 20 15 14
;
run;
proc freq data=ac;
weight f;
tables y*x/chisq;
run;
proc corr data=ac spearman;
freq f;
var y;
with x;
run;
proc reg data=ac;
freq f;
model y=x;
run;
但是在proc reg 中算出来的回归系数b与课本上的不一样,课本上的b=0.4053,标准差为0.0026,sas做出的为b= 0.56457
标准差为 0.06238,这是为什么,有谁有更好的程序,分享一下!谢谢!
还有文献中的例子就是用这个程序来做的,有结果!而为什么我用这个程序,换了数据就没有结果!是不是并非所有的双向有序属性不同的资料都能用SAS作出线性趋势性检验,上述的那个例子是我在《医学统计学》(研究生第二版)教材中找到的