bs.hqm {HQM} | R Documentation |
Calculates the Brier score for the HQM estimator.
bs.hqm(xin, est, landm, th, event_time_name, status_name)
xin |
A data frame containing event times and the patient status. |
est |
The HQM estimator values, typically the output of |
landm |
Landmark time. |
th |
Time horizon. |
event_time_name |
The column name of the event times in the data frame |
status_name |
The column name of the status variable in the data frame |
The function bs.hqm
implements the Brier score calculation for the HQM estimator estimator.
Scalar: the Brier score of the HQM estimator.
library(pec)
library(survival)
Landmark <- 2
pbcT1 <- pbc2[which(pbc2$year< Landmark & pbc2$years> Landmark),]
timesS2 <- seq(Landmark,14,by=0.5)
b=0.9
arg1<- get_h_x(pbcT1, 'albumin', event_time_name = 'years',
time_name = 'year', event_name = 'status2', 2, 0.9)
br_s2 = seq(Landmark, 14, length=99)
sfalb2<- make_sf( (br_s2[2]-br_s2[1])/4 , arg1)
tHor <- 1.5
bs.use<-bs.hqm(pbcT1, sfalb2, Landmark,tHor,
event_time_name = 'years', status_name = 'status2')
bs.use