R/Q2vect.n.R

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

"Q2vect.n" <-
function(XtX,xbar,E.S,E0l,E1,E2,MT0,MS0,q1,beta,q5) {
#i1   <- 0 #integrate(s1phi.n, lower=-u,upper=u)$value
#i2   <- 0 #integrate(s1s2phi.n, lower=-u,upper=u)$value
Ex    <- as.matrix(xbar)
#Q2   <- Ex%*%(i2+q4*t(E1)%*%t(MT0)%*%Ex+q5*MS0%*%E1-beta*i1)
#tmp  <- E2 + q4*E.S%*%t(MT0)%*%Ex +q5*E.S%*%t(MS0) - beta*E0l
#Q2   <- Q2 + q1*XtX%*%MT0%*%tmp + q2*Ex%*%MS0%*%tmp
Q2    <- q5*Ex%*%MS0%*%E1 + q1*XtX%*%MT0%*%(E2 + q5*E.S%*%t(MS0) - beta*E0l)   
Q2}