R/hvalues.toureg.R

In touchard: Touchard Model and Regression

hvalues.toureg <- function(object, ...){
  w <- object$var
  W <- diag(w)
  k <- length(object$coefficients)
  V <- object$vcov[1:k,1:k]
  
  terms <- delete.response(terms(object))
  x <- if (is.null(object$x)){
	   model.matrix(terms, model.frame(object))
	   }else{
		    object$x
		}
    
  h <- W %*% x %*% V %*% t(x) 
  return(diag(h))
}