使用Coxモデルは、私が最初にRにCoxモデルを確立し、生存

Question

を推定する：

test1<- test[1:20,] 
model.1 <- coxph(Surv(test1$days,test1$status==1) ~ test1$MTT+test1$ADC,data=test1) 

をし、私はこのような次の患者の生存を予測しようとしたとき：

covs1 <- data.frame(test[21,]$MTT,test[21,]$ADC) 
summary(survfit(model.1, newdata= covs1, type ="aalen")) 

は、それは私にあまりにも多くの生存率を与えました結果と警告は "newdataは1行ありましたが、変数には20行あり" fyi、イベントは20件あり、結果には20の生存結果が含まれています。

Answer 1

予測の基礎として指定されているデータフレーム内の列の名前は、モデル式のRHSと同じ列名でなければなりません。

test1<- test[1:20,] 
model.1 <- coxph(Surv(days, status==1) ~ MTT + ADC, data=test1) 

covs1 <- test[21, c("MTT", "ADC")] 
# then do your prediction 

あなたはSurvに引数を指定する$を使用しないでください。私はあなたがこのような何かをしない限り、あなたがqualifiyうとは思いません。モデルがデータフレームの環境で構築されることが重要です。