Làm thế nào một trước không thích hợp có thể dẫn đến một phân phối sau thích hợp?


22

Chúng tôi biết rằng trong trường hợp phân phối trước thích hợp,

P(θX)=P(Xθ)P(θ)P(X)

P(Xθ)P(θ) .

Sự biện minh thông thường cho bước này là phân phối biên của , , không đổi đối với và do đó có thể bị bỏ qua khi lấy phân phối sau.P ( X )XP(X)θ

Tuy nhiên, trong trường hợp không đúng trước, làm thế nào để bạn biết rằng phân phối sau thực sự tồn tại? Dường như có một cái gì đó còn thiếu trong cuộc tranh luận dường như tròn này. Nói cách khác, nếu tôi giả sử hậu sinh tồn tại, tôi hiểu cơ chế làm thế nào để lấy được hậu thế, nhưng dường như tôi đang thiếu lý do biện minh cho lý do tại sao nó tồn tại.

Tái bút Tôi cũng nhận ra rằng có những trường hợp trước đó không đúng dẫn đến hậu thế không đúng.

Câu trả lời:


16

Chúng tôi thường chấp nhận hậu thế từ các linh mục không phù hợp nếu tồn tại và là phân phối xác suất hợp lệ (nghĩa là, nó tích hợp chính xác đến 1 trên hỗ trợ). Về cơ bản, điều này sôi xuống là hữu hạn. Nếu đây là trường hợp, thì chúng ta gọi số lượng này là và chấp nhận nó là phân phối sau mà chúng ta muốn. Tuy nhiên, điều quan trọng cần lưu ý rằng đây KHÔNG phải là phân phối sau, cũng không phải là phân phối xác suất có điều kiện (hai thuật ngữ này đồng nghĩa trong ngữ cảnh ở đây).π(θ)

π(X|θ)π(θ)π(X)
π(X)= =π(X|θ)π(θ)dθπ(θ|X)

Bây giờ, tôi đã nói rằng chúng tôi chấp nhận phân phối 'hậu thế' từ các linh mục không phù hợp được đưa ra ở trên. Lý do chúng được chấp nhận là vì vẫn sẽ cung cấp cho chúng tôi 'điểm số' tương đối trên không gian tham số; tức là tỷ lệ mang lại ý nghĩa cho phân tích của chúng tôi. Ý nghĩa chúng ta nhận được từ các linh mục không phù hợp trong một số trường hợp có thể không có sẵn trong các linh mục thích hợp. Đây là một biện minh tiềm năng cho việc sử dụng chúng. Xem câu trả lời của Sergio để kiểm tra kỹ hơn về động lực thực tế cho các linh mục không đúng.π ( θ 1 )π(θ)π(θ1)π(θ2)

Điều đáng chú ý là số lượng này cũng có các đặc tính lý thuyết mong muốn, Degroot & Schervish :π(θ|X)

Các linh mục không đúng không phải là phân phối xác suất thực sự, nhưng nếu chúng ta giả vờ như vậy, chúng ta sẽ tính toán các phân phối sau gần đúng với các hậu thế mà chúng ta sẽ có được bằng cách sử dụng các linh mục liên hợp thích hợp với các giá trị cực trị của siêu âm trước.


Tôi bối rối bởi một vài điều trong câu trả lời của bạn. Bạn nói rằng chúng tôi chấp nhận hậu thế nếu những điều trên là hữu hạn. Điều đó có nghĩa là nếu tích phân đó không hữu hạn thì hậu thế sẽ không hữu hạn? Ngoài ra, bạn dường như ngụ ý rằng chúng tôi sử dụng hậu thế trong trường hợp này, nhưng nó không phải là một bản phân phối thực sự - đúng không? không có trường hợp nó là một phân phối thực sự? Ngoài ra, tỷ lệ của các linh mục phải làm gì với điều này? Tôi không thấy kết nối.
Phường Ben Elizabeth

