Proc Luaの世界①-こんにちはからif文あたりまで

Proc Luaは間違いなく、相当に便利な機能です。今までSASマクロが担ってきた機能の多くはProc Luaでより簡易に記述することができます。
前の記事を読んでいない方は先にご一読ください

「Proc Luaの衝撃」
http://sas-tumesas.blogspot.jp/2016/08/proc-lua.html


今後は、具体的かつ実践的な例をあげる記事と、基礎から一緒に勉強する記事を織り交ぜていきたいと思います。

僕も超初心者なので、多分間違ったことも書いてしまうと思いますが、一緒に切磋琢磨していきたいので適宜指摘や情報交換いただけますと幸いです。
Luaそのものについては、既に良質な情報が公式、非公式を問わずに出揃っているのですが、Proc LuaはあくまでSASの中で動くものなので、ちょっとSASプログラマーにしかわかりにくところがあります。
なので僕が一から解説することも無意味ではないはず、、。ということでいきましょう！

まず最初は以下のコードです

proc lua;
 submit;
print('こんにちは世界')
 endsubmit;
quit;

実行するとログにこのようにでます。

まず、Luaプロシジャはproc luaで始まりquitで終わります。
そしてsubmit;からendsubmit;の間がLuaの世界で、ここにはLuaの文法で記述されたコードが
実行されます。

print関数はログ画面に引数を出力する関数になります。
Luaでは文字列はコーテーションで括りますが、シングルでもダブルでもOKです。
そして、SASと違ってセミコロン「;」は要りません。まあ、つけてもLuaでは空文の意味なので
影響ないですが。

つづいてコメントのやり方について、以下をみてください
proc lua;
/*submitまでのこのスペースはまだSASの世界なのでSASのコメントができる*/
submit;
  --ここからはLuaのやり方でコメントする
  --1行コメントはこのようにバー2連続で書く
print('こんにちは世界')
  --[[
    複数行コメントは
このように書きます
  --]]
 endsubmit;
quit;

とまあ、こんな感じですが、SASの拡張エディタだと色変わらないんでちょっと嫌ですね。

次は変数についてです。

proc lua;
 submit;

  local a =1
  local b='あああ'

print("a=" ..a)
print("b=" ..b)

 endsubmit;
quit;

結果は

local a=1でローカル変数aを定義しそのまま値を代入しています。
local aとして値を入れない宣言もありです。
後で説明しますが、基本localでの定義を推奨します。
これも後で説明しますが、Luaには変数の型というものはなく、
入っている値によって型が動的に決まります。
入っている値が数値なら数値型、文字なら文字型ということです。
print("a=" ..a)の..ですが、これは文字列の連結を表します。SASでいうところの||と同じですね


変数についてもう少し見てみます

proc lua;
 submit;

    --こうも書ける
  local c ,d =1 , 2

print("c=" ..c)
print("d=" ..d)

--こうやって中身を入れ替えることもできる
local c ,d = d,c

print("c=" ..c)
print("d=" ..d)

--型はtype関数で調べれる
print("cの型は"..type(c).."です")

--Luaは動的型つき言語→つまり変数は型を持たず、値のみが型を持つ。
c="A"
print("c=" ..c)
print("cの型は"..type(c).."です")

--未初期化の変数には「nil」が入っている
print(zzz)

 endsubmit;
quit;

結果は

です。
 local 変数,変数,変数... =値,値,値...といった書き方がかけるということ
local c ,d = d,cで値が入れ替えられるということ。
またtype関数で型が調べられること。そして型が値によって変わるということを押さえてください。

またnilという概念があって、要するにnullみたいなもんなんですが、未初期化の変数にはnilが入ってます。
いずれまた詳しく掘り下げますが、Luaでは、論理演算において「false」と「nil」だけが偽それ以外は全て真というルールで動きます。つまり0だろうが、長さ0の""のような文字列だろうがそれらを真偽値では真になるので注意です。また関数などでも、指定しなかった引数にはnilがはいるので注意です。

