R/call_estimating_equations.R

In jaredhuling/aftiv: Instrumental Variable Estimation for Time-to-Event Data

AFTScorePre.cpp <- function(beta, survival, X, ZXmat) 
{
  
  if (is.null(survival$log.t)) 
  {
    log.t <- log(survival$t)
  } else 
  {
    log.t <- survival$log.t
  }
  
  stopifnot(length(beta) == NCOL(X))
  stopifnot(NROW(X) == length(survival$delta))
  
  err       <- log.t - X %*% beta
  order.idx <- order(err)
  X         <- as.matrix(X[order.idx,, drop = FALSE]) 
  delta     <- survival$delta
  delta     <- delta[order.idx]
  
  .Call("AFTScorePre", XX = X, delta_vec = delta, PACKAGE = "aftiv")
}

AFTivScorePre.cpp <- function(beta, survival, X, ZXmat) 
{
  
  if (is.null(survival$log.t)) 
  {
    log.t <- log(survival$t)
  } else 
  {
    log.t <- survival$log.t
  }
  
  stopifnot(length(beta) == NCOL(X))
  stopifnot(NROW(X) == NROW(ZXmat))
  stopifnot(NROW(X) == length(survival$delta))
  
  err       <- log.t - X %*% beta
  order.idx <- order(err)
  ZXmat     <- as.matrix(ZXmat[order.idx,, drop = FALSE]) 
  delta     <- survival$delta
  delta     <- delta[order.idx]
  
  .Call("AFTScorePre", XX = ZXmat, delta_vec = delta, PACKAGE = "aftiv")
}


AFTivIPCWScorePre.cpp <- function(beta, survival, X, ZXmat, GC) 
{
  
  if (is.null(survival$log.t)) 
  {
    log.t <- log(survival$t)
  } else 
  {
    log.t <- survival$log.t
  }
  
  
  stopifnot(length(beta) == NCOL(X))
  stopifnot(NROW(X) == NROW(ZXmat))
  stopifnot(NROW(X) == length(survival$delta))
  
  err       <- log.t - X %*% beta
  order.idx <- order(err)
  ZXmat     <- as.matrix(ZXmat[order.idx,, drop = FALSE]) 
  delta     <- survival$delta
  delta     <- delta[order.idx]
  
  #.Call("AFTivIPCWScorePre", XX = ZXmat, 
  #                           delta_vec = delta, 
  #                           GCT_ = GCT, 
  #                           icpw_num_denom_ = icpw.num.denom,
  #                           conf_x_loc_ = conf.x.loc,
  #      PACKAGE = "aftiv")
}