Grayの検定 グラフ下に number at risk 表を 書く方法について

[daicon]

お世話になっています。Sas Studio 3.8　および　SAS9.4使用しております
ノンパラメトリックな生存時間を用いて
ある疾病の累積発生率について　時間（TIME)　打ち切り（０）発症（１）死亡（２）として
死亡(2)を競合リスク、　関心のあるイベントを発症(1)とし
Grayの検定で計算しています。
グラフ下に number at risk 表を　書く方法について　教えて頂けますでしょうか？
どうかよろしくお願いします

proc lifetest data=WORK.TEST plots=(cif(test));
time TIME*Event_012(0) / eventcode=1;
strata Group;

[yutaka]

たぶん　通常のカプランマイヤ曲線と違って，plotステートメントのcifにはまだ現時点でnumber at riskを出す機能はないと認識してます
そのためoutcif= でだしてからちょっと加工して，SGPLOTなどで頑張らないといけないと思われます
ただ，同様の問題にあたって，一連の処理についてマクロ化して，githubにコードも含めて公開されているものがありました
検証はしてませんが，Paperとコードを見ていただいて，利用・カスタマイズするのはどうかなと思いました

・paper
https://www.hoajonline.com/journals/pdf/2053-7662-7-2.pdf

・program
https://github.com/chaozhangEmory/CompetingRiskSAS
Plots_CIF V1.sas

[daicon]

ご教示ありがとうございます。早速paper見て試みてみたいと思います。

また分からない点が出てきた際は　ご教示いただけると幸いです。