Czy ktoś mógłby mi powiedzieć, dlaczego otrzymuję różne wyniki od R
ważonego najmniejszego kwadratu i ręcznego rozwiązania przez operację macierzy ?
W szczególności próbuję ręcznie rozwiązać , gdzie to macierz diagonalna na wagach, to macierz danych, to odpowiedź wektor.
Próbuję porównać wyniki z R
lm
funkcją za pomocą weights
argumentu.
Odpowiedzi:
Jak widać z wyrażeń matematycznych do obliczeń, otrzymujesz
Najwyraźniej twoje ciężary sąW.2) , nie W. . Dlatego powinieneś porównać swoją odpowiedź z wynikami
Zgoda jest idealna (w ramach błędu zmiennoprzecinkowego - wewnętrznie
R
wykorzystuje bardziej stabilny numerycznie algorytm.)źródło