Displaying 1 result from an estimated 1 matches for "pred_score".
Did you mean:
pre_store
2012 Jan 17
1
Scoring using cox model: probability of survival before time t
..., close the
earlier session and run the below script in the new #R session, it gives
error.
library(survival)
library(peperr)
load(file = file.path("C:/Desktop","fit_coxph.RData"))
n = 1000
set.seed(1)
x1 = rnorm(n,0)
x2 = rnorm(n,0)
score_data <- data.frame(x1,x2)
pred_score <- predictProb.coxph(fit_coxph, Surv(time, event), score_data,
0.04)
#Error in Surv(time, event) : Time variable is not numeric
#After creating dummy place holder for Surv(time, event), it gives another
error:
time <- rep(2, n)
event <- rep(1, n)
pred_score <- predictProb.coxph(fit...