I was trying to create a function to calculate leverages for some data, and I came to this:
leverages = function(x,n){
n = length(x)
leverage = array(0,length(x))
for (i in 1:length(x)){
leverage[i] = round(1/n + ((mean(x)-x[i])^2)/(n-1)*var(x), d = 4);
cat(paste("h_ii number", c(i,i), "equals", round(print(leverage[i]),d = 4)))
}
}
Here's the output, for some leverages(X,length(X)):
[1] 11.3596
h_ii number 1 equals 11.3596[1] 1.9143
h_ii number 2 equals 1.9143[1] 0.0577
How can I make this output better? Thank you!