Adakah yang bisa memberi tahu saya mengapa saya mendapatkan hasil yang berbeda dari R
kuadrat terkecil tertimbang dan solusi manual dengan operasi matriks ?
Secara khusus, saya mencoba untuk secara manual menyelesaikan , di mana W adalah matriks diagonal pada bobot, A adalah matriks data, b adalah vektor respons.
Saya mencoba membandingkan hasilnya dengan R
lm
fungsi menggunakan weights
argumen.
Jawaban:
Seperti yang dapat Anda lihat dari ekspresi matematika untuk perhitungan Anda, Anda memperolehnya
Perjanjian ini sempurna (untuk kesalahan floating point - secara internal,
R
menggunakan algoritma numerik yang lebih stabil.)sumber