ループにフィットしたcoxphモデルの100個の複製があります。データ フレームまたはリスト内の各レプリケートの p 値を使用して、ログランク スコア テスト結果を抽出しようとしています。私は以下を使用しています。ただし、p 値ではなく、ログ ランク スコアのみが表示されます。どんな助けでも大歓迎です。
データセットを共有できますが、ここに添付する方法がわかりません。
ありがとう、クリナ
Repl_List <- unique(dat3$Repl)
doLogRank = function(sel_name) {
dum <- dat3[dat3$Repl == sel_name,]
reg <- with(dum, coxph(Surv(TIME_day, STATUS) ~ Treatment, ties = "breslow"))
LogRank <- with(reg, reg$score)
}
LogRank <- t(as.data.frame(lapply(Repl_List, doLogRank)))