R/get_hatvalues.R

In spmodel: Spatial Statistical Modeling and Prediction

#' Get the (whitened) hat matrix (leverage)
#'
#' @param cov_betahat Covariance of fixed effects
#' @param cholprods A \code{cholprods} object
#'
#' @return Leverage (hat) values
#'
#' @noRd
get_hatvalues <- function(cov_betahat, SqrtSigInv_X) {
  # the hat matrix of the whitened residuals
  diag(SqrtSigInv_X %*% tcrossprod(cov_betahat, SqrtSigInv_X))
}