### In the example(TML.censored), the control argument for the refinement
### algorithm can be built using this function:
if (FALSE) {
data(MCI)
attach(MCI)
# Robust Accelerated Failure Time Regression with Gaussian errors
ctrol.S <- list(N=150, q=5, sigma0=1, MAXIT=100, TOL=0.001,seed=123)
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.S=ctrol.S,control.ref=ctrol.ref,
control.tml=ctrol.tml)
summary(WML)
}
Run the code above in your browser using DataLab