Розрахунок канонічної функції зв'язку в 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(μ)

Відповіді:


14

Функцією дисперсії для змінної Бернуллі є . Ми легко перевіримо, що за допомогою канонічного посилання то V(μ)=μ(1μ)g(μ)=logμ1μ=logμlog(1μ)

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

У загальному випадку випливає з визначення, що див., наприклад, сторінку 28-29 у McCullagh та Nelder . З канонічного зв’язку маємо , а функцію дисперсії визначаємо як , яка з точки зору стає Диференціюючи тотожність отримуємо

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

При побудові функцій квазі-правдоподібності природно почати з співвідношенням між середнім і дисперсією, заданої в термінах функції дисперсії . У цьому контексті анти-похідне від можна інтерпретувати як узагальнення функції зв'язку, див., Наприклад, визначення квазіімовірності (log) на сторінці 325 (формула 9.3 ) у McCullagh та Nelder . VV(μ)1


Дякую @NRH Насправді я знаю еквівалентність розподілу Бернуллі. Мені цікаво загальний випадок. І дякую за вашу довідку, я перевірю це :)
ziyuang

@ziyuang, тепер загальна справа включена.
НРХ

1
@NRH - просто для додання цієї відповіді середні та дисперсійні формули можна отримати, диференціюючи рівняння з обох сторін щодо (або еквівалентно ). Перша похідна дає середнє значення, друга - дисперсія. 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.