私は 2 つの独立したデータセットを持っています。私がやりたいことは、外部データセット (test_dat、n=1000) を使用してメイン データセット (main_dat、n=5421) から取得した Cox モデルを検証することです。しかし、Hmisc パッケージの rcorr.cens() と rms の val.surv の両方を使用すると、エラー メッセージが表示されます。これが私がやっていることです：
library(rms)
surv.obj=with(main_dat,Surv(survival,surv_cens)) ## to use with rcorr.cens
phmodel=cph(surv.obj~sex+age+treatment, x=TRUE, y=TRUE, surv=T, time.inc=10, data=main_dat, se.fit=T)
estimates=survest(phmodel, newdata=test_dat, times=10)
rcorr.cens(x=estimates, S=surv.obj)
rcorr.cens(x = 推定値、S = surv.obj) のエラー: y は x と同じ長さでなければなりません
w=val.surv(phmodel ,newdata=test_dat, u=10)
val.surv(phmodel, newdata = test_dat, u = 10) のエラー: ディム [製品 1000] はオブジェクト [5421] の長さと一致しません さらに: 警告メッセージ: est.surv + S[, 1] : 長いオブジェクトの長さが短いオブジェクトの長さの倍数ではありません
何か間違ったことをしていますか、それとも 2 つのデータセットの観測数が同じでなければなりませんか?
どんな助けでも大歓迎です。