In celehs/DPH: Discrete Proportional Hazards Models

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(DPH)

Here is a simulated example which shows you how to use DPH:

sim_data

time <- sim_data$time
status <- sim_data$status
pred <- sim_data[, -(1:2)]

Outcome summary of the simulated data:

tbl <- table(time, status)
addmargins(tbl)

Standard DPH Model

dph(
  time = time, 
  status = status, 
  pred = pred
)

ADPH Model

Case 1

adph(
  time = time, 
  status = status, 
  pred = pred,
  sens = 0.8,
  spec = 0.9
)

Case 2

adph(
  time = time, 
  status = status, 
  pred = pred,
  sens = 0.8,
  spec = 0.9,
  sens_known = FALSE
)

Case 3

adph(
  time = time, 
  status = status, 
  pred = pred,
  sens = 0.8,
  spec = 0.9,
  spec_known = FALSE
)

Case 4

adph(
  time = time, 
  status = status, 
  pred = pred,
  sens = 0.8,
  spec = 0.9,
  sens_known = FALSE,  
  spec_known = FALSE
)

ADPH2 Model

Model Fitting

system.time(
  adph2_model <- adph2(
    time = time, 
    status = status, 
    pred = pred
  )
)
adph2_model

Perturbation Resampling

system.time(
  adph2_model <- adph2(
    time = time, 
    status = status, 
    pred = pred,
    n_ptb = 10
  )
)
adph2_model

proc.time()

celehs/DPH documentation built on April 2, 2021, 7:27 p.m.