@BenElizabethWard Nếu tồn tại, thì tích phân phải tồn tại (và do đó là hữu hạn). Giao tiếp cũng đúng: nếu không tồn tại (là vô hạn), thì không tồn tại. Khi nó tồn tại và là phân phối xác suất hợp lệ, là phân phối xác suất. Tuy nhiên, đây không phải là bản phân phối sau cho với khả năng dữ liệu đã cho . Các hậu thế cho trước đó không tồn tại. Chúng tôi chấp nhận trong phân tích của chúng tôi vì đó là một xấp xỉ. π ( X ) π ( X ) π ( θ | X ) π ( θ | X ) π ( θ ) π ( X | θ ) π ( θ | X )π(θX)π(X)π(X)π(θX)π(θX)π(θ)π(Xθ)π(θX)

1
@BenElizabethWard Tỷ lệ đã được sử dụng để chứng minh rằng trước đó vẫn chứa thông tin hữu ích mà chúng tôi có thể không thể tải vào một ưu tiên thích hợp. Tôi sẽ chỉnh sửa câu trả lời của tôi để bao gồm điều này.

2
@jsk không phải là phân phối xác suất, nhưng định nghĩa phân phối sau yêu cầu là phân phối xác suất, do đó, thật gian lận khi gọi là phân phối sau khi nó là một phân phối xác suất. Degroot & Schervish nói '.. chúng tôi sẽ tính toán các phân phối sau đó ..' mà họ cho rằng bạn đã đồng ý 'giả vờ rằng họ [các linh mục không đúng] là [các linh mục phù hợp]' như đã trình bày trước đó trong đoạn trích dẫn. π ( θ ) π ( θ X )π(θ)π(θ)π(θ|X)

1
Để làm cho câu trả lời của bạn hoàn chỉnh và khép kín để người đọc trong tương lai không phải đọc qua trao đổi nhận xét này, bạn có muốn cập nhật câu trả lời của mình không?
JSK

9

Có một câu trả lời "lý thuyết" và một câu trả lời "thực dụng".

Từ quan điểm trị liệu, khi một ưu tiên không phù hợp thì hậu thế không tồn tại (tốt, hãy xem câu trả lời của Matthew cho một tuyên bố âm thanh hơn), nhưng có thể được xấp xỉ bằng một hình thức giới hạn.

Nếu dữ liệu bao gồm một mẫu có điều kiện iid từ phân phối với tham số Bernoulli θ có sự phân bố beta với các thông số αβ , sự phân bố sau của θ là sự phân bố beta với các thông số α + s , β + n - s ( n quan sát, thành công) và trung bình của nó là ( α + s ) / ( α + β + n )θθαβθα+S,β+n-SnS(α+S)/(α+β+n). Nếu chúng ta sử dụng phân phối không đúng cách (và không có thật) beta trước với hypeparameters trước , và giả vờ rằng π ( θ ) α θ - 1 ( 1 - θ ) - 1 , chúng ta có được một sau đúng tỉ lệ với θ s - 1 ( 1 - θ ) n - s - 1 , tức là pdf của bản phân phối beta với các tham số sα= =β= =0π(θ)αθ-1(1-θ)-1θS-1(1-θ)n-S-1Sα 0 β 0n-Sngoại trừ một yếu tố không đổi. Đây là hình thức giới hạn của hậu thế cho bản beta trước với các tham số và (Degroot & Schervish, Ví dụ 7.3.13).α0β0

Trong một mô hình bình thường với giá trị trung bình , phương sai đã biết và phân phối trước cho , nếu độ chính xác trước, , nhỏ so với độ chính xác của dữ liệu, , khi đó phân phối sau xấp xỉ như thể : tức là phân phối sau là xấp xỉ kết quả từ giả định tỷ lệ với hằng số choσ 2 N ( μ 0 , τ 2 0 ) θ 1 / τ 2 0 n / σ 2 τ 2 0 = p ( θ x ) N ( θ ˉ x , σ 2 / n ) p ( θ ) θ ( - , ) τ 2 0θσ2N(μ0,τ02)θ1/τ02n/σ2τ02= =

