サブセット化IFでありがちな落とし穴

SASプログラマ歴が長いと、みんなこれ経験してるんじゃなかろうか。

「サブセット化IF」自体は以下記事で解説しています。

【サブセット化IF】条件に一致するオブザベーションのみ処理を継続する。

😕失敗例1: 「_N_」と組み合わせて失敗しがち

各オブザベーションに1～連番をふった変数を作りたくて、以下のように書いたのですが、、

data dt1;  set sashelp.class;  if age=12;  x = _n_;  keep age x; run;

変数Xに1～連番を入れたいのに、なんかおかしい。

解説

• まずデータステップは内部の動きとして、

オブザベーションの数だけ、データステップをぐるぐるとループ（反復）させています。

ループのたびに1オブザベーションずつ読み込んでいるイメージで、「_N_」はそのループ毎に「+1」した番号が入ります。

• 次に「_N_」は内部の動きとして、

まずSET文・WHERE文で読み込むオブザベーションに対して、データステップがループ（反復）されますが、自動変数「_N_」は、各反復内でいち早く番号がふられます。

番号をふった後で、サブセット化IFにより処理が継続されずに出力もされなかったオブザベーションがあると、番号が連番にならなくなります。

📝特に「IF _N_=1 THEN なんかの処理;」みたいな書き方で失敗しがち。

例えば、以下プログラム（中身はめっちゃ意味のないことやってるのであんまり見ないでください）

data dt2;  set sashelp.class;  if age=12;  if _n_=1 then do;   dcl hash hs(dataset:"sashelp.class");   hs.definekey("name");   hs.definedone();  end;  if hs.check()=0 then x2=1; run;

_N_=1 のオブザベーションがたまたまサブセット化IFによって「処理継続・出力」の対象外になっているので、「IF _N_=1 THEN DO; ～END;」の中の処理も動いていません。

上の例ではエラーが出て分かりやすいですが、書き方によってはエラーが出ない場合もあって失敗に気づかない可能性もあるので要注意です。

😕失敗例2: 「END=オプション」と組み合わせて失敗しがち

最後のオブザベーションにフラグを立てたくて、以下のように書いたのですが、、

data dt3;  set sashelp.class end=_eof;  if age=12;  if _eof then y=1;  keep age y; run;

フラグ立ってないじゃないか！

これは「END=オプション」の挙動を勘違いしていると起こりやすいです。

「END=オプション」はSET文・WHERE文で読み込まれる最後のオブザベーションに作用しますが、その最後のオブザベーションがサブセット化IFによって「処理継続・出力」の対象外になっているためです。

あと思いつくのが、以下で紹介している「FIRST.BY変数」「LAST.BY変数」も、サブセット化IFと組み合わせると、同様の原理で意図しない結果になりやすいですね。

「FIRST.BY変数」と「LAST.BY変数」で、グループ毎の最初と最後のオブザベーションを特定する。

SASで年賀状2023

 

SASで年賀状をつくりました（ほぼ「SASで年賀状」の焼き直し）
2023年もどうぞよろしくお願いいたします！

年賀状をつくるプログラム (SAS9.4M7で動作確認)

*** 描画用のパーツ ; data rabbit;   x1= 0;     y1= 0;   x2= 0;     y2=-6.5;   x3=-0.31;  y3= 1.2;   x4= 0.31;  y4= 1.2;   x5= 0;     y5= 0.15;   x6= 0;     y6= 0.7;   x7=-0.5;   y7= 5;   x8= 0.5;   y8= 5;   x9=-0.5;   y9= 5.5;   x10=0.5;   y10=5.5;   x11=-0.8;  y11=0.1;   x12=0.8;   y12=0.1; run; *** 描画 ; title ; ods graphics on / height=10cm width=14.8cm; proc sgplot data=rabbit noautolegend;  styleattrs wallcolor=khaki;      /* Text */  inset "HAPPY NEW YEAR"  /  position=top textattrs=(color=brown size=31cm);  inset "20"  /  position=bottomleft textattrs=(color=brown size=140cm );  inset "23"  /  position=bottomright textattrs=(color=brown size=140cm );      /* rabbit */  scatter x=x1 y=y1  /  markerattrs=(symbol=circlefilled  size=3cm color=white);  scatter x=x2 y=y2  /  markerattrs=(symbol=circlefilled  size=4cm color=white);  scatter x=x3 y=Y3  /  markerattrs=(symbol=circlefilled  size=0.2cm color=black);  scatter x=x4 y=Y4  /  markerattrs=(symbol=circlefilled  size=0.2cm color=black);  scatter x=x5 y=y5  /  markerattrs=(symbol=ArrowDown  size=0.5cm color=black);  scatter x=x6 y=y6  /  markerattrs=(symbol=TriangleDownFilled  size=0.2cm color=black);  scatter x=x7 y=y7  /  markerattrs=(symbol=CircleFilled  size=1cm color=white);  scatter x=x8 y=y8  /  markerattrs=(symbol=CircleFilled  size=1cm color=white);  scatter x=x9 y=y9  /  markerattrs=(symbol=CircleFilled  size=1cm color=white);  scatter x=x10 y=y10  /  markerattrs=(symbol=CircleFilled  size=1cm color=white);  scatter x=x7 y=y7  /  markerattrs=(symbol=CircleFilled  size=0.7cm color=mistyrose);  scatter x=x8 y=y8  /  markerattrs=(symbol=CircleFilled  size=0.7cm color=mistyrose);  scatter x=x9 y=y9  /  markerattrs=(symbol=CircleFilled  size=0.7cm color=mistyrose);  scatter x=x10 y=y10  /  markerattrs=(symbol=CircleFilled  size=0.7cm color=mistyrose);  scatter x=x11 y=y11  /  markerattrs=(symbol=CircleFilled  size=0.4cm color=pink);  scatter x=x12 y=y12  /  markerattrs=(symbol=CircleFilled  size=0.4cm color=pink);  xaxis min=-3 max=3  display=none;  yaxis min=-5 max=1  display=none; run;

