非商用目的で 無料で使えるSAS OnDemand for Academicsについて，始め方がわからないという方は，あるいはインストール版のSAS9.4について，インストールの仕方にとまどうって方が多いということで，利用手順をスライドにしてアップしてくれた方がいるので紹介させてくださいませ
https://www.docswell.com/s/NorihiroSuzuki/KNLD4E-SAS_ODA_9.4_installguide#p1
SASユーザー総会 論文集
1982年～2024年までのSASユーザー総会論文集について，公開場所が分散していたりしていましたので世話人会のほうで一つ所にアーカイブしました
SASユーザー総会論文集アーカイブ
これは高橋 行雄 先生が2014年に当時の世話人会で行った，電子化・公開の取り組みを引き継いだものになります
目録のページには，全演題と，SAS社で公開されている発表スライドへのリンクをつけていってますが，確認しながら手作業でつけていっているので，まだ2000～2016年ぐらいまでです
また目録データはエクセルファイルでダウンロード可能です．
OCR由来なことと，名寄せ処理はおこなっておりませんが，簡単な分析はできると思います
エクセル取り込んでもらって
proc sgplot data=sugi;
where author_cat="筆頭筆者";
vbar year;
xaxis values=(1982 to 2024) valueattrs=(size=7);
run;
とかで 演題数の推移など
さて，2025年もユーザー総会について演題募集が始まっておりますのでこちらもぜひ
https://sas-user2025.ywstat.jp/
REGプロシジャでCLASSステートメントがなかったり，交互作用の事前処理が必要な点はHPREGで改善されているよという話
REGプロシジャにはCLASSステートメントがないため，カテゴリ変数をモデルに含める場合，事前に数値型でダミー変数を用意しないといけなかったり X*Yといった交互作用の記述をmodelステートメント内にかけないことから，例えば
data class1;
set sashelp.class;
sexn = choosen(whichc(sex,"男子","女子"),0,1);
age_height=age*height;;
run;
proc reg data=class1;
model weight = age height sexn age_height;
run;
quit;
などとして
こういう感じに解析したりしてるわけですが
HPREGだと
proc hpreg data=sashelp.class;
class sex;
model weight = age height sex age*height;
run;
quit;
基本的にはREGプロシジャで，できることを押さえた上で，機能拡張されているような印象ですが，REGにあってHPREGにない機能もあり，似たことのできるGLM,GLMSELECTとの機能比較などは公式リファレンスにあるのでそちらを参照ください．
「PROC HPREG Contrasted with Other SAS Procedures」
https://documentation.sas.com/doc/en/stathpug/15.2/stathpug_hpreg_overview02.htm
テキストをモールス信号に変換する関数
やっぱり モールス信号って業務で良く使いますよね．
プログラマ派遣先にいったら，全ての通信機器を取り上げられ，鉄格子窓のある部屋に閉じ込められて，SASの入ったPCだけ与えられて労働させられたみたいなあるあるの状況の際に
電気の消灯点灯タイミングで外部に助けを求めることができるかどうかが大事です
いちいち調べるの面倒なので，変換関数作ってみました
proc fcmp outlib=work.functions.common;
function text2morse (text $) $ ;
length key1 $1 text _rvalue rvalue $2000 ;
_text=upcase(compress(text));
declare dictionary d1 ;
d1["A"] = "・－";
d1["B"] = "－・・・";
d1["C"] = "－・－・";
d1["D"] = "－・・";
d1["E"] = "・";
d1["F"] = "・・－・";
d1["G"] = "－－・";
d1["H"] = "・・・・";
d1["I"] = "・・";
d1["J"] = "・－－－";
d1["K"] = "－・－";
d1["L"] = "・－・・";
d1["M"] = "－－";
d1["N"] = "－・";
d1["O"] = "－－－";
d1["P"] = "・－－・";
d1["Q"] = "－－・－";
d1["R"] = "・－・";
d1["S"] = "・・・";
d1["T"] = "－";
d1["U"] = "・・－";
d1["V"] = "・・・－";
d1["W"] = "・－－";
d1["X"] = "－・・－";
d1["Y"] = "－・－－";
d1["Z"] = "－－・・";
do i = 1 to length(_text);
_key1=char(_text,i);
_rvalue= d1[_key1];
rvalue=catx(" ",rvalue,_rvalue);
end;
put _text= rvalue=;
return(rvalue);
endsub;
run;
で，使ってみます
options cmplib=work.functions ;
data test;
length x $200.;
x = "SOS";output;
x = "SAS";output;
x = "VVV";output;
run;
data a;
set test;
y = text2morse(x);
run;
proc print data=a noobs;
run;
こっちのけんと さんのおかげで誰でもわかるようになったSOSとかちゃんと変換できてますね
VVVはムスカ大佐が通信してたやつで，本日は晴天なり みたいな試験信号です
数字や記号は今回無視しています
こっちのけんと さんのおかげで誰でもわかるようになったSOSとかちゃんと変換できてますね
VVVはムスカ大佐が通信してたやつで，本日は晴天なり みたいな試験信号です
数字や記号は今回無視しています
Excelのガントチャートでイベントの持続期間を可視化
👆
こういうのを簡単に作る方法．DDEとかは使わない．SASのRWI(Report Writing Interface)とODS EXCELでやります
テストデータとして，
data wk1;
SUBJID="AAA";AETERM="事象1"; AESTDY=1;AEENDY=5;output;
SUBJID="AAA";AETERM="事象2"; AESTDY=3;AEENDY=8;output;
SUBJID="BBB";AETERM="事象1"; AESTDY=2;AEENDY=30;output;
SUBJID="BBB";AETERM="事象2"; AESTDY=-2;AEENDY=3;output;
SUBJID="BBB";AETERM="事象3"; AESTDY=4;AEENDY=.;output;
SUBJID="CCC";AETERM="事象1"; AESTDY=15;AEENDY=20;output;
run;
要は，開始と終了の相対日付（DY変数）をもっていればなんでもOK
終了日がnullの場合は，継続として最大までバーをひっぱる仕様を想定
以下 コードですが，面白いのがexcelに出力するSAS RWIの部分での指定が，信じられないぐらいシンプルなのと，セル幅が自動でlengthに連動してくれるので細かいメッシュにたいして何の範囲指定や命令を書かなくても全部えー感じにしてくれるところ．
この技術はわりと使える
proc sql noprint;
select min(min(AESTDY,1)) into:min trimmed from wk1;
select max(max(AESTDY,AEENDY)) into:max trimmed from wk1;
quit;
%put &=min;
%put &=max;
%let range =%eval( %sysfunc(range(&min,&max))+ (0<&min));
%put &=range;
data wk2;
set wk1;
array ar{&range.} $1.;
if &min < 1 then do;
st = AESTDY - &min + (AESTDY<1);
if ^missing(AEENDY) then en = AEENDY - &min + (AEENDY<1);
end;
else do;
st=AESTDY;
en=AEENDY;
end;
do i = 1 to ⦥
if st <=i <=en then ar{i}="Y";
else if st=i and missing(en) then ar{i}="Y";
else if st<=i and missing(en) then ar{i}="O";
end;
run;
%macro hed();
%do i = &min %to &max;
if &i ne 0 then ob.format_cell(data: "&i", style_attr: "background=blue color=white");
%end;
%mend;
%macro row();
%do i = 1 %to %eval(&range);
if ar&i="Y" then ob.format_cell(data: ar&i, style_attr: "background=red color=red");
else if ar&i="O" then ob.format_cell(data: ar&i, style_attr: "background=pink color=pink");
else ob.format_cell(data: ar&i, style_attr: "background=white color=white");
%end;
%mend;
ods excel options( sheet_name= "GANT" );
data _NULL_;
set wk2 end=EOF;
if _N_=1 then do;
dcl odsout ob();
ob.table_start();
*** header ;
ob.head_start();
ob.row_start();
ob.format_cell(data: "ID", style_attr: "background=darkblue color=white");
ob.format_cell(data: "Event", style_attr: "background=darkblue color=white");
ob.format_cell(data: "STDY", style_attr: "background=darkblue color=white");
ob.format_cell(data: "ENDY", style_attr: "background=darkblue color=white");
%hed();
ob.row_end();
ob.head_end();
end;
*** report ;
ob.row_start();
ob.format_cell(data: SUBJID);
ob.format_cell(data: AETERM);
ob.format_cell(data: AESTDY);
ob.format_cell(data: AEENDY);
%row();
ob.row_end();
if EOF then ob.table_end();
run;
ods excel close;
Proc ODSTABLEでRTFを作る際に任意の位置で改ページをする方法
Proc OdstableはProc Reportでできる機能をほぼすべて実装可能なうえに，構文がよりシンプルであり,個人的には完全な上位互換だと思っています
ただ，RTFを作成する際の大きな弱点として，Proc ReportのbreakがOdstableにはなく任意の位置での改ページができないことがありました．
2018年にSASユーザー総会で「PROC ODSTABLEを用いた帳票作成」という論文が発表されていますが，そちらでもその点書かれていました．
8年間色々，チャレンジしてきましたが，最近また向き合う機会があり，一応実装できたので共有．
data output;
set sashelp.class;
if mod(_N_,3)=0 then break ="Y";
keep Name Age break;
run;
↑
break="Y"の直前のレコードまで出力して，Yのレコードから次ページにうつるようにしてみます
ods rtf ;
proc odstable data=output ;
column Name Age break ;
cellstyle break="Y" as data{PRETEXT="(*ESC*)R'\pagebb' " };
define Name; end;
define Age; end;
define break; print=no; end;
run;
ods rtf close;
なんとたったのこれだけで
みたいに，改ページできていることがわかります
結局種がわかってしまえば簡単で，インラインフォーマット関数に改ページがあることと
print=noでも，RTFに埋め込みさえすれば作動するのでcellstylesのPRETEXTで入れてしまえばいいということでした
開始終了日の情報から持続期間をまとめる話
まあ，医薬なら有害事象の持続とか，薬剤の継続使用を特定のカテゴリでまとめて
切れ間なく続いていた期間をまとめたレコードを合成するってことはあるのですが
これが結構SASだとやりにくい．SQLでもやれるんだけど，レコード一時的にかさむから
環境のパワーないとちょっと辛い時がある
なので結局，ハッシュオブジェクトを使ってしまう
👆適当なサンプルですが，症例ごとにカテゴリで，切れ間なく持続してる期間でまとめて
レコードをおこす
data wk1;
length ID cat term $30.;
format st en yymmdd10.;
ID="001";cat="痛み";term="頭痛";st="01JAN2025"d;en="05JAN2025"d;output;
ID="001";cat="痛み";term="筋肉痛";st="06JAN2025"d;en="10JAN2025"d;output;
ID="001";cat="痛み";term="関節痛";st="12JAN2025"d;en="13JAN2025"d;output;
ID="001";cat="痛み";term="神経痛";st="14JAN2025"d;en="19JAN2025"d;output;
ID="001";cat="胃腸障害";term="胃痛";st="01JAN2025"d;en="08JAN2025"d;output;
ID="001";cat="胃腸障害";term="腸炎";st="06JAN2025"d;en="12JAN2025"d;output;
ID="002";cat="感染症";term="インフルエンザ";st="02JAN2025"d;en="08JAN2025"d;output;
ID="002";cat="感染症";term="感冒";st="01JAN2025"d;en="10JAN2025"d;output;
run;
ここで，setしつつ，setしたもの自体をハッシュオブジェクトにマルチキー許容でいれて
自症例の全検索をかけちゃいます．ハッシュオブジェクトはメモリ内走査なので，特にアルゴリズム工夫なく力技で十分いけます
proc sort data=wk1;
by ID cat st en;
run;
data wk2;
length ID cat $30. cat_st cat_en 8.;
format cat_st cat_en yymmdd10.;
set wk1;
if 0 then set wk1(rename=(st=_st en=_en));
if _N_ = 1 then do;
declare hash h1(dataset:"wk1(rename=(st=_st en=_en))",multidata:"Y");
h1.definekey("ID","cat");
h1.definedata("_st","_en");
h1.definedone();
end;
cat_st = st;
cat_en = en;
do while(h1.do_over() = 0);
if n(_st,_en) = 2 then do;
if _st-1<=cat_en<=_en and cat_en <= _en then do;
cat_en = _en;
if _st < cat_st then cat_st = _st;
end;
end;
end;
keep ID cat cat_st cat_en;
run;
do_overメソッドはキーがヒットする限り0を返し続けハッシュオブジェクトの中を探しつづけます．なのでdo whileと相性がいいですね．
そのままだと，元レコード数と変わらず集約されてないので，一番包括的なものだけを残します
proc sort data=wk2 ;
by ID cat cat_en cat_st;
run;
data wk3;
set wk2;
by ID cat cat_en cat_st;
if first.cat_en;
run;
どうでしょう？そんな 難しくないと感じるのではないでしょうか
どうでしょう？そんな 難しくないと感じるのではないでしょうか
登録:
投稿 (Atom)