TML.censored.control.ref: Control parameters for the refinement IRLS algorithm of the...

In RobustAFT: Truncated Maximum Likelihood Fit and Robust Accelerated Failure Time Regression for Gaussian and Log-Weibull Case

Control parameters for the refinement IRLS algorithm of the TML.censored initial S-estimates

Description

Auxiliary function for TML.censored. Typically only used internally by TML.censored, but may be used to provide a control argument. This function provides default values.

Usage

TML.censored.control.ref(maxit.sigma=2, tol.sigma=0.0001, maxit.Beta=2, 
      tol.Beta=0.0001, Maxit.S=50, tol.S.sigma=0.001, tol.S.Beta=0.001, 
      alg.sigma=1, nitmon = FALSE)

Arguments

maxit.sigma	Maximum number of iterations in scale step.
tol.sigma	Tolerance for sigma in scale step.
maxit.Beta	Maximum number of iterations in coefficient step.
tol.Beta	Tolerance for coefficients in coefficient step.
Maxit.S	Maximum number of iterations in global cycle.
tol.S.sigma	Tolerance for sigma in global cycle.
tol.S.Beta	Tolerance for coefficients in global cycle.
alg.sigma	Type of algorithm in scale step: 1: fixed point algorithm. 2: regula falsi.
nitmon	Set to TRUE if iteration monitoring is desired. Default=FALSE.

Value

A list with components named as the arguments.

See Also

TML.censored, TML.censored.control.S, TML.censored.control.tml

Examples

     ### In the example(TML.censored), the control argument for the refinement 
     ### algorithm can be built using this function:
## Not run: 
     data(MCI)
     attach(MCI)
     
     # Robust Accelerated Failure Time Regression with Gaussian errors

     ctrol.ref <- TML.censored.control.ref(maxit.sigma=2,tol.sigma=0.0001,
                  maxit.Beta=2,tol.Beta=0.0001, Maxit.S=50, tol.S.sigma=0.001, 
                  tol.S.Beta=0.001,alg.sigma=1,nitmon=FALSE)

     ctrol.tml <- list(maxit.sigma=50,tol.sigma=0.0001,maxit.Beta=50,
                  tol.Beta=0.0001, Maxit.TML=50, tol.TML.sigma=0.001, 
                  tol.TML.Beta=0.001, alg.sigma=1,nitmon=FALSE)
     
     WML<-TML.censored(log(LOS)~TypAdm*Age,data=MCI,delta=Dest,otp="adaptive",
                  control.ref=ctrol.ref,control.tml=ctrol.tml)

     summary(WML)

## End(Not run)

RobustAFT documentation built on Aug. 21, 2023, 5:13 p.m.