SQLプロシジャのINTO句の注意点

SQLのINTO句を使用する際の注意点を2つ紹介したいと思います。

 

1. 変数にフォーマットが割り当てられている場合

data DT1;   format X yymmdd10.;   X=1; run; proc sql noprint;   select X into :X trimmed   from DT1; quit; %put &X; ログ 1960-01-02

マクロ変数にはフォーマットをあてた値が格納されるようです。この挙動、SASのリファレンスに書いてない気がする。。

もし、フォーマット値ではなく、変数値をマクロ変数に格納したい場合、

事前にフォーマットを除いておくか、もしくは数値変数ならデータに応じた適切な長さのBESTフォーマットを「FORMAT=」で便宜的にあてるなどの対応が必要になります。

proc sql noprint;   select X format=best16. into :X trimmed   from DT1; quit; %put &X; ログ 1

2. 数値変数かつフォーマットが割り当てられていない場合

data DT2;   X=123.456789; run; proc sql noprint;   select X into :X trimmed   from DT2; quit; %put &X; ログ 123.4568

デフォルトの挙動として「BEST8.」の出力形式で格納されるようです。

状況によっては、もとの数値がだいぶ削られちゃいますね（これはINTO句に限ったことではないですが。。）

【PROC FORMAT】CNTLIN=オプションの落とし穴

まず、前提知識として「PROC FORMAT入門9 : CNTLIN=オプション 」をご参照下さい。

そして、ここから本題。以下の例をご覧ください。

失敗例

①まず「データセット」から「フォーマット」を生成

data FMT; input FMTNAME:$10. START:$10. LABEL:$20.; cards; $TEST_ OTHER その他 ; proc format cntlin=FMT; run;

②実データに先ほどのフォーマットを当ててみる。

data TEST; input X:$20.; cards; OTHER . ; data TEST2;  set TEST;  format X $TEST_.; run;

あれ？！

フォーマット当てた結果、2行目の欠損値に「その他」が当てられてしまってる？！

解説

「CNTLIN=」に使用するフォーマット生成用のデータセットに注目

実は上の "OTHER" が、範囲指定としての「OTHER」として解釈されてしまってます（範囲指定？という方は「PROC FORMAT入門3 : 範囲の指定」を参照）

要は、以下のようにFORMATプロシジャで①みたいな定義をしたつもりが、②の定義をしてしまった、ということ。

proc format; /* ① */  value $TEST_  "OTHER" = "その他"  ; run; proc format; /* ② */  value $TEST_  OTHER = "その他"  ; run;

余計なことをしてますね。。

ちなみに、"OTHER" 以外にも、FORMATプロシジャの予約語的なキーワードは気を付けた方が良いです（色々実験してみたところ、"LOW", "HIGH" なども同様に範囲指定の意味として解釈されてしまいました。たぶん他にもあると思います）

解決策

今回の例では "OTHER" を範囲指定の意味ではないと明示してあげる必要があります。

なので、以下のようにHLOという変数を作ってNullに設定しておけば良いはず。

（HLO？という方は「PROC FORMAT入門9 : CNTLIN=オプション」を参照）