さて、localでローカル変数が定義できるということはグローバルも当然あります。これは何もつけずにそのまま値を代入するとグローバル変数になります。グローバル変数になるとproc lua; quit;が終了してもその変数が残ります。つまり

proc lua;
 submit;
  --local宣言せずに変数をつくるとGlobalになっちゃう。多分それは避けた方がいい
  e = 100
 endsubmit;
quit;

一旦プロシジャを閉じても

proc lua;
 submit;
  --global変数はずっと生きちゃうから
  print ("e="..e)
 endsubmit;
quit;

と変数が生きちゃってます。nilを入れない限り、ずっと参照できます。
それをうまく使うことも当然できますが、しばらくはグローバルで定義せずにひとつのproc luaからquitの間で完結するコードを書いていこうと思います。

つづいて、算術演算子に関しては、まあ普通な感じなので、さっと見て覚えてください

/*算術演算子*/
proc lua;
submit;
 local a ,b, c, d;
 print(3+2)
 print(3-2)
 print(3/2)
 print(3%2) --剰余
 print(3^2)
endsubmit;
quit;

関係演算子については、ちょっと間違いやすいです

/*関係演算子*/
proc lua;
submit;
 local a =1
 local b =2
 local c ='1'
 print(a==1) --★SASと違うから要注意！=1だとエラーになる。
 print(a>b)
 print(a>=b)
 print(a<b)
 print(a<=b)
 print(a~=b) --★^=じゃなくて~=なので注意
 print(a==c) --★型が違うと不等とみなされるから注意！
 endsubmit;
quit;

SASは割り当ての際も=をつなぐし、等しいことも=で表現しますがLuaの等しいは==で表現することに要注意です。僕は既にこれで10回はエラーだしました


最後にIF文について紹介します。
基本的に
IF 条件式 THEN 処理　END　と書きます。
SASだとendがつくのは処理が複数で then doの場合のみですがLuaは必ずENDがいるので注意です。
例をみます。

/*if文*/
proc lua;
submit;
if sas.weekday(sas.today()) == 2 then
print('今日は月曜だ')
end  --★if文の終わりに絶対endが必要なことに要注意
endsubmit;
quit;

結果は日によって変わるので省略。

また後日説明しますが、sas.SAS関数名でSASの関数を使うことができます。

proc lua;
submit;
if sas.weekday(sas.today()) == 2 then
print('今日は月曜だ')
print('もう一度言おう今日は月曜だ')
end
endsubmit;
quit;

then からendの間にはいくらでもかけます。

proc lua;
submit;
if sas.weekday(sas.today()) == 1 then
 print('今日は日曜だ')
elseif sas.weekday(sas.today()) == 7 then --★else ifじゃなくてelseifなことに注意
          print('今日は土曜だ')
else
 print('今日は平日だ')
end
endsubmit;
quit;

elseifで条件分岐をいくらでも深くできます。elseifとスペースをいれずに1キーワードで記述することに注意です。

今回はここまで。

Proc Luaの衝撃

先日、下記のようなコードを見かけて、久々に震えるほどの感激を受けました。

proc Lua;
 submit;
  if sas.exists("WORK.CLASS") then
   print("NOTE:データセットCLASSは既にあるので何もしません")
  else
   print("NOTE:データセットCLASSがなかったから作成しました")
   sas.submit([[data CLASS; set sashelp.CLASS ; run;]])
  end
 endsubmit;
run;

実行されるコードの分岐がマクロを使わずに記述されてる......。

今までなら例えば

%macro test;

 %if %sysfunc(exist(WORK.CLASS)) eq 1 %then %do;

  %put NOTE:データセットCLASSは既にあるので何もしません;

 %end;

 %else %do;

  %put NOTE:データセットCLASSがなかったから作成しました;

   data CLASS; set sashelp.CLASS; run;

  %end;

%mend test;

%test

って書いてた処理ですね。