p(θ|x)N(θ|x¯,σ2/n)
p(θ)θ(-,), một phân phối không hoàn toàn có thể, nhưng hình thức giới hạn của hậu thế như cách tiếp cận không tồn tại ( Gelman et al. , trang 52).τ02

Từ quan điểm "thực dụng" của xem, khi p ( x | θ ) = 0 bất cứ điều gì p ( θ ) là, vì vậy nếu p ( x | θ ) 0 trong ( một , b ) , sau đó - p ( x | θ ) p ( θp(x|θ)p(θ)= =0p(x|θ)= =0p(θ)p(x|θ)0(một,b) . Linh mục không đúng có thể được sử dụng để đại diện chohành viđịa phươngcủa phân phối trước trong khu vực nơi khả năng là đáng kể, nói ( a , b ) . Bằng cách giả sử rằng với một xấp xỉ đủ, một dạng trước theo sau như f ( x ) = k , x ( - , ) hoặc f-p(x|θ)p(θ)dθ= =mộtbp(x|θ)p(θ)dθ(một,b)f(x)= =k,x(-,) chỉ hơn ( a , b ) , rằng nó phù hợp với 0 bên ngoài phạm vi đó, chúng tôi đảm bảo rằng các linh mục thực sự được sử dụng là phù hợp (Box và Tiao, trang 21 ). Vì vậy, nếu phân phối trước của θ U ( - , ) nhưng ( một , b ) được bao bọc, nó là như nếu θ ~ U ( một ,f(x)= =kx-1,x(0,)(một,b)θBạn(-,)(một,b) , tức là p ( x | θ ) p ( θ ) = p ( x | θ ) k α p ( x | θ ) . Đối với một ví dụ cụ thể, đây là những gì xảy ra trongStan: nếu không có chỉ định trước cho một tham số, thì nó được mặc đồng phục trước sự hỗ trợ của nó và điều này được xử lý như một phép nhân của hằng số.θ~Bạn(một,b)p(x|θ)p(θ)= =p(x|θ)kαp(x|θ)


Bạn có thể nói thêm về lý do tại sao nó không tồn tại từ quan điểm lý thuyết?
JSK

Tôi không thể giải thích tốt hơn Matthew trong câu trả lời của anh ấy và trong những bình luận của anh ấy.
Sergio

Trong phần thực dụng, y là gì? Cũng trong phần đó, một số nên thuật ngữ được khả năng p ( x | q ) ? p(θ|x)p(x|θ)
JSK

Cảm ơn. Tôi nghĩ đó có thể là một sai lầm nữa ... Bạn viết , nhưng lon trước không phụ thuộc vào x . Bạn có nghĩa là P ( θ ) = k θ - 1 ? P(θ)= =kx-1xP(θ)= =kθ-1
JSK

Đúng! Tôi đã viết lại công thức đó khi chúng ở trong Box & Tiao. Tôi đã cố gắng chọn một ký hiệu đồng nhất (ví dụ Gelman sử dụng thay vì x , DeGroot sử dụng ξ ( . ) Cho priors và posteriors vv), nhưng tôi đã kết thúc trong một mớ hỗn độn ... Cảm ơn! yxξ(.)
Sergio

2

Tuy nhiên, trong trường hợp không đúng trước, làm thế nào để bạn biết rằng phân phối sau thực sự tồn tại?

Các hậu thế cũng có thể không thích hợp. Nếu ưu tiên là không phù hợp và khả năng là phẳng (vì không có quan sát có ý nghĩa), thì sau đó bằng với trước và cũng không đúng.

Thông thường bạn có một số quan sát, và thường thì khả năng là không bằng phẳng, do đó, hậu thế là phù hợp.

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.