dados<-read.table("d:\\Dados\\combustível.txt", header=TRUE)

attach(dados)

X <- cbind(1, X1, X2, X3, X4)

p <- ncol(X) - 1    # número de covariáveis adicionadas

n <- nrow(X)

t(X) %*% X

beta.chapeu = solve( t(X) %*% X) %*% t(X) %*% Y

## Alavanca

H = X %*% solve( t(X) %*% X) %*% t(X)

plot(diag(H), pch=16)

abline(2*(p+1)/n, 0)

identify(diag(H))