どうでしょう？SAS初心者、或いは全く知らない人の気持ちになってコードを見た時に
どっちが見やすいですか？わかりやすいと感じますか？

僕は前者だと思います。

proc luaというものについて何も知らない状態で最初のコードを見ましたが
何の処理をしているのか、実行した時にどうなるかが一瞬で理解できました。

で、肝心のProc Luaって何ですか？って話ですが、Luaっていうプログラミング言語があるんですね。
スクリプトタイプの言語で簡易な記述と、高速な処理、他の言語への組み込みのしやすさが特徴みたいです。

これは、正直惚れてしまいました。ちょっと真剣にLuaをゼロから勉強してみようと思います。
SASプログラミングは、なまじSASマクロがなんでもできるだけに、SASマクロに頼り過ぎていて、簡易な記述、合理的な構造から遠ざかっている側面があるのではないか？
それがSASを苦手な人、SASアレルギーの一因になっているのでは？って思ったりもするので
proc Luaを知ったのは僥倖でした。

まだ勉強し始めですが、やはり結構簡単な上に、処理をわかりやすくできるので
実際の業務にも取り入れやすいんじゃないかと思います。マクロでどうしょうもなくなったプロジェクトの銀の弾丸になりえるんじゃなかなぁって淡い希望を抱いてます。

ただ惜しむらくは、9.4のメンテナンスリリース3からの正式導入なことですね（とりあえず9.4なら試用版で使えますが）。

でも、コレはこれからかなりくる技術だと思います。だって普通に便利だから。いや、正直DS2より効果がわかりやすい気が…（オイッ）

というわけで暫くproc Luaの話が増えると思います。9.4の最新なんか入ってね～しって人が多いとは思いますが、未来のため、後ろに続くこれからのSASプログラマー達が仕事しやすくなるようにみんなで勉強してみませんか？

RTF出力の新視点－STREAMプロシジャの応用力は侮れないという話

RTF出力、実は未だあんまり経験ないんですが、色々大変なのはわかります。

多分おおまかに2つ方法があって、一つはODS RTFを利用する方法、もう一つはputでタグ情報と共にデータを出力して作ってしまう方法だと思います。
どちらも過去のユーザー総会等で素晴らしい発表が多くでているので、そちらを参照してください。

まあ、それで終わりの話なんですが、僕はひねくれているので、何か他に面白い方法はないかなと考えてました。

上記の2つの方法もファイルの新規作成なので、書式等の調整も全てコードでやる必要があります。
要するに、あらかじめ用意したRTFファイルに値を差し込む方法ではないということです。
（一応DDEを使った既存WORDファイルへの値入れする方法もあるみたいですが、少なくとも日本ではあんまり主流ではない気がします）

そこで、書式設定等を済ませたテンプレートRTFにSASから出力する新しい方法について紹介します。

ズバリproc streamを使います! 9.4からの新プロシジャです。一応9.3でも正式版ではないですが実行できます。

proc stream？って方は過去記事
「proc stremの世界」
http://sas-tumesas.blogspot.jp/2015/07/proc-strem.html

を見てください。

要はマクロ変数を展開したテキスト生成ができるプロシジャで、普通の平テキストはもとより、HTMLやXML、RTFなどタグ言語の生成に強みがあるんですね。

さあ、前置きが長くなりましたが、やっと本題に行きます。


まず、ノートパッドでもMS WORDでもいいので、以下の様な出力テンプレートを拡張子RTFで作成、保存してください。(例は「テンプレート.rtf」という名前で保存)

ちなみに値を差し込みたい部分をマクロ変数で記述します。

で、例えば以下の様なデータセットがあった場合

data _NULL_;
set SALES end=eof;
call symputx(ITEM,SALE);
TOTAL+SALE;
if eof then call symputx('TOTAL',TOTAL);
run;

とでも書いて、とりあえずテンプレートのマクロ変数名に合わせて、マクロ変数を作ってやります。

そうしたら、いよいよ次のように書きます

