Блог

Я хочу автоматизировать какой-то процесс. Я использовал SAS proc compare, и я могу получить список переменных с неравными значениями. Далее я хочу извлечь эти переменные (с неравными значениями) и сравнить их средние / медианные / минимальные / максимальные и т. Д., Используя либо proc means, либо proc uni-variate. Мой вопрос в том, как я могу сохранить результат сравнения процесса в виде таблицы и извлечь оттуда переменные? Спасибо.

OUTSTATS= Опция выведет базовую статистику для переменных. Статистика такова N, MEAN, STD, MIN, MAX, STDERR, T, PROBT, NDIF, DIFMEANS, and R,RSQ

Если вам нужно вычислить статистику, отличную от этой, вы можете дополнительно обработать таблицу OUTSTATS, чтобы создать список переменных, которые имели некоторые различия (per NDIF ).

Пример:

 data have1 have2;
  do row = 1 to 100;
    array x(100);

    do _n_ = 1 to dim (x);
      x(_n_) = _n_ * 1000   floor(50*ranuni(123)) - 25;
    end;
    output have1;

    * every 5th row in every 5th column have2 could be different;
    if mod(row,5) = 0 then
      do _n_ = 1 to dim (x);
        if mod(_n_,5) = 0 and ranuni(123) < _n_ / 100 then x(_n_)   _n_;
      end;
    output have2;
  end;
run;

proc compare noprint 
  base=have1 
  compare=have2 
  out=differences
  outstats=summary_stats
  outnoequal
  ;
run;
* review summary_stats;

* need more stats than in summary_stats ?
* get list of variables have some differences;
proc sql;
  reset noprint;
  select _var_
  into :vars_that_differed separated by ' '
  from summary_stats
  where _TYPE_ = 'NDIF' and (_BASE_ ne 0 or _COMP_ ne 0)
  ;
quit;

* show the variables that would used in VAR statement of subsequent MEANS or UNIVARIATE;
%put NOTE: amp;=vars_that_differed;
----- LOG -----
NOTE: VARS_THAT_DIFFERED=x5 x10 x20 x25 x30 x35 x40 x45 x50 x55 x60 x65 x70
x75 x80 x85 x90 x95 x100
  

Если ваш вопрос касается сохранения выходных данных Proc compare в другом наборе данных / таблице, вы можете использовать опцию out:

 proc compare base=old compare=new
out=Out_ds outnoequal outbase outcomp outdif noprint;
id code;
run;
  

out_ds сохранит результаты.

См. Ниже, чтобы сохранить только разные имена переменных в выходном наборе данных:

 data old;
input code A B C D;
datalines;
101 1 a 1 100
102 2 b 1 104
103 3 c 1 54
104 4 d 1 87
105 5 e 1 201
;
run;

data new;
input code A B C D;
datalines;
101 1 a 1 100
102 2 b 1 13
103 3 c 1 54
104 4 d 2 87
105 5 e 1 12
;
run;

proc sort data=old; by code; run;
proc sort data=new; by code; run;

/*suppresses the writing of an observation to the output data set when all values in the observation are judged equal. 
  In addition, in observations containing values for some variables judged equal and others judged unequal, 
  the OUTNOEQUAL option uses the special missing value ".E" to represent differences and percent differences for variables judged equal.*/
proc compare base=old compare=new
out=Out_ds outnoequal;
id code;
run;

/*Get all the variable names from output dataset which you are comparing*/
proc sql ;
  select strip(name)
    into :vnames
    separated by " "
    from dictionary.columns
   where libname="WORK" and
         upcase(memname)="OUT_DS" and
         upcase(strip(name)) not in('_TYPE_','_OBS_','CODE')
  ;
quit;

/*This macro will loop through every value and will keep only those variables in keep_vars, which have unequal values*/
options merror mprint nosymbolgen mlogic;
%macro keepv(varlist);
data out_ds1;
length keep_vars $100.;
     set out_ds;
     retain keep_vars;
     %let var_c=%sysfunc(countw(amp;varlist));
     %do i=1 %to amp;var_c;
     %let var_t=%sysfunc(scan(amp;varlist,amp;i));
        if strip(amp;var_t) ne 'E' and findc(keep_vars,"amp;var_t") ne 1 then 
        do;
            keep_vars=catx(",",keep_vars,"amp;var_t");
        end;
     %end;
    run;
%mend keepv;

%keepv(amp;vnames);

/*keep the last obs -  keep_vars have the required variable names only*/
data out_ds_final;
if 0 then set out_ds1 nobs=nobs end=eof;
set out_ds1 point=nobs;
output;
stop;
keep keep_vars;
run;

proc print data=out_ds_final; run;