源数据如下表(从2007年1月4日到2007年6月x日，还有其他代码的数据未列出），现想依据3个月的日数据来循环统计lnri的均值、分位数、偏度等统计量，编了错误的code如下：%macro mad;
%do i=1 %to 4;
   if mm>=&i. and mm<=2+&i.;
   proc univariate data=data.sj_9 mean skewness kurtosis p1 noprint;
   var r;
   by code;
   output out=data.sjcl1 mean=mean skewness=skewness kurtosis=kurtosis p1=p1;
run;
%end;
%mend mad;
%mad;
望指正

code	date	month	rm	lnri
000001	2007/1/4	1	-0.0016	-0.025194
000001	2007/1/5	1	0.0272	-0.07733
000001	2007/1/8	1	0.04	0.0234593
000001	2007/1/9	1	0.0215	0.0316552
000001	2007/1/10	1	0.0273	0.0200868
000001	2007/1/11	1	0.009	0.0485196
000001	2007/1/12	1	-0.0336	0.0488546
000001	2007/1/15	1	0.0505	0.0490356
000001	2007/1/16	1	0.0268	0.043224
000001	2007/1/17	1	-0.0061	0
000001	2007/1/18	1	0.0243	-0.050613