filename in "/home/sasyamasasyama0/テンプレート.rtf" lrecl = 32755;;
filename out "/home/sasyamasasyama0/アウトプット.rtf" lrecl = 32755;;
proc stream outfile=out resetdelim="rdlm" quoting=single;
begin
rdlm; %include in;
;;;;

これだけです。

実行すると、「アウトプット.rtf」ができて、それを開いてみると

と、値が展開されていることがわかります。

どうでしょう？
結構面白いですよね。

resetdlm=で、リセットデリミタを指定しています。リセットデリミタの後に%includeを記述しないと、ただのテキストとして「%include」ってでちゃいます。
%include　ファイル名で、指定したファイルをproc streamの中に展開します。そこでその中に含まれるマクロ変数が展開されるわけです。

こういうレイアウトが固定されたrtfに値を出すのには結構向いていると思います。
まあ、多分ご想像されているとおり、表がデータに応じて可変だったり、出力する値の数が多いものには正直オススメできないですね。
（できなくはないけど、それに応じて表を構成するためのタグを作って一緒に出す必要があるので、putで出す際のコードをマクロ化してstreamの中で展開するイメージになるので、お手軽さがなくなる）

定型のちょっとしたレポート向きだと思いますが、知っていて損はない使い方だと思います。

以上です。

monotonic関数の性質。_N_とはちょっと違うという話

monotonic関数は色々注意点もありますが、中々面白い関数です。

例えば以下のデータセットがあります

data Q1;
X=1;output;
X=1;output;
X=2;output;
X=3;output;
X=2;output;
X=2;output;
X=4;output;
X=3;output;
X=3;output;
run;

それで、偶数のオブザベーション番号だけproc printしろっていわれたら、

proc print data=Q1;
　where mod(_N_,2)=0;
run;

なんて書きたくなりますが、もちろん_N_はwhereでは使えません。

ところが

proc print data=Q1;
　where mod(monotonic(),2)=0;
run;

ならいけます

これは何回も紹介している話なんですが、上記のことから
monotonic関数はwhereでも使える_N_みたいなものと認識してしまいがちなのですが
それはちょっと違います

例えば

data A1;
set Q1;
Y=monotonic();
run;

は

となり

Y=_N_;とした場合と全く同じですが

以下のコードはどうなるでしょう？

data A2;
set Q1;
if x=1 then Y=monotonic();
if x=2 then Y=monotonic();
if x=3 then Y=monotonic();
if x=4 then Y=monotonic();
run;

実は

と、Xの値ごとに連番が振られるんですね。

つまりmonotonic関数は、行番号を返しているのではなく、
monotonic関数を通って処理される回数をカウントして返しているようなイメージということなんですね。（理解が間違っていたらご指摘ください）

これは非常に面白い性質ですね。

パッケージ内で別パッケージのメソッドを使用するための方法-コンポジション

今回の記事については説明や用語・概念の説明についてちょっと自信がないです。
SASメインで、他のオブジェクト指向の言語をあんまりやったことがないので、もしもっと上手く説明できる方がいらっしゃれば是非どこかにアップして欲しいです。

さて、苦手ながらも大事な部分なので頑張って書いてみます。

一般的なオブジェクト指向言語でのクラスという概念が、DS2ではパッケージになると思います。

DS2のパッケージには、いわゆる「インヘリタンス（継承）」という機能がありません。

継承は、既存クラスの全機能を引き継いで新たなクラスを定義することで、オブジェクト指向の3大要素とかって言われるぐらい重要なものです。親のメソッドを全部引き継いだ子を作るイメージです

それがないってどうなの？ってことですが、代わりにすべて、同じくオブジェクト指向の「コンポジション（集約）」という概念でまかないます。

インヘリタンスとコンポジションの違いは、ちょっと検索して貰えればたくさん説明がでてくるので、そちらを参考にしてください。

僕がわかりやすかったのは以下のとか？

http://promamo.com/?p=4291

