r - rmsでのcr.setup

Question

パッケージcr.setup内の関数がどのように機能するか説明してもらえますか？rms初期データを再マッピングする方法と、この再マッピングが継続率モデルにどのように役立つかを理解できないようです。ヘルプと例はそれほど役に立ちません。ネット上でも他の説明は見つかりません。

Answer 1

彼の優れたテキスト「回帰モデリング戦略」の中で、ハレルは、継続率モデルに専念する通常のロジスティック回帰に関する章の最後に3ページあります。cr.setupは、特定の行を複製し、さまざまな比較のための層マーカーを作成することによって「通常のロジスティック回帰をトリックする」プロセスをサポートしていますY >= 0; Y>=1, ... Y>=K-1。また、特定の層の「結果」を表す適切な応答変数を作成します。cr.setupの彼の最初の例を見てください。

y <- c(NA, 10, 21, 32, 32)
> cr.setup(y)
$y
[1] NA  1  0  1  0  0  0  0

$cohort
[1] <NA>  all   all   y>=21 all   y>=21 all   y>=21
Levels: all y>=21

$subs
[1] 1 2 3 3 4 4 5 5

$reps
[1] 1 1 2 2 2

NA以外のYのレベルが3つある場合、ネオアウトカムのレベルは2つだけになります。yベクトルはネオアウトカムです。ベクトル要素は、元のデータのsubsインデックスです。repsベクトルは、必要な複製の数をソフトウェアに通知します。例に従うことで、これが実際にどのように使用されているかを確認できます。

combinations <- expand.grid(cohort=levels(cohort), sex=levels(sex))
combinations