R/residLCM.R

In TimeWz667/LeeCarteR: A Likelihood-based Lee Carter Modelling Tool

#' Title
#'
#' @param obj
#' @param type
#'
#' @return
#' @export
#'
#' @examples
residuals.LCM <- function(obj, type=c("deviance", "pearson")) {
  type <- match.arg(type)
  dhat <- obj$model$fitted
  ddat <- obj$data$Event/obj$data$Population

  if (obj$link=="log") {
    if (type=="deviance") {
      rd <- sign(ddat-dhat)*sqrt(ddat*(log(ddat)-log(dhat)) - (ddat-dhat))
    } else {
      rd <- (log(ddat) - log(dhat))*sqrt(ddat)
    }
  } else {
    if (type=="deviance") {
      # todo
    } else {
      # todo
    }
  }

  rd[is.na(rd)] <- 0
  rd[is.infinite(rd)] <- 0
  dimnames(rd) <- dimnames(ddat)

  rd
}