Làm thế nào là phân phối gamma nghịch đảo liên quan đến


8

Cho rằng dự toán sau của của một khả năng bình thường và một gamma nghịch đảo trước trên σ 2 là:σ2σ2

σ2IG(α+n2,β+i=1n(yiμ)22)

tương đương với

σ2IG(n2,nσ22)

kể từ khi một yếu trước trên σ 2 để loại bỏ việc alphaβ từ eqn 1:IG(α,β)σ2αβ

σ2IG(n2,i=1n(yiμ)22)

Rõ ràng là các ước tính sau của là một hàm của kích thước mẫu và tổng bình phương của các khả năng. Nhưng điều này có nghĩa là gì? Có một dẫn xuất trên Wikipedia mà tôi không hoàn toàn theo dõi.σ2

Tôi có những câu hỏi sau

  1. Tôi có thể đến phương trình thứ hai này mà không cần gọi quy tắc của Bayes không? Tôi tò mò liệu có một cái gì đó vốn có trong các tham số của một IG có liên quan đến giá trị trung bình và phương sai không phụ thuộc vào khả năng bình thường.
  2. Tôi có thể sử dụng kích thước mẫu và độ lệch chuẩn từ một nghiên cứu trước đó để ước tính thông báo trước , và sau đó cập nhật trước với dữ liệu mới không? Điều này có vẻ đơn giản, nhưng tôi không thể tìm thấy bất kỳ ví dụ nào về việc đó, hoặc lý do tại sao đây sẽ là một cách tiếp cận hợp pháp - ngoài những gì có thể nhìn thấy ở phía sau.σ2
  3. Có một xác suất phổ biến hoặc sách giáo khoa thống kê mà tôi có thể tham khảo để giải thích thêm?

Bạn không có nghĩa là một khả năng nghịch đảo gamma và một nghịch đảo gamma trước?
Neil G

1
Trước hết, tôi thấy trong câu hỏi của bạn có một số hiểu lầm: từ định lý Bayes bạn không nhận được ước tính sau, mà toàn bộ phân phối sau. Điểm thứ hai là phân phối sau này không phụ thuộc vào "tổng bình phương khả năng". Nó chỉ đơn giản phụ thuộc vào kích thước mẫu của bạn (cụ thể là n) và các giá trị mẫu, điều này là hoàn toàn tự nhiên và hợp lý. Các phụ thuộc này ảnh hưởng đến ước tính trung bình, phương sai của bạn, v.v. Ví dụ: thông số phương sai trung bình sau của bạn bằng
1n2(yiμ)2
Tomas

@thomas theo ước tính, tôi có nghĩa là ước tính phân phối sau;. Không phải tổng số bình phương ở phía sau là phép tính chính xác giống như thuật ngữ ss trong khả năng bình thường sao?
Abe

1
@ Gần đây tôi đã hỏi (và trả lời) một câu hỏi liên quan đến câu hỏi của bạn. 2. Người ta đã cho SD và SD của SD cách tính gamma tương ứng trước độ chính xác của phân phối bình thường: Câu hỏi có tại đây: stats.stackexchange.com/questions/41187/iêu
Rasmus Bååth 12/11

Câu trả lời:


9

Tôi nghĩ rằng đó là chính xác hơn để nói về sự phân bố sau của tham số của bạn chứ không phải là ước tính sau của nó. Để làm rõ các ký hiệu, tôi sẽ thả thủ trong σ ' 2 trong những gì sau.σ2σ2

Giả sử rằng được phân phối như N ( 0 , σ 2 ) , - Tôi thả μ cho bây giờ để làm một ví dụ dựa trên kinh nghiệm - và 1 / σ 2 = σXN(0,σ2)μ được phân phối nhưΓ(α,β)và không phụ thuộc vàoX.1/σ2=σ2Γ(α,β)X

Pdf của cho σ - 2 là Gaussian, nghĩa làXσ2

f(x|σ2)=12πσ2exp(x22σ2).

(X,σ2)f(x,σ2)f(x|σ2)g(σ2)σ2

f(x,σ2)=12πσ2exp(x22σ2)βαΓ(α)exp(βσ2)1σ2(α1).

Chúng ta có thể nhóm các thuật ngữ tương tự và viết lại như sau

f(x,σ2)σ2(α1/2)exp(σ2(β+x2/2)).

σ2σ2xf(x,σ2)/f(x)f(σ2|x)f(x,σ2)Γσ2f(x)

f(x)(β+x2/2)(α+1/2),

vì vậy bằng cách chia chúng ta nhận được

f(σ2|x)(β+x2/2)(σ2(β+x2/2))α1/2exp(σ2(β+x2/2))(σ2(β+x2/2))α1/2exp(σ2(β+x2/2)).

Γ(α+1/2,β+x2/2)

((x1,σ12),...,(xn,σn2))σi2f(x1,...,xn)f(σ12,...,σn2|x1,...,xn)

f(σ12,...,σn2|x1,...,xn)i=1n(σi2(β+xi2/2))α1/2exp(σi2(β+xi2/2)),

Γσi2Γ

xiσ2σ2Γ(α,β)xiσ2

f(x1,...,xn,σ2)σ2(α+n/2)exp(σ2(β+12i=1nxi2)),

σ2

σ2Γαβnσ2α/βα/β2α=βσ2bởi vì phương sai trở nên rất lớn Các giá trị là nhỏ, bạn có thể thả chúng từ các phương trình trên và bạn kết thúc với phương trình 3 của mình.

ΓS2σ2

Về câu hỏi của bạn 2. tất nhiên bạn có thể sử dụng các giá trị thu được trong một thử nghiệm trước đó với tư cách là linh mục của bạn. Vì chúng tôi đã thiết lập song song giữa giải thích Bayes và người thường xuyên ở trên, chúng tôi có thể giải thích và nói rằng nó giống như tính toán một phương sai từ một cỡ mẫu nhỏ và sau đó thu thập nhiều điểm dữ liệu hơn: bạn sẽ cập nhật ước tính của phương sai thay vì vứt bỏ các điểm dữ liệu đầu tiên.

Về câu hỏi của bạn 3. Tôi thích Giới thiệu về Thống kê toán học của Hogg, McKean và Craig, thường đưa ra chi tiết về cách rút ra các phương trình này.


1

Đối với câu hỏi 1, phương trình thứ hai tuân theo quy tắc của Bayes khi bạn chỉ ra và tôi không thấy cách nào để tránh điều đó.

Đối với câu hỏi 2, vâng, bạn có thể làm điều này. Chỉ cần sử dụng một hình thức trước giống như phương trình thứ hai của bạn.

Đối với câu hỏi 3, tôi sẽ tìm kiếm một cái gì đó về các gia đình theo cấp số nhân. Có thể ai đó sẽ giới thiệu một nguồn lực tốt.

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.