การคำนวณฟังก์ชั่นลิงก์แบบบัญญัติใน GLM


12

ฉันคิดว่าฟังก์ชั่น canonical linkมาจากพารามิเตอร์ธรรมชาติของตระกูล exponential พูดดูครอบครัว ดังนั้นคือฟังก์ชันลิงก์แบบบัญญัติ ใช้การกระจาย Bernoulliเป็นตัวอย่างเรามี ดังนั้นฟังก์ชันลิงก์แบบบัญญัติ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

ฟังก์ชั่นความแปรปรวนตัวแปร Bernoulli เป็นV(μ)=μ(1μ)MU) เราตรวจสอบได้ง่ายว่าด้วยการเชื่อมโยงที่ยอมรับได้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(μ),

ในการก่อสร้างของฟังก์ชั่นกึ่งโอกาสมันเป็นธรรมชาติที่จะเริ่มต้นด้วยความสัมพันธ์ระหว่างค่าเฉลี่ยและความแปรปรวนที่ได้รับในแง่ของฟังก์ชั่นความแปรปรวนVในบริบทนี้การต่อต้านอนุพันธ์ของสามารถตีความได้ว่าเป็นลักษณะทั่วไปของฟังก์ชั่นการเชื่อมโยงดูตัวอย่างคำจำกัดความของโอกาสเสมือน (บันทึก) ในหน้า 325 (สูตร 9.3) ) ในMcCullagh และ Nelder VV(μ)1


ขอบคุณ @NRH ที่จริงฉันรู้ถึงความเท่าเทียมกันของการแจกแจงเบอร์นูลลี ฉันสงสัยกรณีทั่วไป และขอขอบคุณสำหรับการอ้างอิงของคุณฉันจะตรวจสอบ :)
ziyuang

@ziyuang กรณีทั่วไปรวมอยู่ในขณะนี้
NRH

1
@NRH - เพียงเพิ่มคำตอบนี้สูตรค่าเฉลี่ยและความแปรปรวนสามารถหาได้จากการหาสมการทั้งสองด้านด้วยความเคารพ (หรือเทียบเท่า ) อนุพันธ์อันดับแรกให้ค่าเฉลี่ยกับคุณอันดับที่สองให้ความแปรปรวน f(y,θ,ψ)dy=1θμ
ความน่าจะเป็นทางการ

ขอบคุณ. และฉันได้พบลิงค์อ้างอิงอื่น: fedc.wiwi.hu-berlin.de/xplore/ebooks/html/spm/…
ziyuang
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.