R/ResExpW.r

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

ResExpW <-
function(r) {zero <- 1e-6
nr <- length(r); num <- rep(0,nr)
for (i in 1:nr) num[i] <- intxlwbl(r[i])
den <- 1-plweibul(r)
ind <- (den < zero)
val <- num/den; val[ind] <- r[ind]
val}