retainステートメントに配列指定できる話

もしかしたら、初等技術なのかもですが、最近初めて知った話です。

例えば

があって、
Xで重複を消して、Yを使って、Zを横展開しろと言われたら
迷うことなく

proc transpose data=Q1 out=A1(drop=_NAME_) prefix=Z_;
 var Z;
 by X;
 id Y;
run;

で

ですね。

一番好きなプロシージャはと訊かれたらtransposeと答える僕なので

（訊かれたこと一度もありませんが）、これ以外のやり方がぱっと思いつかないんですよね。

ただ、最近知ったのが

proc sql noprint;

 select max(Y) into:MY

 from Q1;

quit;

data A2;

 set Q1;

 by X;

 array Z_{&MY};

 retain Z_;

 if first.X then call missing(of Z_{*});

  Z_{Y}=Z;

  if last.X then output;

 drop Y Z;

run;

で、この場合結果は同じです。
Yの値が例えば1と4だけなどで飛び値になっている場合、transposeだとZ_1とZ_4だだけが作られますが、こちらの方法だとZ_2とZ_3もnullで作成されるので、ケースによって使い分けれそうです。

何が面白いって、retain 配列名が通るんですね！確かに、通りそうですけど、意外と

その発想がなかったです。

あとついでに関係ないけど、これまた最近知った

data A3;

array AX{2,3,4}(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1);

 D1=dim(AX);

 D2=dim2(AX);

 D3=dim3(AX);

keep D1 D2 D3;

run;

へ～dim関数って　dim2 dim3とかにして、次元ごとの数とれるんですね！

いつも1次元でdim使ってたけど、dim1ってことだったのか～

詰めSAS1回目_全変数全obsで最大の値をとる

詰めSASの1回目のテーマは、データセット中の最大の値を1変数1obsのデータセットに格納する最善手を考えたいと思います。1手？詰めです。

まず、問いの内容は

data Q1;

input X Y Z;

cards;

1 4 3

2 9 8

7 6 5

;

run;

のデータセットから最大の値、この場合、Yの2obs目の9という値を

変数Mのみの新しいデータセットに格納するという目的です。

目的局面図は

です。

以下、解法です。

【解法１】

data A1;

 set Q1 end=eof;

 retain M;

 if max(X,Y,Z)>M then M=max(X,Y,Z);

 if eof;

 keep M;

run;

最大の値をだす場合に、第一感で思い浮かぶのはproc univariateやmeans等の

利用かもしれませんが、その場合は変数ごとの最大をだした後にデータステップで

その中からの最大をだすことになり2ステップになってしまうと思います。

なので、データステップ1回でやっつけたいのですが、SASのmax関数は横（行）方向の最大値

をとる関数なので、そこで出した値を次のobsに持ち込まないといけません。

なのでｒｅｔａｉｎを使って、持ち越し、endオプションで最終obsだけ残します。

【解法2】

proc sql noprint;

 create table A2 as

  select max(M) as M

   from 

  (select X as M from Q1

    union

   select Y as M from Q1

    union

   select Z as M from Q1)

;

quit;

SQL内でmax関数は１変数内を縦にみて最大値を返すので、じゃあ全部の変数を縦に

つないでからかけてやればいいんじゃないかと思って書いてみると

思ったより頭悪い感じのコードになりました。また、サブクエリを使うと、あまり一手といえなくなって

しまいますが。これならサブクエリ内でＸ、Y、Ｚで最大をとってからCASE文で比較して最大を

とる方がスマートだったかも。

【解法3】

data A3;

 obs1=1;

 obs2=2;

 obs3=3;

  set Q1(rename=(X=X1 Y=Y1 Z=Z1)) point=obs1;

  set Q1(rename=(X=X2 Y=Y2 Z=Z2)) point=obs2;

  set Q1(rename=(X=X3 Y=Y3 Z=Z3)) point=obs3;

   M=max(of X1--Z3);

  keep M;

  output;

 stop;

 run;

これは悪ふざけです。横に最大を返すなら、全部横にしてから、かければ

それで決着という方向で書きました。

1回目はここまで。


====================================================
後日、コメントをいただいて追記
====================================================
SQLでの記述について、突っ込みをいれていただきました。
以下のコードで詰みですね！スマート！！
思いつきもしませんでした！

標準SQLのmax関数は引数が1つなのですが、複数の引数を指定した場合
SASはSQLプロシジャ内でもSASのMAX関数と認識するそうです。
奥が深い、、

proc sql;
create table A4 as
select max(max(X,Y,Z)) as M
from Q1;
quit;