※ 引述《haoyunwhite (毛毛)》之铭言： : [软体程式类别]: : SAS : [程式问题]: : 资料处理 : [软体熟悉度]: : [问题叙述]: : 变数有异动之时间点(change_date)以及各月份的指标 如下面的例子 : change_date 9701_data 9702_data 9703_data 9704_data 9705_data 9706_data : __________________________________________________________________________ : 970401 1 0 1 0 4 1 : 970204 1 1 2 2 0 3 : 想要利用异动时间点判定异动的月份撷取相对应的变数， : 也就是说，第一笔资料要抓的是9704_data的值; : 第二笔要抓的是9702_data的值。 : 更进一步要做的是取得异动时间点後面两个时间点的值， : 也就是，第一笔资料要抓的是9705_data以及9706_data的值; : 第二笔资料要抓的是9703_data以及9704_data的值; : 资料型态变成 : change_date Time_0 Time_1 Time_2 : __________________________________________ : 970401 0 4 1 : 970204 1 2 2 : 想法 : 目前只有想到土法炼钢的方法 : if substr(change_date,1,4)=9702 then do; : Time_0=9702_data ; Time_1=9703_data ; Time_2=9703_data ; end; : else if substr(change_date,1,4)=9706 then do; : Time_0=9706_data ; Time_1=9707_data ; Time_2=9708_data ; end; : 如果只有一年的资料还好 : 但是如果延伸至五年，可能程式中月份重复的部分会增加很多 : 而且未来不一定会只抓後面2期，可能会前後的数值都要撷取， : 1.想要请问如果撷取到年月(9706)可以转而撷取相对应数值的栏位(9706_data)吗 : 2.是否有其他的作法，得到相同的结果呢 : 谢谢大家有耐心看完 假设一个资料 data time;input change_date a9701_data a9702_data a9703_data a9704_data a9705_data a9706_data a9707_data a9708_data a9709_data a9710_data a9711_data a9712_data; cards; 970401 1 0 1 0 4 1 1 0 2 1 3 2 970204 1 1 2 2 0 3 2 2 0 2 1 2 970105 2 1 0 3 2 1 2 0 3 3 0 0 970311 1 0 2 0 3 1 2 3 1 1 1 2 970912 1 2 3 1 0 2 1 1 3 0 2 1 971002 2 1 3 3 2 3 1 2 0 1 2 2 971101 1 2 3 3 0 0 1 1 2 2 3 1 971225 2 2 1 1 3 3 1 2 1 0 2 3 ; run; 在这边change_date我是直接用数字型的,因为你下面的程式是根据数字型写的, 如果是日期格式的话,那就要修改 /**************************************** indata=原始资料 start=起始日期(起始月份) end=终止日期(终止月份) m1=往後观察的个数(月份) m2=往前观察的个数(月份) out=输出资料 *****************************************/ %macro time(indata=,start=,end=,m1=,m2=,out=); data &out.;set _null_; %do i=&start. %to &end.; data temp;set &indata.;if substr(left(change_date),1,4)=&i.; if substr(left(change_date),1,4)=&i. then do;Time_0=a&i._data ; %do j=1 %to &m1.; %if %substr(%eval(&i.+&j.),3,2) <= 12 %then %do; Time_&j.=a%eval(&i.+&j.)_data ; %end; %else %do; Time_&j.=.; %end; %end; %do k=1 %to &m2.; %if %substr(&i.,3,1)=0 %then %do; %if %substr(&i.,4,1)> &k. %then %do; Time_bef_&k.=a%eval(&i.-&k.)_data ; %end; %else %do; Time_bef_&k.=.; %end; %end; %else %do; %if %substr(&i.,3,2)> &k. %then %do; Time_bef_&k.=a%eval(&i.-&k.)_data ; %end; %else %do; Time_bef_&k.=.; %end; %end; %end; end;run; data &out.;set &out. temp(keep=change_date Time_0 Time_1-Time_&m1. Time_bef_1-Time_bef_&m2.);run; %end; %mend; %time(indata=time,start=9711,end=9712,m1=3,m2=4,out=out); 不过这只能抓一年的资料,如果跨年度,例如97年11月,想向後追踪3个月,只能抓到12月 而已,除非change_date改成日期格式,才能调整%macro成可以跨年度去追踪 --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 114.43.119.2 ※ 编辑: imaltar 来自: 114.43.119.2 (06/21 15:24) ※ 编辑: imaltar 来自: 114.43.119.2 (06/21 15:47)