GLM中规范链接函数的计算


12

我认为规范链接函数来自指数族的自然参数。说,考虑族 则是规范的链接函数。以伯努利分布为例,我们有 因此,规范链接函数g()

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

但是,当我看到这张幻灯片时,它声称 尽管可以很容易地针对此特定分布(以及其他一些分布,例如泊松分布)进行验证,我看不到一般情况的等效性。谁能给出提示?谢谢〜

g(μ)=1V(μ)

Answers:


14

Bernoulli变量的方差函数为。我们可以通过规范链接轻松地检查然后 V(μ)=μ(1μ)g(μ)=logμ1μ=logμlog(1μ)

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

对于一般情况,可以从以下定义得出: 请参见例如McCullagh和Nelder中的第28-29页。使用的规范链接,我们有,并且方差函数定义为,根据变为 通过标识微分,我们得到

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

拟似然函数的构造中,自然的是从均值和方差之间的关系开始,以方差函数。在这种情况下,的反导数可以解释为链接函数的一般化,例如,请参见第325页上的(log)拟似然的定义(公式9.3 )在McCullagh和Nelder中VV(μ)1


谢谢@NRH。实际上,我知道伯努利分布的等价性。我想知道一般情况。感谢您的参考,我将对其进行检查:)
ziyuang 2011年

@ziyuang,现在包括一般情况。
NRH

1
@NRH-只需添加此答案,就可以通过相对于(或等价在两侧微分方程来得出均值和方差公式)。一阶导数为您提供平均值,二阶为您提供方差。f(y,θ,ψ)dy=1θμ
概率

谢谢。而且我发现了另一个参考链接:fedc.wiwi.hu-berlin.de/xplore/ebooks/html/spm/...
ziyuang
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.