R/cooks.dist.toureg.R

cooks.dist.toureg <- function(object, ...){
	h <- gleverage(object)   # former hatvalues
	rp <- residuals(object)
	k <- length(object$coefficients)
	cd <- (h*rp^2) / ( k * (1-h)^2 )
	return(cd)
}

Try the touchard package in your browser

Any scripts or data that you put into this service are public.

touchard documentation built on May 31, 2019, 5:04 p.m.