全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
3624 5
2007-12-05

从一开头学统计学到现在好几年了,我发现不论是我国现行的教材、期刊都存在大量的统计学错误,呈在cos论坛上探讨过。

我在温习秩和比分析并参照另外两个例题作sas编程的时候发现这两个例题都有不同程度的错误。

例1  来自于《医学综合评价方法及其应用》化学工业出版社2006年2月第一版P44例题4-6,作者在排序分档的表4-14中把观测7错分到第4档,致结果异常。

本人的编程如下:

data w;
input no x1-x5;
cards;
1.00    6.00 .03 .04 54.00 108.00
2.00 7.00 .10 .06 42.00 165.00
3.00 5.00 .26 .48 110.00 180.00
4.00 8.00 .09 .32 85.00 231.00
5.00 4.00 .05 .70 78.00 144.00
6.00 10.00 .19 .62 97.00 240.00
7.00 7.00 .01 .41 39.00 130.00
8.00 8.00 .15 .55 65.00 144.00
;
proc rank descending out=a;var x1-x5 ;ranks r1-r5;/*递降排秩的变量*/
data a1;set a;rs=mean(of r1-r5)/8;/*观测的个数10*/
proc sort data=a1 out=a2; by rs;
proc rank data=a2 out=a3;var rs;ranks rr;
data a4;set a3;sr=lag7(1-1/4/8);/*校正*/rrr=rr/8;if sr=. then r=rrr; if sr^= . then r=sr;
if dif(r)=0 then delete;p=probit(r)+5;if p<= 3.5 then g=1;if p >3.5 & p<=5 then g=2;if p>5 & p<=6.5
then g=3; if p>6.5 then g=4;/*分档情况*/drop rr sr rrr;
proc print data=a4;
proc reg;
model rs=p;
proc glm;
class g;
model rs=g/ss3;
means g/snk;
run;
例2 来自于《中国医院统计》2006年6月13卷2期,在P118的表1诊断符合率、危重病人抢救成功率两列的手工排序上出现严重混乱,致使其一切结果不可信。

编程:略

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2007-12-5 21:15:00
医院都不可信了呀,哎呀,那得死多少人啊
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2007-12-12 22:20:00
hahaha, 你才发现阿?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2008-8-21 09:58:00

ding

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2008-8-28 09:23:00
可以呀 小问题大家没注意的
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2008-8-28 16:30:00
请提高看问题的水平!统计界?民族?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群