求高手写代码 - 经管之家

2810

收藏 2011-07-03

求高手写代码，附件是数据集。

要求：FLAG= "Yes" , in the following  cases:

(1) TYPE="Vomiting"  ;

(2) TYPE="Retching",  and its duration >= 5 minutes ;

(3) TYPE="Retching",  and it is preceded or followed by a "Vomiting", and the time  between "Retching" and "Vomiting" <= 1 minute ;

(4) TYPE="Retching",  and it is preceded or followed by other "Retching" events, and  Ending time of the last event - Beginning time of the first event >= 5  minutes, and the time between any two events <= 1 minute ;

我实在是写不出来了，看过别人写的代码，没看懂，怎么看怎么像是错的。由于是真实的数据，所有要求中的（1）（2）两条就把所有的FLAG都置成yes了，无法根据结果判断（3）（4）是否正确。
所以恳求各位高手有时间的话帮我看看。

附件列表

demo.rar

大小:1 KB

马上下载

本附件包括：

demo.sas7bdat

扫码加我拉你入群

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

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

全部回复

yugao1986

2011-7-3 10:15:09

lz可以把看到的代码贴出来吗？

扫码加我拉你入群

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

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

nomad5

2011-7-3 18:24:46

这是整篇代码中关于要求的一部分，单独执行这一部分得出的结果和整篇代码得出的结果不一样，但是只有这一部分是符合要求的。

data demo ;  set demo ;  by id stdt descending endt descending type ;

length flag $3 ;  label flag="Emetic episode? (Yes/No)" ;

length pr_type pr1_type pr2_type pr3_type pr4_type pr5_type $8 ;

retain pr_endt pr_stdt pr_type pr1_endt pr1_stdt pr1_type pr2_endt pr2_stdt pr2_type pr3_endt pr3_stdt pr3_type pr4_endt pr4_stdt pr4_type pr5_endt pr5_stdt pr5_type ;

*format pr_endt pr_stdt pr1_endt pr1_stdt pr2_endt pr2_stdt  pr3_endt pr3_stdt pr4_endt pr4_stdt pr5_endt pr5_stdt datetime20.0 ;

if first.id then do ;
pr_endt=. ; pr_stdt=. ; pr_type="" ; pr1_endt=. ;  pr1_stdt=. ;  pr1_type="" ;
pr2_endt=. ;  pr2_stdt=. ;  pr2_type="" ;  pr3_endt=. ;  pr3_stdt=. ;  pr3_type="" ;
pr4_endt=. ;  pr4_stdt=. ;  pr4_type="" ;  pr5_endt=. ;  pr5_stdt=. ;  pr5_type="" ;
end ;

if stdt>pr1_stdt>. then do ; dif_flag=1 ;  end ;

if type="Vomiting" then flag="Yes" ;
if type="Retching" then do ;
      if endt-stdt>=60*5 then flag="Yes" ;
      else if pr_type="Vomiting" and .<stdt-pr_endt<=60 then flag="Yes" ;
      else if  .<stdt-pr1_endt<=60 then do ;
         if pr1_type="Retching" then do ;
            if endt-pr1_stdt>=60*5 then flag="Yes" ; else do ;
                  if .<pr1_stdt-pr2_endt<=60 then do ;
                     if pr2_type="Retching" then do ;
                        if endt-pr2_stdt>=60*5 then flag="Yes" ; else do ;
                              if .<pr2_stdt-pr3_endt<=60 then do ;
                                 if pr3_type="Retching" then do ;
                                    if endt-pr3_stdt>=60*5 then flag="Yes" ; else do ;
                                          if .<pr3_stdt-pr4_endt<=60 then do ;
                                             if pr4_type="Retching" then do ;
                                                if endt-pr4_stdt>=60*5 then flag="Yes" ; else do ;
                                                      if .<pr4_stdt-pr5_endt<=60 then do ;
                                                         if pr5_type="Retching" then do ;
                                                            if endt-pr5_stdt>=60*5 then flag="Yes" ;
                                                         end ;
                                                      end ;
                                                end ;
                                             end ;
                                          end ;
                                    end ;
                                 end ;
                              end ;
                        end ;
                     end ;
               end ;
            end ;
         end ;
   end ;
end ;

pr_type=type ;  pr_stdt=stdt ;  pr_endt=endt ;

