Nothing
one <- rep(1,4)
X <- cbind(one,x); X
x <- as.matrix(x);x
y <- as.matrix(y); y
v0 <- as.matrix(one); v0
u0 <- v0/sqrt(sum(v0^2)); u0
v1 <- x - mean(x); v1
u1 <- v1/sqrt(sum(v1^2)); u1
#
# the hat matrix
H <- u0 %*% t(u0) + u1 %*% t(u1)
H
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.