R/FD_lexc_d_alphak.R

In curesurv: Mixture and Non Mixture Parametric Cure Models to Estimate Cure Indicators

FD_lexc_d_alphak <- function(z_tau, z_alpha, x, theta) {
  n_z_tau <- ncol(z_tau)
  n_z_tau_ad <- n_z_tau - 1
  n_z_alpha <- ncol(z_alpha)
  n_z_alpha_ad <- n_z_alpha - 1
  alpha0 <- theta[1]
  if (n_z_tau > 0) {
    alpha_k <- theta[2:(n_z_alpha + 1)]
    beta <- (theta[n_z_alpha + 2])
    tau0 <- theta[n_z_tau + 3]
    tau_z <- theta[(n_z_tau + 4):(n_z_tau + 4 + n_z_tau_ad)]
    alpha <- exp(alpha0 + z_alpha %*% alpha_k)
    tau <- exp(tau0 + z_tau %*% tau_z)
    u <- x / (tau)

  }else {
    beta <- (theta[2])
    tau0 <- theta[3]
    alpha <- exp(alpha0)
    tau <- exp(tau0)
    u <- x / (tau)
  }
  uprime <- ifelse((u <= 1 & u >= 0), u, 0.99999999)

  res <- (z_alpha*c((uprime ^ ((alpha) - 1) * (1 - uprime) ^ (exp(beta - 1)))*(log(uprime)* (alpha - 1))))


  return(res)
}