R/Q3sca.R

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

"Q3sca" <-
function(xbar,E.S,E0l,E2,MT0,MS0,l,u,beta,q4,q5) {
i1   <- beta^2*(pezez(u)-pezez(l))
i2   <- integrate(s2s2phi.w,lower=l,upper=u)$value
i3   <- integrate(s2phi.w,lower=l,upper=u)$value
Ex   <- as.matrix(xbar)
Q3   <- i1+i2+q4^2*t(Ex)%*%MT0%*%E.S%*%t(MT0)%*%Ex+q5^2*MS0%*%E.S%*%t(MS0)
Q3   <- Q3 - 2*beta*i3 - 2*beta*q4*t(Ex)%*%MT0%*%E0l - 2*beta*q5*MS0%*%E0l
Q3   <- Q3 + 2*q4*t(Ex)%*%MT0%*%E2 + 2*q5*MS0%*%E2 + 2*q4*q5*t(Ex)%*%MT0%*%E.S%*%t(MS0)
Q3}