作者imaltar (..)
站內Statistics
標題Re: [程式] SAS資料計算
時間Fri Jun 25 21:46:16 2010
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)
1F：→ imaltar:不過日期有重複的話 就可能會有錯 如果有重複日期 06/25 22:18
2F：→ imaltar:那就要再想一下 06/25 22:19
※ 編輯: imaltar 來自: 123.193.72.239 (06/25 22:44)
※ 編輯: imaltar 來自: 123.193.72.239 (06/25 23:11)
3F：→ imaltar:真糟 一直改 不知道原PO有沒有看到>< 06/25 23:11
4F：推 tew:不要用巨集，這個可以簡單處理 06/26 02:43
5F：→ imaltar:更簡單的方式可能是不要用迴圈 我用具集的用意是像 06/26 09:01
6F：→ imaltar:indata我的例子是survival 原po的data名稱也許不一樣 06/26 09:02
7F：→ imaltar:我的time(日期名稱)是令為date 但原po的資料變數名會不同 06/26 09:03
8F：→ imaltar:所以巨集的用意是讓使用者可以根據他原始資料 資料名稱 06/26 09:03
9F：→ imaltar:變數名稱 去彈性輸入 而不用當資料不同時 還要一個個去改 06/26 09:04
10F：→ imaltar:當然能迴圈處理會比較省時 所以我說寫的不夠精簡就是在此 06/26 09:06
11F：→ imaltar:說錯 是迴圈比較費時 06/26 09:07
12F：→ ASD332:巨集跟迴圈... 應該是不一樣的東西吧 @@ 06/26 17:20
13F：→ imaltar:我沒說一樣阿.... 06/26 21:12
14F：→ ASD332:抱歉 眼殘漏看了 06/26 22:05