おはようございます、
ケース/コントロール スタディのマッチング後に R で層化ロジスティック回帰を実行しようとしましたが、予想外の間違いに遭遇しました。mtcars データ セットを使用してエラーを再現しました。
test=mtcars
test$am=as.factor(test$am)
test$cyl=as.factor(test$cyl)
test$vs=as.factor(test$vs)
library(survival)
clogit(am~vs+strata(cyl),data=test)
coxph(formula = Surv(rep(1, 32L), am) ~ vs + strata(cyl), data = test, のエラー: Cox モデルは "mright" 生存データをサポートしていません
私の理解では、clogit-function は時間パラメーターを作成しますが、R には問題があるようです。am を因子として使用しないと、間違いはなくなりますが、それが全体ですが、ロジスティック回帰を行うために因子としてラベルを付ける必要はありませんか?
ちなみに私はR3.2.2とサバイバルパッケージ2.41-3を使っていますが、どちらもサバイバルパッケージの機能のようですので、これが根本的な問題ではないでしょうか？それとも、この間違いは新しい R バージョンでは再現できませんか?