### In the example(TML.censored), the control argument for the final estimates
### can be built using this function:
if (FALSE) {
data(MCI)
attach(MCI)
# Robust Accelerated Failure Time Regression with Gaussian errors
ctrol.ref <- list(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 <- TML.censored.control.tml(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)
}
Run the code above in your browser using DataLab