Re: [程式] SAS資料計算

[imaltar (..)]

data survial;input date date10. response;
format  date date10. ;
cards;
09AUG07        1
23AUG07        1
06SEP07        2
20SEP07        2
05OCT07        1
16OCT07        1
30OCT07        2
29NOV07        1
01JAN08        2
10JAN08        2
19JAN08        2
30JAN08        2
10JUN08        1
;
run;

在巨集中,indate是放原始資料,本例就是survial
time就是日期變數,本例是date
var就是判斷的變數,本例是response
value2是欲觀察練續情形的值,本例是2
value1就是另外的值,本例是1
out就是最終的資料,本例是survial5
我除了試你的例子,還有再多加一些變化
目前是都可以正確跑出結果,你可以試試看你的data執行有沒有錯
應該還有更有效率的作法,只是目前還沒想到~~~

%macro search(indata=,time=,var=,value1=,value2=,out=);
proc sort data=&indata. out=survial1;by descending &time.;run;
data survial1;set survial1;r_right=lag(&var.);
proc sort data=survial1;by &time.;run;
data survial2;set &indata.;r_left=lag(&var.);
data  survial2;merge survial1 survial2;run;
data  survial3;set  survial2(where=(&var.=&value2.));
data  survial4;set survial3(where=( r_right ne  r_left));
if r_right=. then  r_right=&value1.;if r_left=. then r_left=&value2.;n+1; run;
proc sql noprint;select max(n) into:m_n from survial4;
data _null_;set  survial4(where=(n=1));
if r_right=&value1. then call symput('start',2);
if r_right=&value2. then call symput('start',1);run;
data s;set _null_;
%do i=&start. %to &m_n. %by 2;
data s1_&i.;set survial4;if n=&i.;
data s2_&i.;set survial4;if n=%eval(&i.+1);
data s&i.;set  s1_&i. s2_&i.;
data _null_;set s&i.(where=(r_left=&value1.));call symput('date_start',&time.);
data _null_;set s&i.(where=(r_left=&value2.));call symput('date_end',&time.);
data s&i.;set s&i.;days=intck('day',&date_start.,&date_end.)+1;call
symput('days',days);
data s&i.;merge survial3 s&i.;by &time.;if r_left=&value1. &
r_right=&value1. then delete;
if &date_start. <= &time.<= &date_end.;if days=. then days=&days.; run;
data s;set s s&i.;run;
%end;
proc sort data=s;by &time.;
data &out.;merge &indata.(in=a) s(keep=&time. days);by &time.;if a;run;
%mend;
%search(indata=survial,time=date,var=response,value1=1,value2=2,out=survial5);



--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 123.193.72.239
※ 編輯: imaltar         來自: 123.193.72.239       (06/25 21:48)
※ 編輯: imaltar         來自: 123.193.72.239       (06/25 22:07)
※ 編輯: imaltar         來自: 123.193.72.239       (06/25 22:17)
→ imaltar:不過日期有重複的話 就可能會有錯 如果有重複日期
→ imaltar:那就要再想一下
※ 編輯: imaltar         來自: 123.193.72.239       (06/25 22:44)
※ 編輯: imaltar         來自: 123.193.72.239       (06/25 23:11)
→ imaltar:真糟 一直改  不知道原PO有沒有看到><
推 tew:不要用巨集，這個可以簡單處理
→ imaltar:更簡單的方式可能是不要用迴圈 我用具集的用意是像
→ imaltar:indata我的例子是survival 原po的data名稱也許不一樣
→ imaltar:我的time(日期名稱)是令為date 但原po的資料變數名會不同
→ imaltar:所以巨集的用意是讓使用者可以根據他原始資料 資料名稱
→ imaltar:變數名稱 去彈性輸入 而不用當資料不同時 還要一個個去改
→ imaltar:當然能迴圈處理會比較省時 所以我說寫的不夠精簡就是在此
→ imaltar:說錯  是迴圈比較費時
→ ASD332:巨集跟迴圈... 應該是不一樣的東西吧 @@
→ imaltar:我沒說一樣阿....
→ ASD332:抱歉 眼殘漏看了