if first.id or dif_flag=1 then do ;
pr5_type=pr4_type ;  pr5_stdt=pr4_stdt ;  pr5_endt=pr4_endt ;
pr4_type=pr3_type ;  pr4_stdt=pr3_stdt ;  pr4_endt=pr3_endt ;
pr3_type=pr2_type ;  pr3_stdt=pr2_stdt ;  pr3_endt=pr2_endt ;
pr2_type=pr1_type ;  pr2_stdt=pr1_stdt ;  pr2_endt=pr1_endt ;
pr1_type=type ; pr1_stdt=stdt ; pr1_endt=endt ;
end ;

/*drop dif_flag pr: ;*/

run ;

扫码加我拉你入群

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

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

honghejing

2011-7-5 11:08:55

这么多end？看得头晕。。。

扫码加我拉你入群

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

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

SAS程序员

2011-7-5 11:58:25

头晕！！！！！

扫码加我拉你入群

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

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

nomad5

2011-7-5 12:01:31

我已经头晕了好多天了

扫码加我拉你入群

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

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

点击查看更多内容…

nomad5

2011-7-5 12:11:51

重新编辑了一下谁帮我看看啊

扫码加我拉你入群

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

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

coolsky2004

2011-7-7 22:41:42

/**************** begin*****************/
求高手写代码，附件是数据集。

要求：FLAG= "Yes" , in the following  cases:

(1) TYPE="Vomiting"  ;

(2) TYPE="Retching",  and its duration >= 5 minutes ;

(3) TYPE="Retching",  and it is preceded or followed by a "Vomiting", and the time  between "Retching" and "Vomiting" <= 1 minute ;
  可以对此进行分解: R and (lagR is V or leadR is V) and duration<=1  等于 R and lagR is V or R and leadR is V and abs(dif(duration))<=1

(4) TYPE="Retching",  and it is preceded or followed by other "Retching" events,
         and  Ending time of the last event - Beginning time of the first event >= 5  minutes,
         and the time between any two events <= 1 minute ;

/**************** end*****************/

/**************** 代码开始*****************/
data aa; set sasuser.demo;duration=(endt-stdt)/60;run;
data bb; set aa;
if  TYPE="Vomiting" then FLAG= "Yes";
if  TYPE="Retching" and duration>=5
         then FLAG= "Yes";
if  ((TYPE="Retching" and lag(TYPE)="Vomiting") or  (TYPE="Vomiting" and lag(TYPE)="Retching")) and       abs(dif(duration))<=1
         then FLAG= "Yes";
if  (TYPE="Retching" and lag(TYPE)="Retching") ) and (endt-lag(stdt))/60>=5 and       abs(dif(duration))<=1
         then FLAG= "Yes";

run;

扫码加我拉你入群

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

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

coolsky2004

2011-7-7 22:54:32

将上面的数据集以excel的形式附上。
ID       STDT       ENDT       TYPE       duration
33798       05APR2010:23:00:00       05APR2010:23:05:00       Retching       5
33798       06APR2010:01:30:00       06APR2010:01:35:00       Retching       5
33798       06APR2010:09:15:00       06APR2010:09:20:00       Retching       5
33798       07APR2010:06:30:00       07APR2010:06:40:00       Retching       10
33798       07APR2010:08:30:00       07APR2010:08:45:00       Vomiting       15
33798       07APR2010:08:30:00       07APR2010:08:45:00       Retching       15
33798       07APR2010:11:00:00       07APR2010:11:05:00       Retching       5
33798       07APR2010:14:30:00       07APR2010:14:35:00       Retching       5
33798       08APR2010:06:30:00       08APR2010:06:35:00       Retching       5
33798       08APR2010:14:30:00       08APR2010:14:40:00       Vomiting       10
33798       08APR2010:14:30:00       08APR2010:14:40:00       Retching       10
33798       09APR2010:07:30:00       09APR2010:07:40:00       Retching       10
33798       09APR2010:16:00:00       09APR2010:16:05:00       Retching       5
34018       05APR2010:13:40:00       05APR2010:13:42:00       Vomiting       2
34018       05APR2010:17:30:00       05APR2010:17:32:00       Vomiting       2
34018       05APR2010:19:01:00       05APR2010:19:03:00       Vomiting       2
34018       06APR2010:14:25:00       06APR2010:14:27:00       Vomiting       2
34018       06APR2010:18:40:00       06APR2010:18:42:00       Vomiting       2
34018       06APR2010:20:05:00       06APR2010:20:07:00       Vomiting       2
34018       07APR2010:10:35:00       07APR2010:10:38:00       Vomiting       3
34018       07APR2010:13:05:00       07APR2010:13:08:00       Vomiting       3
34018       08APR2010:10:30:00       08APR2010:10:32:00       Vomiting       2
65397       07JUL2010:21:00:00       07JUL2010:21:05:00       Vomiting       5
65397       10JUL2010:20:45:00       10JUL2010:20:50:00       Vomiting       5
91184       14AUG2010:22:30:00       14AUG2010:22:35:00       Vomiting       5
142715       10OCT2010:05:00:00       10OCT2010:05:10:00       Retching       10
142715       11OCT2010:05:00:00       11OCT2010:05:40:00       Retching       40
142715       12OCT2010:05:00:00       12OCT2010:05:05:00       Vomiting       5
142715       12OCT2010:10:00:00       12OCT2010:10:30:00       Retching       30
147547       20OCT2010:08:00:00       20OCT2010:08:15:00       Vomiting       15
148906       07NOV2010:03:00:00       07NOV2010:03:10:00       Vomiting       10
148906       07NOV2010:03:00:00       07NOV2010:03:10:00       Retching       10
148906       07NOV2010:09:40:00       07NOV2010:10:00:00       Retching       20
148906       07NOV2010:11:30:00       07NOV2010:11:50:00       Retching       20
149058       23OCT2010:16:30:00       23OCT2010:22:00:00       Vomiting       330

