Tertimbang bobot persegi definisi setidaknya: R fungsi lm vs

9

Adakah yang bisa memberi tahu saya mengapa saya mendapatkan hasil yang berbeda dari Rkuadrat 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. WAx=WbWAb

Saya mencoba membandingkan hasilnya dengan R lmfungsi menggunakan weightsargumen.

masukkan deskripsi gambar di sini

Haitao Du
sumber
Saya mengedit tag: ini jelas bukan [belajar sendiri]. Ini juga bukan benar-benar tentang GLS (tetapi tentang kasus yang sangat khusus), jadi saya menghapus yang juga.
amoeba

Jawaban:

13

Seperti yang dapat Anda lihat dari ekspresi matematika untuk perhitungan Anda, Anda memperolehnya

((WA)(WA))1((WA)(Wb))=(AW2A)1(AW2b).

W2W

> lm(form, mtcars, weights=w^2)
Coefficients:
      wt        hp      disp  
14.12980   0.08391  -0.16446 

Perjanjian ini sempurna (untuk kesalahan floating point - secara internal, Rmenggunakan algoritma numerik yang lebih stabil.)

whuber
sumber
1
WW2
Ax=bWAx=Wb
8
W