HQM (version 1.0)

bs.hqm: Brier score for the High Quality Marker estimator

Description

Calculates the Brier score for the HQM estimator.

Usage

bs.hqm(xin, est, landm, th, event_time_name, status_name)

Value

Scalar: the Brier score of the HQM estimator.

Arguments

xin

A data frame containing event times and the patient status.

est

The HQM estimator values, typically the output of get_h_x.

landm

Landmark time.

th

Time horizon.

event_time_name

The column name of the event times in the data frame xin.

status_name

The column name of the status variable in the data frame xin.

Details

The function bs.hqm implements the Brier score calculation for the HQM estimator estimator.

See Also

auc.hqm

Examples

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