http://www.kab-studio.biz/Programing/JavaA2Z/Word/00000070.html

要するに、単純に他のパッケージの必要なメソッドだけもってきて新規パッケージの中で使えるようにしちゃうってイメージで僕は理解しました。

例えば以下のような、数値に1を加えるメソッド(tasu1メソッド)と2倍するメソッド(kake2)をもったパッケージ(p1パッケージ)があるとします

proc ds2 libs=work;

package p1 /overwrite=yes;

declare double v;

method tasu1(double v) returns double;

return v+1;

end;
 method kake2(double v) returns double;
  return v*2;
 end;

endpackage;

run;

quit;

次に新しく作成するパッケージ(p2パッケージ)の中で、p1のtasu1メソッドの機能を使いたい場合は以下のように書けます

proc ds2 libs=work;

package p2 /overwrite=yes;

declare package p1 pp1();

declare double v;

method p2();

pp1=_new_ p1();

end;

method tasu1(double v) returns double;

return pp1.tasu1(v);

end;

method tasu2(double v) returns double;

return v+2;

end;

endpackage;

run;

quit;

最初のdeclareでパッケージp1をpp1という名前で呼び出している（パッケージp1からインスタンスpp1を生成）のは、普通にパッケージを利用する場合と同様なので自然にかけると思います。
ただ、それだけだと不十分で、

method p2();

pp1=_new_ p1();

end;

の部分、パッケージ名と同じなのでコンストラクタになるわけですが、
ここで謎の式を書いてますが、これも実はパッケージの呼び出し方（インスタンスの生成）の書き方なんですね。p1パッケージからpp1インスタンスを生成するという意味になります。

declareステートメントと同じことしてんじゃんって感じで、なんでこんなのが必要なのかなんですが
呼び出しのコードをみるとわかりやすいかもです

proc ds2 libs=work;

data A1(overwrite=yes);

declare double x y;

declare package p2 p2();

method run();

x=p2.tasu1(1);

y=p2.tasu2(1);

end;

enddata;

run;

quit;

ちなみに結果は

なんですが、ようするにp2パッケージしか呼び出してないのに、内部でp1パッケージの機能も呼び出されている、つまりp1のインスタンスも生成されているのがミソなんです。
パッケージが呼び出された際に、最初に実行されるコンストラクタの部分でp1のインスタンスを作らないと実際使えないよってことなんですね。パッケージを呼び出した際に内部のdecalreステートメントが実行されるわけではないので。


はぁ、自分の説明力のなさに嫌気がさします、、。
ほんとすみません、どなたか、もっとちゃんと説明しなおしてください

DS2でも外部ファイル出力はできます！！な話

ちょっとソース見つけられなかったんですが、だいぶ昔、海外のSASブログだったか何かのサイトで、「DS2にはfileステートメントがないみたいだぜ！」「マジかよ！ヤベェな！ハンドリング専門かよ！DS2はXXXX（なんか汚い言葉）だな！」みたいな外人同士のやりとりがあって、それをボーっと見てて、なんとなく頭に残ってました。

そうなんだぁ～、アウトプット部分は通常のSASコードで書かなきゃいけないんだ、そりゃXXXXだなぁって思ってたんですが、流石にそんなわけなかったですね。
普及し始めはガセじゃないですけど、みんな知らないんで結構間違ったコード上がるんで注意です。僕の記事に関しても間違いが含まれていることは多々あるはずなので鵜呑みは注意です。

話を戻しまして、実際問題ちょっと気づきにくいところなので説明します

じゃあまずはテストデータ

data Q1;
x=1;Y='A';output;
x=2;Y='B';output;
x=3;Y='C';output;
run;

そして、手始めに以下のコードを見てください

proc ds2 libs=work;
data _null_;
method run();
set Q1;
put X Y;
end;
enddata;
run;
quit;

こうすると、ログに

って出力されます。

うん、そりゃそうだろうね。データステップもそうじゃん。

