R/E0vect.n.R

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

"E0vect.n" <-
function(xbar) {
#xk   <- 1.5477
#E01p <- integrate(Psiphi.n, lower=-xk,upper=xk)$value
#E0pp <- integrate(Chiphi.n, lower=-10,upper=10)$value
i1   <- 0; i2  <-  -1.337277e-5
E0   <- matrix(c(i1*xbar,i2),ncol=1)
E0}