生存パッケージからのclogit回帰の強固な標準エラーR

Question

Rのサバイバルパッケージからclogit回帰の堅牢な標準エラーを取得しようとしています。その際、Stataによって報告された標準エラーを再現しようとしていますclogitコマンドはvce(robust)オプション付きです。

Rでの私の式は、関数にrobust = TRUE引数を追加

conditional_logit <- clogit(dependent_variable ~ independent_variable + some_controls + strata(year), method= "exact", data = data_frame) 

がエラーで失敗している：

Error in residuals.coxph(fit2, type = "dfbeta", weighted = TRUE) : 
score residuals are not available for the exact method 

提案されるようにサンドイッチやPLMパッケージを経由して堅牢な標準誤差を抽出しようとする試みをhere,here,here、およびhereは同じエラーで失敗します。同様に、clogit functionには、exactメソッド（44行目）を使用しているときの堅牢な標準エラーの計算を停止する条件が含まれています。ただし、clogit回帰オブジェクトには、conditional_logit $ residualsおよびconditional_logit $ scoreが存在します。

誰かが次の質問に答えることができれば、私は感謝するでしょう：

• は、「正確な」条件付きロジスティック回帰のための堅牢な標準誤差を計算すること、一般的に不可能または「間違っている」となっていますか？もしそうなら、Stataはなぜそうすることを許可しますか？
• そうでない場合：Rのclogit regressionの堅牢な標準誤差をどのように計算できますか？
• clogit回帰オブジェクトのデータに基づいて堅牢な標準誤差を計算することができない場合：サバイバルパッケージのclogit関数によって生成された条件付きロジスティック回帰モデルを生成する別のRパッケージがありますか？堅牢な標準エラーを計算するために必要なデータ

Answer 1

オリジナルのclogit関数のタイを扱う方法を変更します。

methodコマンドのデフォルトはmethod="exact"です。たとえば、method="efron"を使用すると、robust=TRUEコマンドが機能します。