Расчет канонической функции связи в GLM


12

Я думал, что каноническая функция связи происходит от естественного параметра экспоненциального семейства. Скажем, рассмотрим семейство f ( y , θ , ψ ) = exp { y θ - b ( θ )g(⋅) тогдаθ=θ(µ)- каноническая функция связи. Возьмемв качестве примера распределение Бернулли:P(Y=y)=μy(1-μ)1-y=exp{ylogμ

f(y,θ,ψ)=exp⁡{yθ−b(θ)a(ψ)−c(y,ψ)}
θ=θ(μ) Итак, каноническая функция связиg(μ)=logμ
P(Y=y)=μy(1−μ)1−y=exp⁡{ylog⁡μ1−μ+log⁡(1−μ)}
g(μ)=log⁡μ1−μ

Но когда я вижу этот слайд , он утверждает, что Хотя это может быть легко проверено для этого конкретного распределения (и некоторых других распределений, таких как распределение Пуассона), я не вижу эквивалентности для общего случая. Кто-нибудь может дать подсказки? Спасибо ~

g′(μ)=1V(μ)

Ответы:


14

V(μ)=μ(1−μ)g(μ)=log⁡μ1−μ=log⁡μ−log⁡(1−μ)

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

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

VV(μ)−1


Спасибо @NRH. На самом деле я знаю эквивалентность для распределения Бернулли. Мне интересно общее дело. И спасибо за вашу ссылку, я проверю это :)
— ziyuang

@ziyuang, общий случай теперь включен.
— NRH

1
∫f(y,θ,ψ)dy=1θμ

Спасибо. И я нашел другую ссылку: fedc.wiwi.hu-berlin.de/xplore/ebooks/html/spm/…
— ziyuang
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.