sas单变量的特征分析.doc

sas单变量的特征分析 尽管抢不到回家的高铁票，但是我依然还在。大炮，我有个烦恼，我领导最近老叫我单变量结合因变量分析，但是都是分段分析，我总是写proc sql然后group by ,但是这个过程好无聊啊，有木有什么新的代码，让我可以分析的快点啊。最近写了个宏，刚好可以解决你这个问题 在上代码之前，先来个结果图詹大炮这个结果对于分析来说是不好的，因为这个结果没啥实际意义，说白了就是跟因变量没关系，但是这个图我们不是要来讲变量怎么有用，我们要介绍的是这段代码最后呈现的一个结果是怎样的。代码：%macro ChcAnalysis(DSin, DVVar, VarX, NBins, Method, DSChc);proc sort data=&DSin;by &VarX;run;Data temp;set &DSin ;by &VarX;_Obs=_N_;keep &DVVAr &VarX _Obs; run;proc sql noprint;%if &Method=1 %then %do;select count(&DVVar) into :N from temp; select max(_Obs), min(_Obs) into :Vmax, :Vmin from temp; %let BinSize=%sysevalf(&Vmax)/&Nbins);%let LB_1=0;%do i=1 %to %eval(&Nbins-1);%let LB_&i=%sysevalf(&LB_1+(&i-1)*&BinSize);%let UB_&i=%sysevalf(&LB_&i + &BinSize);select sum(&DVVar) , count(*) into :Sum_&i , :N_&i from temp where _obs=&Lb_&i and _obs=LB_&NBins and _obs=&Lb_&i and &VarX=&Lb_&i and &VarX=&Lb_&i and _obs=LB_&i. and _obs=&Lb_&i and &VarX=&Lb_&i and &VarX=&Lb_&i and _obs=&Lb_&i and _obs=LB_&i. and _obs=UB_&i.;%end;这就是第四步啦，%let LB_&NBins=%sysevalf(&LB_1+(&NBins-1)*&BinSize); %let UB_&NBins=&Vmax;产生最后的区间，这里的i是5了，其实我一直很不能理解到这步，为什么还可以输出i=5呢，不是i只循环到4吗？但是执行的时候就是这样子的，这个套路是仿照之前的等高度分变量区间的那个代码写的。我还特地看了日志也循环了：我觉得应该是sas的处理流程，在pdv层面应该可以解释，跪求大神在留言区解释。万分感激。那么&Method=2的部分就留给你自己去看啦，还是&Method=1的那种套路，只是等区间变量等量而已。04data &DSChc;%do i=1 %to &NBins;Bin=&i;LowerBound=&LB_&i;UpperBound=&UB_&i;if (&sum_&i =. ) then N_1=0; else N_1=&Sum_&i;if &N_&i=. then BinTotal=0; else BinTotal=&N_&i;N_0 = BinTotal-N_1;Percent_1=100*N_1/BinTotal;Percent_0=100*N_0/BinTotal;output;%end;Run;然后这部分就是以上产生的宏，拼接成结果数