SAS使汇总统计数据在过程平均值中不可用

我有一个有很多列的表，但是为了解释我的 问题我将使用这个简单的表

data test;
    input a b c;
    datalines;
    0 0 0
    1 1 1   
    . 4 2
    ;        
run;

我需要计算常见的汇总统计信息，如最小值、最大值和缺失数量。但我还需要计算一些特殊的数字，如本例中高于某个级别的值的数量>0和>1

我可以使用proc-mean，但它只给我一些正常的结果，比如min，max等

我想要的是以下格式的结果：

var minval maxval nmiss n_above1 n_above2
a    0     1   1     1        0
b    0     4   0     2        1
c    0     2   0     2        1   

我已经能够用这个相当于 愚蠢的代码：

data result;
    set test(keep =b) end=last;
    variable = 'b';
    retain minval maxval;
    if _n_ = 1 then do;
        minval = 1e50;
        maxval = -1e50;
    end;
    if minval > b then minval = b;
    if maxval < b then maxval = b;
    if b=. then nmiss+1;
    if b>0 then n_above1+1;
    if b>2 then n_above2+1;
    if last then do;        
        output;
    end;
    drop b;
run;

我知道必须有更好的方法来做到这一点。我习惯了蟒蛇和熊猫。在这里，我将只循环遍历每个变量，计算不同的摘要statistick，并将结果附加到每个变量的新数据帧中

我可能还可以使用procsql。下一个例子

proc sql;
    create table res as
    select count(case when a > 0 then 1 end) as n_above1_a,
           count(case when b > 0 then 1 end) as n_above1_b,
           count(case when c > 0 then 1 end) as n_above1_c
    from test;
quit;

这给了我：

n_above1_a n_above1_b n_above1_c
1          2          2

---

但这并不能解决我的问题

如果您为每一行添加一个唯一标识符，那么您可以使用PROC TRANSPOSE和PROC SQL来获得结果

data test;
  input a b c;
  id+1;
datalines;
0 0 0
1 1 1   
. 4 2
;        
proc transpose data=test out=tall ;
 by id ;
run;
proc sql noprint ;
  create table want as 
    select _name_
         , min(col1) as minval
         , max(col1) as maxval
         , sum(missing(col1)) as nmiss
         , sum(col1>1) as n_above1
         , sum(col1>2) as n_above2
    from tall
    group by _name_
  ;
quit;

结果

Obs    _NAME_    minval    maxval    nmiss    n_above1    n_above2

 1       a          0         1        1          0           0
 2       b          0         4        0          1           1
 3       c          0         2        0          1           0

---

您熟悉SQL吗？SAS有一个PROC SQL，这是实现这一点的一种方法。如果你想坚持使用PROC-MEANS，你可以使二进制变量大于0和大于2，并使用PROC-MEANS对它们求和。我使用过PROC-sql，但我不知道如何解决我的问题。我只使用它来求和、最大和最小。但我看不出如何计算大于x的列。现在无法测试。你不能选择mina、maxa、summissinga、suma>0吗？不知道。我明天试试。

Obs    _NAME_    minval    maxval    nmiss    n_above1    n_above2

 1       a          0         1        1          0           0
 2       b          0         4        0          1           1
 3       c          0         2        0          1           0