这段代码的问题到底出在哪呢？？？？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

这是我写的一段代码，但是在我的电脑上运行出来，status一列最后显示的全是healthy，实在是找不出问题出在哪，麻烦哪位给指点一下：

PROC IML;
bmi={15,22,27,19,32,17,23,31,20};
status=J(nrow(bmi),1,"xxxxxxxxx");
do i=1 to nrow(bmi);
if 18.5<=bmi<=20 then
status="healthy";
ELSE status="unhealthy";
end;
print bmi status;
QUIT;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Healthy status Health 我的电脑 BMI status 电脑

这个还真帮不了你。因为根本看不懂

一切重来 发表于 2012-9-25 09:46
这个还真帮不了你。因为根本看不懂

没有关系，我现在还在初步学习中，问题好多

PROC IML;
        bmi={15,22,27,19,32,17,23,31,20};
        status=J(nrow(bmi),1,"xxxxxxxxx");
        do i=1 to nrow(bmi);
                if 18.5<=bmi & bmi<=20 then status="healthy";
                ELSE status="unhealthy";
        end;
        print bmi status;
QUIT;

PROC IML;
        bmi={15,22,27,19,32,17,23,31,20};
        status=J(nrow(bmi),1,"xxxxxxxxx");
        do i=1 to nrow(bmi);
                if bmi[i] >= 18.5 & bmi[i] <= 20 then ;
                        status[i]="healthy";
                ELSE status[i]="unhealthy";
        end;
        print bmi status;
QUIT;

pobel 发表于 2012-9-25 11:40
PROC IML;
        bmi={15,22,27,19,32,17,23,31,20};
        status=J(nrow(bmi),1,"xxxxxxxxx");

嗯嗯，出来结果了，谢谢
可是为什么原来的if 条件不行呢，书上的一段代码就是这么用的，结果却可以运行

henn2010 发表于 2012-9-25 14:48
PROC IML;
        bmi={15,22,27,19,32,17,23,31,20};
        status=J(nrow(bmi),1,"xxxxxxxxx");

嗯嗯，出来结果了，谢谢
可是为什么原来的if 条件不行呢，书上的一段代码就是这么用的，结果却可以运行

jolterheadmmtt 发表于 2012-9-25 18:20
嗯嗯，出来结果了，谢谢
可是为什么原来的if 条件不行呢，书上的一段代码就是这么用的，结果却可以运行

如果是在data步中，a<b<c 的条件就相当于 a<(b<c), 而不是判断 b是否介于a和c之间。
IML中的 a <b<c，既不相当于 a<(b<c), 也不是 a<b & b<c， 可能是 （a<b） <c

data  chk;
    input bmi;
        if 1.1<=bmi<=20 then status="Healthy   ";
        else status="Unhealthy";
        put bmi=  status=;
        cards;
1
2
199
90000
.
;
run;

jolterheadmmtt 发表于 2012-9-25 18:21
嗯嗯，出来结果了，谢谢
可是为什么原来的if 条件不行呢，书上的一段代码就是这么用的，结果却可以运行
...

It is really inefficient to operate on index when one can vectorized an operation. Here is a vectorized one.

PROC IML;
bmi={15,22,27,19,32,17,23,31,20};
status=J(nrow(bmi),1,"unhealthy");
status[loc(bmi>=18.5 & bmi<=20)]="healthy";
print status;
quit;

one can see it is 10 times slower as tested below!!!!!!!!



PROC IML;
*****1)********************;
bmi=j(1e6,1,.);
status=J(nrow(bmi),1,"unhealthy");

call randseed(12345);
call randgen(bmi, 'UNIFORM');           

time=time();
status[loc(bmi>=0.4 & bmi<=0.6)]="healthy";
time=time()-time;
print 'time of vectorized' time;


*****2)********************;
bmi=j(1e6,1,.);
status=J(nrow(bmi),1,"unhealthy");

call randseed(12345);
call randgen(bmi, 'UNIFORM');  


time=time();
do i=1 to nrow(bmi);
                if bmi >= 0.4 & bmi <= 0.6 then ;
                        status="healthy";
                ELSE status="unhealthy";
        end;

time=time()-time;
print 'time of index' time;

quit;

pobel 发表于 2012-9-26 07:46
如果是在data步中，a

原来是这样，谢谢了