R/genfrail.control.R

In frailtySurv: General Semiparametric Shared Frailty Model

genfrail.control <- function(censor.reltol = 1e-4,
                             censor.subdivisions = 1000L,
                             crowther.reltol = 1e-4,
                             crowther.subdivisions = 1000L
                             ) {
  
  if (censor.reltol < 0) 
    stop ("Invalid censor relative tolerance")
  if (censor.subdivisions < 0) 
    stop ("Invalid censor subdivisions")
  if (crowther.reltol < 0) 
    stop("Invalid Crowther relative tolerance")
  if (crowther.subdivisions < 0) 
    stop ("Invalid Crowther integration subdivisions")
  
  list(censor.reltol=censor.reltol,
       censor.subdivisions=censor.subdivisions,
       crowther.reltol=crowther.reltol,
       crowther.subdivisions=crowther.subdivisions
       )
}