R/CIF2.CR.R

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

CIF2.CR <- function(data, H01, H02, s, u, bl) {
  a <- nrow(H01)
  b <- nrow(H02)
  CH02 <- cumsum(H02[, 3])
  CH01 <- rep(0, b)
  count <- 1
  i = 1
  while (count <= a & i <= b) {
    if (H01[count, 1] < H02[i, 1]) {
      CH01[i] <- CH01[i] + H01[count, 3]
      count <- count + 1
    } else {
      i <- i + 1
    }
  }
  CH01 <- cumsum(CH01)
  
  CIF2 <- 0
  p1a <- nrow(data$Sig) - 1
  bb <- as.vector(bl[1:p1a])
  w <- bl[p1a+1]
  for (i in 1:b) {
    if (s < H02[i, 1] & u >= H02[i, 1]) {
      if (i >= 2) {
        CIF2 <- CIF2 + H02[i, 3]*exp(data$X2%*%data$gamma2 + data$alpha2%*%bb + data$nu2*w)*
          exp(-CH01[i-1]*exp(data$X2%*%data$gamma1 + data$alpha1%*%bb + data$nu1*w)-
                CH02[i-1]*exp(data$X2%*%data$gamma2 + data$alpha2%*%bb + data$nu2*w))
      } else {
        CIF2 <- CIF2 + H02[i, 3]*exp(data$X2%*%data$gamma2 + data$alpha2%*%bb + data$nu2*w)
      }

      
    } else next
  }
  CIF2
}