Trợ giúp Tối đa hóa kỳ vọng từ giấy: làm thế nào để bao gồm phân phối trước?


9

Câu hỏi dựa trên bài báo có tiêu đề: Tái tạo hình ảnh trong chụp cắt lớp quang học khuếch tán bằng mô hình khuếch tán vận chuyển bức xạ kết hợp

Liên kết tải xuống

Các tác giả áp dụng thuật toán EM với độ chính thưa của của một vectơ không xác định để ước tính các pixel của hình ảnh. Mô hình được đưa ra bởitôi1μ

(1)y= =Mộtμ+e
Ước tính được đưa ra trong biểu thức (8) là

(2)μ^= =tranh luậnmmộtxlnp(y|μ)+γlnp(μ)

Trong trường hợp của tôi, tôi đã coi là một bộ lọc có độ dài và là các vectơ đại diện cho các bộ lọc. Vì thế,L μ L × 1μLμL×1

Mô hình có thể được viết lại thành

(3)y(n)= =μTmột(n)+v(n)

Câu hỏi: Xây dựng vấn đề: (n by 1) là đầu vào không quan sát được và là giá trị trung bình bằng 0 với phương sai không xác định nhiễu phụ gia. Giải pháp MLE sẽ dựa trên Tối đa hóa kỳ vọng (EM). { e ( n ) } σ 2 eμ(n){e(n)}σe2

Trong bài báo Eq (19) là hàm - khả năng ghi nhật ký hoàn chỉnh nhưng đối với trường hợp của tôi, tôi không hiểu làm thế nào tôi có thể đưa phân phối của vào biểu thức khả năng ghi nhật ký hoàn chỉnh. Một , μMộtMột,μ

Điều gì sẽ là khả năng đăng nhập hoàn chỉnh sử dụng EM của bao gồm phân phối trước?y


Bạn có thực sự muốn khả năng đăng nhập hay thay vào đó bạn muốn đăng nhập sau? Chỉ sau này sẽ bao gồm Laplacian trước. Cái trước chỉ có thể có được bằng cách lấy nhật ký của khả năng, mà dường như bạn đã viết ra

Có hai biểu thức mà tôi muốn - (1) Một biểu thức sẽ được sử dụng để tìm Ma trận thông tin Fisher và (2) khác sẽ là pdf của tập dữ liệu hoàn chỉnh bao gồm biến ẩn và các quan sát là khớp mật độ xác suất của dữ liệu quan sát là một hàm của tham số . Bản pdf mà tôi đã viết có thể áp dụng cho mô hình MA để ước tính mù . Nhưng, nó sẽ khác nhau như thế nào đối với ràng buộc về độ thưa thớt = Laplacian trước để Ma trận thông tin Fisher từ các dẫn xuất một phần của khả năng đăng nhập có thể được tìm thấy. θ θZθθ
SKM

@ Tây An: Tôi không hiểu làm thế nào để cắm 3 pdf, trong đó bao gồm trước khi xây dựng khả năng đăng nhập. Tôi có thể tìm ra mức tối đa hóa để lấy đạo hàm riêng và bằng 0. Bạn có thể vui lòng đưa ra một câu trả lời với biểu thức khả năng được viết rõ ràng. Điều này thực sự có ích
SKM

Câu trả lời:


3

log L ( θ | x ) = E [ log L ( θ | x , Z ) | x , θ ] - E [ log q

tranh luậntối đaθL(θ|x)π(θ)= =tranh luậntối đaθđăng nhậpL(θ|x)+đăng nhậpπ(θ)
đăng nhậpL(θ|x)= =E[đăng nhậpL(θ|x,Z)|x,θ]-E[đăng nhậpq(Z|x,θ)|x,θ]
θ
q(z|x,θ)= =f(x,z|θ)/g(x|θ)
g(x|θ)= =f(x,z|θ)/q(z|x,θ)
zZ
đăng nhậpg(x|θ)= =đăng nhậpf(x,z|θ)-đăng nhậpq(z|x,θ)= =E[đăng nhậpf(x,Z|θ)-đăng nhậpq(Z|x,θ)|x]
Zđã cho , ví dụ . Do đó, nếu chúng tôi tối đa hóa trong với giải pháp chúng tôi có trong khi bằng các đối số chuẩn của EM. Do đó, và sử dụng làm E bước vào mục tiêu X= =xq(z|x,θ)θ
E[đăng nhậpL(θ|x,Z)|x,θ]+đăng nhậpπ(θ)
θ1
E[đăng nhậpL(θ1|x,Z)|x,θ]+đăng nhậpπ(θ1)E[đăng nhậpL(θ|x,Z)|x,θ]+đăng nhậpπ(θ)
E[đăng nhậpq(Z|x,θ)|x,θ]E[đăng nhậpq(Z|x,θ1)|x,θ]
E[đăng nhậpL(θ1|x,Z)|x,θ]+đăng nhậpπ(θ1)E[đăng nhậpL(θ|x,Z)|x,θ]+đăng nhậpπ(θ)
E[đăng nhậpL(θ|x,Z)|x,θ]+đăng nhậpπ(θ)
dẫn đến sự gia tăng về phía sau ở mỗi bước M, có nghĩa là thuật toán EM đã sửa đổi hội tụ đến một MAP địa phương.

Cảm ơn bạn đã trả lời của bạn. Có đại diện cho pdf của không? Bạn có thể vui lòng tại sao có 2 kỳ vọng với Bị trừ trong phương trình được đề cập trong dòng thứ hai không? Z E [ l o g q ( . ) ]q()ZE[tôiogq(.)]
SKM

Tôi đã thêm một số giải thích, nhưng bạn nên kiểm tra trong sách giáo khoa về đạo hàm của thuật toán EM vì đây là tài liệu chuẩn.
Tây An

1

Tôi không nghĩ rằng việc hiển thị log-post tăng đơn điệu (hoặc khả năng đăng nhập cho MLE) là đủ để hiển thị sự hội tụ đến điểm dừng của ước tính MAP (hoặc MLE). Ví dụ, số gia có thể trở nên nhỏ tùy ý. Trong bài báo nổi tiếng của Wu 1983 , một điều kiện đủ để hội tụ đến điểm dừng của EM là sự khác biệt trong cả hai đối số của hàm ràng buộc dưới.

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.