扫码加我拉你入群

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

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

nomad5

2011-7-10 00:35:45

谢谢coolsky2004肯花费时间写
有一个问题，每一个ID是一个实例，所以应该按ID分组，这个怎么加上去。

扫码加我拉你入群

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

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

ntsean

2011-7-10 10:32:43

................................................

扫码加我拉你入群

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

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

ntsean

2011-7-10 10:42:10

第三第四个if有问题吧

coolsky2004 发表于 2011-7-7 22:41
/**************** begin*****************/
求高手写代码，附件是数据集。

要求：FLAG= "Yes" , in the following  cases:

(1) TYPE="Vomiting"  ;

(2) TYPE="Retching",  and its duration >= 5 minutes ;

(3) TYPE="Retching",  and it is preceded or followed by a "Vomiting", and the time  between "Retching" and "Vomiting"

扫码加我拉你入群

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

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

nomad5

2011-7-14 09:58:13

根据别人写的改的，大家点评一下，合理吗？

proc sort data=test.demo out=demo;
by id descending stdt descending endt;
run;

data aa;
set demo;
by id descending stdt descending endt;
type_1=lag1(type); stdt_1=lag1(stdt); endt_1=lag1(endt);
type_2=lag2(type); stdt_2=lag2(stdt); endt_2=lag2(endt);
type_3=lag3(type); stdt_3=lag3(stdt); endt_3=lag3(endt);
type_4=lag4(type); stdt_4=lag4(stdt); endt_4=lag4(endt);
type_5=lag5(type); stdt_5=lag5(stdt); endt_5=lag5(endt);
if first.id then do;
type_1=''; type_2='';type_3='';type_4='';type_5='';
stdt_1=.; stdt_2=.; stdt_3=.; stdt_4=.; stdt_5=.;
endt_1=.; endt_2=.; endt_3=.; endt_4=.; endt_5=.;
end;

if type='Vomiting' then flag='yes';

if type='Retching' then do;
  if endt-stdt>=5*60 then flag='yes';
  if type_1='Vomiting' and stdt-endt_1<=1*60 then flag='yes';
  if type_1='Retching' and endt_1-stdt<=1*60 then do;          if type_2='Vomiting' and endt-stdt_1>=5*60 then flag='yes';
else if type_2='Retching' and endt_2-stdt_1<=1*60 then do;
   if type_3='Vomiting' and endt-stdt_2>=5*60 then flag='yes';
   else if type_3='Retching' and endt_3-stdt_2<=1*60 then do;
      if type_4='Vomiting' and endt-stdt_3>=5*60 then flag='yes';
      else if type_4='Retching' and endt_4-stdt_3<=1*60 then do;
         if type_5='Vomiting' and endt-stdt_4>=5*60 then flag='yes';
      else if type_5='Retching' and endt_5-stdt_4<=1*60 then flag='yes';
      end;
   end;
end;
  end;
end;

run;

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群