Я пытался создать функцию для расчета рычагов для некоторых данных, и я пришел к этому:
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)))
}
}
Вот вывод для некоторых рычагов (X, длина (X)):
[1] 11.3596
h_ii number 1 equals 11.3596[1] 1.9143
h_ii number 2 equals 1.9143[1] 0.0577
Как мне улучшить этот вывод? Спасибо!