Cálculo de la función de enlace canónico en GLM

12

Pensé que la función de enlace canónico proviene del parámetro natural de la familia exponencial. Digamos, considere la familia then \ theta = \ theta (\ mu) es la función de enlace canónico. Tome la distribución de Bernoulli como ejemplo, tenemos P (Y = y) = \ mu ^ {y} (1- \ mu) ^ {1-y} = \ exp \ left \ {y \ log \ frac {\ mu} {1- \ mu} + \ log {(1- \ mu)} \ right \} Entonces, la función de enlace canónico g (\ mu) = \ log \ frac {\ mu} {1- \ mu}g()

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

Pero cuando veo esta diapositiva , afirma que

g(μ)=1V(μ)
Aunque se puede verificar fácilmente para esta distribución en particular (y algunas otras distribuciones, como la distribución de Poisson), No puedo ver la equivalencia para el caso general. ¿Alguien puede dar pistas? Gracias ~
ziyuang
fuente

Respuestas:

14

La función de varianza para la variable Bernoulli es V(μ)=μ(1μ) . Verificamos fácilmente que con el enlace canónico g(μ)=logμ1μ=logμlog(1μ) luego

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

Para el caso general, se deriva de la definición de que ver p. ej. páginas 28-29 en McCullagh y Nelder . Con el enlace canónico tenemos , y la función de varianza se define como , que en términos de convierte en Por diferenciación de la identidad obtenemos

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

En la construcción de funciones cuasi-verosimilitud es natural comenzar con la relación entre la media y la varianza, dada en términos de la función de la varianza . En este contexto, la anti-derivada de puede interpretarse como una generalización de la función de enlace; véase, por ejemplo, la definición de cuasi-verosimilitud (log) en la página 325 (fórmula 9.3 ) en McCullagh y Nelder . VV(μ)1

NRH
fuente
Gracias @NRH En realidad, sé la equivalencia para la distribución de Bernoulli. Me pregunto el caso general. Y gracias por tu referencia, lo comprobaré :)
ziyuang
@ziyuang, el caso general ahora está incluido.
NRH
1
@NRH: solo para agregar a esta respuesta, las fórmulas de media y varianza se pueden derivar diferenciando la ecuación en ambos lados con respecto a (o equivalentemente ) La primera derivada te da la media, la segunda te da la varianza. f(y,θ,ψ)dy=1θμ
probabilidadislogica
Gracias. Y he encontrado otro enlace de referencia: fedc.wiwi.hu-berlin.de/xplore/ebooks/html/spm/…
ziyuang