function(x)
{
	n <- nrow(x)
	sx <- cbind(jay(n, 1), x)
	hx <- sx %*% inv(t(sx) %*% sx) %*% t(sx)
	lev <- diag(hx)
	lev
}