データステップなら、putの前に、file print;ってつければアウトプットウインドウに
出力するし、先にfilenameで外部ファイルを定義してfileで指定すればそこに出力されます。

でも、DS2はfileステートメントないんだよね？じゃあやっぱ無理じゃん。

と思ったところで、次のコード見てください。

proc ds2 libs=work;
data ;
method run();
set Q1;
end;
enddata;
run;
quit;

？何が違うの？って思った方は、もう一度よく見てください。まずdata の後に_null_がなくなってます。旧データステップで、dataの後何もつけないと、何が起こるかは以下の記事を読んでください

「どうてもいい話　データセット名をかかずに実行」
http://sas-tumesas.blogspot.jp/2013/10/blog-post_1261.html

ようは勝手に連番でDATA１とか2とか作られていきます。
次にもう一点、putで何も指定してない！

これを実行すると何が起きるかというと
アウトプットウインドウに以下のような出力がでます

なるほど！そういうことか！SQLプロシジャでcreate tableつけずにseleｃｔ文流すのとおなじですね。抽出結果がそのままアウトプットになるわけですか。

もうここまできたら読めてきましたね。
つまり、外部ファイルへの出力はアウトプット　デリバリー　システム、そうODSでやれってことですよ。

例えばcsvで欲しければ

filename outf "aa.csv";
ods tagsets.csv file=outf;
proc ds2 libs=work;
data;
method run();
set Q1;
put X Y;
end;
enddata;
run;
quit;
ods tagsets.csv close;

とすればOKなわけです。

9.4からods excelがあるから、拡張子xlsxにしてtagsets.csvのところをexcelにすれば
みんな大好きEXCELファイルも作れます

以上です！

連続・重複した日付期間データをならして最小オブザベーション数にまとめる方法について

イベント開始日と終了日からなるデータがあって、発生ごとに縦積みされるとします。

ただし、諸々の事情で継続途中でオブザベーションが別れたり、重複して記録されてたりするとします。

たとえば

6/1開始　-　6/3終了

6/2開始　-　6/4終了

6/4開始　-　6/5終了

という3行のにわかれたデータを

6/1開始　-　6/5終了

の1行におこしなおすという処理です。

これは結構苦手な人が多いはず。

今回使うのは以下のデータ（ややこいので欠損は含まれないということにしましょう）

でIDとSTDTでソート済みとしましょう

data Q1;

ID=1;STDT='01JUN2016'd;ENDT='03JUN2016'd;output;

ID=1;STDT='02JUN2016'd;ENDT='04JUN2016'd;output;

ID=1;STDT='04JUN2016'd;ENDT='05JUN2016'd;output;

ID=1;STDT='06JUN2016'd;ENDT='09JUN2016'd;output;

ID=1;STDT='12JUN2016'd;ENDT='15JUN2016'd;output;

ID=1;STDT='13JUN2016'd;ENDT='14JUN2016'd;output;

ID=1;STDT='17JUN2016'd;ENDT='17JUN2016'd;output;

ID=2;STDT='01JUN2016'd;ENDT='05JUN2016'd;output;

ID=2;STDT='02JUN2016'd;ENDT='08JUN2016'd;output;

ID=2;STDT='03JUN2016'd;ENDT='04JUN2016'd;output;

ID=2;STDT='10JUN2016'd;ENDT='11JUN2016'd;output;

format STDT ENDT yymmdds10.;

run;

これを

上記のような形にします。

この手の処理を書く場合は、先にガントチャートみたいな図を書いてから考えると楽ですね。

頭の中だけで組むと、結構泥沼になったりします

ちなみにSAS on demandでは、プロダクト「SAS　OR」も使えて、そこに

proc ganttっていうガントチャートかけるプロシジャがあったので、使い方よくわからないけど

適当に指定して流してみました（頑張れば凄い綺麗なのが描けるみたいですけど）。

これのコードは本題じゃないので最後に書きます。

そこで、どうやるかですが、多分以下の感じで書くのが一般的ではないでしょうか？

data _Q1;

