R/GetE.R

In JMH: Joint Model of Heterogeneous Repeated Measures and Survival Data

GetE <- function(beta, tau, gamma1, gamma2, alpha1, alpha2, vee1, vee2, H01, H02, 
                 Sig, Z, X1, W, Y, X2, survtime, cmprsk, mdata, mdataS, xsmatrix, wsmatrix) {
  
  
  n <- nrow(X2)
  p1a <- ncol(Z)

  CUH01 <- rep(0, n)
  CUH02 <- rep(0, n)
  HAZ01 <- rep(0, n)
  HAZ02 <- rep(0, n)
  
  CumuH01 <- cumsum(H01[, 3])
  CumuH02 <- cumsum(H02[, 3])
  
  getHazard(CumuH01, CumuH02, survtime, cmprsk, H01, H02, CUH01, CUH02, HAZ01, HAZ02)
  
  
  status = getEC(beta, tau, gamma1,  gamma2, alpha1, alpha2, vee1, vee2,  H01, 
                 H02, Sig, Z, X1, W, Y, X2, survtime, cmprsk, mdata, mdataS, xsmatrix, 
                 wsmatrix, CUH01, CUH02, HAZ01, HAZ02)
  
  return(status)
  
}