题目很有趣 所以我修正一下网友的程式 只要你的资料是符合这样的格式 (重点是没有缺失时间的a9701_data a9703_data 如果像这样有缺失就不适合) 不需要用到巨集 应该都能跑出来 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; data time; set time; order=_n_; run; proc sort data=time;by change_date order; run; proc transpose data=time out=b;by change_date order; run; proc sort data=b;by order descending _NAME_; run; data b; set b; time_0=col1; time_1=lag(col1); time_2=lag2(col1); if order^=lag(order) then time_1=.; if order^=lag2(order) then time_2=.; if substr(left(change_date),1,4)=substr(left(_name_),2,4) then output; keep change_date time_0 -time_2; run; ※ 引述《imaltar (..)》之铭言： : ※ 引述《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: 134.208.29.183