data FMT; input FMTNAME:$10. START:$10. LABEL:$20. HLO:$1.; cards; $TEST_ OTHER その他 . ; proc format cntlin=FMT; run;

上でフォーマットを作り直したので、再度このフォーマットをあてたデータセットを見てみましょう。

はい！さっきまで2行目の欠損値が「その他」になっちゃってましたが、今度は問題なさそうです。

最後に

私が把握しているのは "OTHER", "LOW", "HIGH" のみで、その他にFORMATプロシジャの予約語的なキーワードに引っかかるものがあるかもしれません。

今回の場合だと"OTHER"に対してHLOという変数を作ってNullに設定することで解決しましたが、その他の私が把握していないキーワードに対して、各自SASのリファレンスや「PROC FORMAT入門9 : CNTLIN=オプション」などを参考にして、他の必要な変数を設定してあげる必要があります。

SASで花火を描く

暑中お見舞い申し上げます。

SASで花火描いてみました。三角関数をつかってるだけの単純なプログラムです。

/* 花火描画用データ */ data fireworks;  rad = constant('pi')/180;  call streaminit(777);  do i = 1 to 5;   xpos = rand("uniform",-30,30);   ypos = rand("uniform",-30,30);   do rmax = 0, 36;   do r = 0 to rmax by 4;   do a = 0 to 360 by 15;     x=cos(a*rad)*r+xpos;     y=sin(a*rad)*r+ypos;     output;   end;   end;   end;  end; run; /* 花火描画 */ title; ods graphics / height=20cm width=20cm; options animation=start animduration=1.2 printerpath=gif; ods noresults; ods printer file="出力するGIFファイルのフルパスをここに指定。ファイル名は「fireworks.gif」とか適当に" ; %macro m_fireworks(i=,color=);  proc sgplot data=fireworks aspect=1 noautolegend;   styleattrs wallcolor=midnightblue;   scatter x=x y=y / markerattrs=(symbol=circlefilled color=&color size=0.1cm );   xaxis min=-50 max=50 display=none;   yaxis min=-50 max=50 display=none;   inset "Summer Greetings to All."/ position=bottom textattrs=(color=lightblue size=25cm family="Arial");   by rmax;   where i=&i;  run; %mend; %m_fireworks(i=1,color=yellow); %m_fireworks(i=2,color=skyblue); %m_fireworks(i=3,color=lightgreen); %m_fireworks(i=4,color=yellow); %m_fireworks(i=5,color=lightpink); ods printer close; ods results; options animation=stop; ods graphics / reset=all;

「PUT "WARNING-"」みたいな書き方。

今回は、以下で紹介した記事のおまけ的内容。

• ログに独自のエラーメッセージを表示する。
• 「PUT」や「%PUT」でログに独自のエラーメッセージを表示するときの定石

最初に

上の記事のおさらい。

以下サンプルのように「NOTE」「WARNING」「ERROR」のいずれかのキーワードとコロン「:」をつけてPUTすると、ログ上に色つきの独自メッセージを表示できる。

data _null_;   put "NOTE:    独自メッセージ1";   put "WARNING: 独自メッセージ2";   put "ERROR:   独自メッセージ3"; run; ログ

本題 

上の例でコロン（:）ではなくハイフン（-）に変えると、、

data _null_;   put "NOTE-    独自メッセージ1";   put "WARNING- 独自メッセージ2";   put "ERROR-   独自メッセージ3"; run; ログ

「NOTE」「WARNING」「ERROR」のキーワード無しで色つきのメッセージを出力できます。

活用例

これは、以下のような使い方が出来ます。

data _null_;   put "WARNING: 独自メッセージ1";   put "WARNING- 独自メッセージ1の続き"; run; ログ

WARNINGメッセージを改行しているように見せることが出来ます。

ただし、1つだけ残念ポイント。

%PUTで同様のテクニックを使おうとすると、、

%put WARNING: 独自メッセージ1; %put WARNING- 独自メッセージ1の続き; ログ

ログを見てみると、間にプログラムの文が挟まっていて、WARNINGメッセージが分割されて見えてしまうところ。残念よ。

「PUT 変数名=」という書き方

PUTステートメントの基本的なテクニックの紹介。

基本的なテクニック

data test;   set sashelp.class;   if age=13 then put height=; run; ログ  height=56.5  height=65.3  height=62.5

「変数名=」という書き方で、例のように先頭に「変数名=」という書式のテキストを追加できます。

以下のように、ちょっと勿体ない書き方をしがち。

ちょっと勿体ない例

data test;   set sashelp.class;   if age=13 then put "height=" height; run; ログ  height=56.5  height=65.3  height=62.5

勿体ない！