format r_STDT r_ENDT yymmdds10.;

set Q1;

by ID;

retain r_STDT r_ENDT ;

if first.ID then do;

r_STDT=STDT;

r_ENDT=ENDT;

end;

if STDT > r_ENDT+1 then do;

r_STDT=STDT;

r_ENDT=ENDT;

end;

if ENDT>r_ENDT then r_ENDT=ENDT;

run;

proc sql noprint;

create table A1 as

select ID,r_STDT as STDT,max(r_ENDT) as ENDT format=yymmdds10.

from _Q1

group by ID ,r_STDT;

quit;

1obs読み込みながらリテインしている開始日と終了日を条件に従って更新していくイメージです。

終了日はリテインしている終了日より後なら置き換える。開始日は、リテインしている終了日+1より

大きければそこで置き換える。

そしてできたデータセットについて、リテインしていた開始日でグループ化して、最も最大のリテインしていた終了日のみ残せばOKです。

sqlの部分に書き方はなんでもよくて、IDとr_STDT　r_ENDTでソートしてlast.r_STDTでしぼってもOKなわけです。

上記の処理は2ステップでやってます。SASは基本的に行の先読みがやりにくい言語なので

上記の考え方を1ステップで表現するのは、かなり難しいです。どこで切ってアウトプットするかが、先まで読まないと確定できないですからね。

ところが、ハッシュオブジェクトなら、簡単にできちゃうんですよね。

前の記事で紹介したmultidataとfind_nextを使えば、先読みみたいなことが。

data A2;

length _STDT _ENDT 8.;

retain r_ENDT;

set Q1;

if _N_=1 then do;

call missing(_STDT,_ENDT);

declare hash h1(dataset:'Q1(rename=(STDT=_STDT

                                                                      ENDT=_ENDT))',multidata:'Y');

h1.definekey('ID');

h1.definedata('_STDT','_ENDT');

h1.definedone();

end;

 l_ID=lag(ID);

 if STDT<r_ENDT and ID=l_ID then delete;

 else do rc = h1.find() by 0 while (rc = 0) ;

  if _STDT<=ENDT+1 and _ENDT>ENDT then do;

  ENDT=_ENDT;

  end;

  rc= h1.find_next() ;

 end ;

 r_ENDT=ENDT;

 keep ID STDT ENDT;

run;

実は僕は、この処理を考えるとき、ハッシュでのやり方は一瞬で考えてすぐに書けたんですが

最初の書き方が中々でてこなくて手こずりました。

慣れたらこっちの方がわかりやすいと思うんですが、どうなんでしょ。

ちなみに、最初から最後までSQL一本でも表現できますが、かんなり難しいコードになります。

開始日（条件付き）と終了日の全組み合わせを作ってから、中に含まれる重複期間を除外するイメージですが、興味のある方は「SQLパズル: プログラミングが変わる書き方/考え方　第2版」に全コード載っているので参考にしてください。

最後にproc ganttのコードですが、僕もよくわかっとりません。as=でスタート、af=でエンド時点を指定します。

チャートに重ねる文字とかの制御をラベルデータセットで制御するみたいなんですが、まあ、あんまりSAS OR使う機会もないし、調べるの面倒だったので適当です。

data label;

 format _LVAR $10.; 

 _Y = -1;

 _LVAR = "STDT";

 _XVAR = "STDT";

 _CLABEL = "BLUE";

 _YOFFSET = .8; OUTPUT;

 _LVAR = "ENDT";

 _XVAR = "ENDT";

 _YOFFSET = .8;

 _XOFFSET = -0.3;OUTPUT;

RUN;

data GQ1;

set Q1;

format STDT ENDT day.;

run;

PROC GANTT DATA = GQ1 LABDATA=LABEL;

 CHART / AS = STDT AF = ENDT NOLEGEND NOJOBNUM NOTNLABEL 

         TIMEAXISFORMAT= day.;

 ID ID STDT ENDT;

 ;

RUN;