data(pbc)
dta = na.omit(pbc[,c('time', 'status', 'age', "edema")])
dta = dta[dta$time>600 & dta$time<2000,]
dta[,'status'] = sign(dta[,'status'])
colnames(dta) = c('time', 'status', 'age', "edema")
res = bandselect(cva = dta[,3:4], delta = dta$status, obstime = dta$time, bandwidth = 700)
Run the code above in your browser using DataLab