SAS忘備録

2023年1月1日日曜日

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;

2022年11月29日火曜日

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のリファレンスに書いてない気がする。。

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

事前にフォーマットを除いておくか、もしくは以下のように「FORMAT=」で別の適切なフォーマットを便宜的にあてるなどの対応が必要になります。

proc sql noprint;

select X format=best16. into :X trimmed

from DT1;

quit;

%put &X;

ログ

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

data DT2;

X=123.456789;

run;

proc sql noprint;

select X into :X trimmed

from DT2;

quit;

%put &X;

ログ

123.4568

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

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

2022年7月9日土曜日

【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入門2 : 範囲の指定」を参照）

要は、以下のように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=オプション」などを参考にして、他の必要な変数を設定してあげる必要があります。

2022年6月29日水曜日

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;

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;

2022年3月31日木曜日

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

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

最初に

上の記事のおさらい。

以下サンプルのように「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メッセージが分割されて見えてしまうところ。残念よ。

2022年2月11日金曜日

「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

勿体ない！

2021年12月5日日曜日

「ODS OUTPUT」で変数属性が自動調整される件

「PROC SQL」と「ODS OUTPUT の PERSIST=RUN」のコンボ技の記事内で

ODS OUTPUTで出力すると変数属性（変数名、length、formatなど）が自動調整される場合がある

と言ったんですが、じゃあ実際どう自動調整されるのか。

結論からいって、

変数属性は自動調整されるけど、データの中身は無事のよう（リファレンスに詳しい挙動が書いてないので、環境やバージョンによって挙動が同じかは不明。。各自で挙動確認してみてください）
私はこのコンボ技を単純な頻度集計で使う程度で支障ないですが、変数属性が自動調整されてしまうのがまずい状況下では、ご注意下さい。

「ODS OUTPUT 」で変数属性がどう自動調整されるのか見てみる。

（実行日: 2021/12/03, 実行環境: SAS Ondemand for Academics [SAS 9,4M6]）

テストデータ

/* d1 */

data d1;

length a b $10. c e 8.;

label c = "test1";

format c yymmdd10.;

a="aaa"; b="bbb"; c=1; e=1;

run;

/* d2 */

data d2;

length a $20. b c d e 8.;

label c = "test2" d="test3";

format c d time5.;

a="aaa"; b=1; c=1; d=1; e=1;

run;

「PROC SQL の SELECTステートメント」の結果を「ODS OUTPUT」を使ってデータセットに出力してみます。

例①「データセットd1」を使ってODS OUTPUTで出力

ods select none;

ods output SQL_Results=out1;

proc sql;

select * from d1;

quit;

ods output close;

ods select all;

各データセットの変数属性

データセットd1

データセットout1 (ODS OUTPUTで出力したやつ)

例②「データセットd2」を使ってODS OUTPUTで出力

ods select none;

ods output SQL_Results=out2;

proc sql;

select * from d2;

quit;

ods output close;

ods select all;

各データセットの変数属性

データセットd2

データセットout2 (ODS OUTPUTで出力したやつ)

例③「データセットd1とd2」を使って「ODS OUTPUT の PERSIST=」でプロシジャ内の結果を結合して出力

ods select none;

ods output SQL_Results(persist=run)=out3;

proc sql;

select * from d1;

select * from d2;

quit;

ods output close;

ods select all;

各データセットの変数属性

データセットout1 (例①で、d1を使ってODS OUTPUTで出力したやつ)

データセットout2 (例②で、d2を使ってODS OUTPUTで出力したやつ)

データセットout3 (例③で、「ODS OUTPUT の PERSIST=」で出力したやつ)

結果の考察

formatを割り当てていない数値変数はなぜかbestフォーマットが割り当てられてた（これはどちらかというとSQLプロシジャ側のODS出力時の設定によるものと思われる。このようにODS出力時に勝手にformatが割り当てられる場合がある。割り当てられるformatはプロシジャによって様々）

例①の変数Eなどを参照

「ODS OUTPUT の PERSIST=」で出力結果を結合すると、

同じ変数名（かつ文字変数）でlengthが違う場合 = 大きい方に合わせられた

例③の変数Aなどを参照

同じ変数名で型が違う場合 = 片方の変数がrenameされた。

例③の変数Bなどを参照（結合後、変数Bと変数B2に分かれてる）

同じ変数名でラベルが違う場合 = 最初のデータセットのラベルになった

例③の変数Cなどを参照

同じ変数名でフォーマットが違う場合 = formatが削除された

例③の変数Cなどを参照

片方のデータセットにしか変数がない場合 = 変数がkeepされた。あとformatを割り当てた変数はちゃんとformatが残った。

例③の変数Dなどを参照

他にも私が把握しきれていないだけで「変数属性が自動調整」されるパターンはあると思います。

同じ変数名（かつ文字変数）のlengthが大きい方に合わせられたのと、同じ変数名で型が違う場合に勝手にrenameして両方保持してくれるのは、便利っていえば便利ですね。

登録: 投稿 (Atom)