R/Q1mat.n.R

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

"Q1mat.n" <-
function(XtX,xbar,E.S,E1,MT0,MS0,u,q1) {
i1    <- integrate(s1s1phi.n, lower=-u,upper=u)$value
Ex    <- as.matrix(xbar)
Q1    <- i1*XtX + q1*Ex%*%t(E1)%*%t(MT0)%*%XtX #+ q2*Ex%*%t(E1)%*%t(MS0)%*%t(Ex)
tmp   <- E1%*%t(Ex) + q1*E.S%*%t(MT0)%*%XtX    #+ q2*E.S%*%t(MS0)%*%t(Ex)
Q1    <- Q1 + q1*XtX%*%MT0%*%tmp
#tmp  <- E1%*%t(Ex) + q1*E.S%*%t(MT0)%*%XtX + q2*E.S%*%t(MS0)%*%t(Ex)
#Q1   <- Q1 + q2*Ex%*%MS0%*%tmp
Q1}