# Script for computing the hat matrix H for
# the CrudeOil dataset.
crude = read.table("c:/datasets/crude-oil.txt", header=TRUE)

# Create datasets
X = cbind(1, as.matrix(crude[ ,1:2])
y = as.matrix(crude[ , 3])

# Compute hat matrix.
H = X %*% solve(t(X) %*% X) %*% t(X)
print(H)