Perhitungan fungsi tautan kanonik dalam GLM

12

Saya pikir fungsi tautan kanonik berasal dari parameter alami keluarga eksponensial. Katakan, pertimbangkan keluarga lalu adalah fungsi tautan kanonik. Ambil distribusi Bernoulli sebagai contoh, kami memiliki Jadi, fungsi tautan kanonikg()

f(y,θ,ψ)=exp{yθb(θ)a(ψ)c(y,ψ)}
θ=θ(μ)
P(Y=y)=μy(1μ)1y=exp{ylogμ1μ+log(1μ)}
g(μ)=logμ1μ

Tetapi ketika saya melihat slide ini , ia mengklaim bahwa Meskipun dapat dengan mudah diverifikasi untuk distribusi khusus ini (dan beberapa distribusi lainnya, seperti distribusi Poisson), Saya tidak bisa melihat kesetaraan untuk kasus umum. Adakah yang bisa memberi petunjuk? Terima kasih ~

g(μ)=1V(μ)
ziyuang
sumber

Jawaban:

14

Fungsi varians untuk variabel Bernoulli adalah . Kami dengan mudah memeriksa bahwa dengan tautan kanonik lalu V(μ)=μ(1μ)g(μ)=logμ1μ=logμlog(1μ)

g(μ)=1μ+11μ=1μ+μμ(1μ)=1μ(1μ)=1V(μ).

Untuk kasus umum, seseorang berasal dari definisi bahwa lihat misalnya halaman 28-29 di McCullagh dan Nelder . Dengan tautan kanonik kita memiliki , dan fungsi varians didefinisikan sebagai , yang dalam istilah menjadi Dengan membedakan identitas kita mendapatkan

E(Y)=μ=b(θ) and Var(Y)=b(θ)a(ψ),
gθ=g(μ)=g(b(θ))b(θ)μ
V(μ)=b(g(μ)).
θ=g(b(θ))
1=g(b(θ))b(θ)=g(μ)V(μ),

Dalam pembangunan fungsi kuasi-kemungkinan itu wajar untuk memulai dengan hubungan antara mean dan varians, diberikan dalam hal fungsi varian . Dalam konteks ini anti-turunan dari dapat diartikan sebagai generalisasi dari fungsi tautan, lihat, misalnya, definisi kemungkinan kuasi (log) pada halaman 325 (rumus 9.3) ) di McCullagh dan Nelder . VV(μ)1

NRH
sumber
@NRH terima kasih. Sebenarnya saya tahu kesetaraan untuk distribusi Bernoulli. Saya ingin tahu kasus umum. Dan terima kasih atas referensi Anda, saya akan memeriksanya :)
ziyuang
@ziyuang, kasus umum sekarang sudah termasuk.
NRH
1
@NRH - hanya untuk menambah jawaban ini, rumus mean dan varians dapat diturunkan dengan membedakan persamaan di kedua sisi sehubungan dengan (atau setara dengan ). Derivatif pertama memberi Anda mean, kedua memberi Anda varians. f(y,θ,ψ)dy=1θμ
probabilityislogic
Terima kasih. Dan saya telah menemukan tautan referensi lain: fedc.wiwi.hu-berlin.de/xplore/ebooks/html/spm/…
ziyuang