Как рассчитать дисперсию оценки наименьших квадратов с помощью QR-разложения в R?
Я пытаюсь изучить QR-декомпозицию, но не могу понять, как получить дисперсию beta_hat, не прибегая к традиционным матричным вычислениям. Я тренируюсь сiris набор данных, и вот что у меня так далеко:
y<-(iris$Sepal.Length)
x<-(iris$Sepal.Width)
X<-cbind(1,x)
n<-nrow(X)
p<-ncol(X)
qr.X<-qr(X)
b<-(t(qr.Q(qr.X)) %*% y)[1:p]
R<-qr.R(qr.X)
beta<-as.vector(backsolve(R,b))
res<-as.vector(y-X %*% beta)
Спасибо за вашу помощь!