フォレストプロットの作成について

[Kazuhiro]

初投稿です。よろしくお願いします。

以下のコードを用いてフォレストプロットを作成しようと考えているのですが

dattrmapを用いる際に変数が見つからないエラーが発生してしまいます。

データセットとして「forestdata」と「attrmap」は読み込ませています。

何か原因があればご教授いただけると幸いです。

proc sgplot data=forestdata noautolegend dattrmap=attrmap nowall noborder nocycleattrs;
styleattrs axisextent=data;
refline ref / axis=y lineattrs=(thickness=13 color=cxf0f0f0);
scatter y=forestval x=hazardratio / errorbarattrs=(thickness=1 color=blue) xerrorupper=HRUpperCL xerrorlower=HRLowerCL markerattrs=(size=0);

/*highlow y=forestval low=HRLowerCL high=HRUpeprCL / lineattrs=(color=blue);*/

bubble y=forestval x=HazardRatio size=N / bradiusmax=4.2 bradiusmin=0.2 fillattrs=(color=blue) lineattrs=(color=blue);
scatter y=forestval x=HazardRatio / markerattrs=(size=0) x2axis;
refline 1 / axis=x lineattrs=(pattern=shortdash) transparency=0.5;
text x=xl y=forestval text=text / position=center contributeoffsets=none;
run;

yaxistable subgroup / location=inside position =left textgroup=valtype labelattrs=(size=8) labeljustify=left titlejustify=left textgroupid=text indentweight=indent;
yaxistable grp1 grp2 / location=inside position=left labelattrs=(size=8) valueattrs=(size=7) valuehalign=center valuejustify=right;
yaxistable HR intpval / location=inside position=right labelattrs=(size=8) valueattrs=(size=7) valuehalign=center;

yaxis reverse display=none offsetmin=0.1 colorbands=odd colorbandsattrs=(transparency);
xaxis type=log value=(0.20.525) minor display=(nolabel) valueattrs=(size=7);
x2axis label='Hazard Ratio (95% IC)' display=(noline noticks novalues) lavelattrs=(size=8);
run;

【参考】

魚住先生の報告

https://www.sas.com/content/dam/SAS/ja_jp/doc/event/sas-user-groups/usergroups2016-d-10.pdf#search='...

[yu_sas]

こんにちは！

変数が見つからないというエラーですが、変数xlやtextが見つからないということであれば、付録Cと付録Dで示されているデータセットを結合していないように思います。付録Bに書かれているように、これらは事前に結合されている必要があると思います。

※付録DでATTRMAPと書かれていますが、内容からこれはxaxisだと思われます。ATTRMAPは付録Bの最後に書かれているものが正しいと思います。

data forestdata;
set forestdata00 xaxis;
run;

あとは途中のRUN;は不要ですので削除するのと、YAXISステートメントのTRANSPARENCYの値を指定すること、いくつかのタイプミスの修正が必要だと思います。

例

colorbandsattrs=(transparency=1)

まずは上記の点を修正してみて、引き続きエラーとなるようでしたらログメッセージをそのままコピーペーストして投稿してみてください😀