Stataを使用して、競合するリスクの存在下で累積発生関数をモデル化しています。各呪文はイベントの発生で終了しますが（ "失敗"）、その後に新しい呪文を入力する可能性があります。したがって、呪文は独立したものではなく、主題に入れ子になっています。Stataのstcrprepを使用して標準エラーをクラスタ化する方法
標準Stataコマンドstcrregは、サブジェクトレベルでクラスタ化された標準エラーをモデリングすることでこの構造を処理できます。しかし、私のデータセットは膨大であり（300万回以上の観測）、計算時間が膨大になります。他のユーザーは、ユーザー作成のプログラムstcrprepの使用を提案しています。stcrprepは追加の機能も備えています。
stcrprepを使用して標準エラーをクラスタ化する方法を知っている人はいますか？
ヘルプファイルの情報がstcrprepであり、コマンドの作者（http://repec.org/usug2013/lambert.uk13.pdf）のプレゼンテーションで解説されています。
クラスタ化された標準エラーは、次のように簡単に達成できます。
stcrprepを用いて生存分析のためのデータを準備する。stsetを使用して正しいフォーマットを作成します。 iウェイトとしてそれらを指定することは、堅牢でクラスターロバストな標準誤差の推定を不可能にする。stcoxを使用して競合リスク回帰を実行し、希望する,vce()オプションを指定します。このようなデータ準備によって、stcrregと全く同じ結果が得られることに注意してください。 use "http://www.stata-journal.com/software/sj4-2/st0059/prostatecancer", clear
stset time, failure(status=1,2,3) id(id) //specify all possible events
stcrprep, events(status) keep(treatment age therapy) trans(1) //trans creates failcode (event of interest)
gen event = failcode == status
stset tstop [pw=weight_c], failure(event==1) enter(tstart)
stcox treatment age, vce(cluster therapy)