Ước tính khả năng tối đa cho phân phối nhị thức âm


11

Câu hỏi như sau:

Một mẫu ngẫu nhiên gồm n giá trị được thu thập từ phân phối nhị thức âm với tham số k = 3.

  1. Tìm ước lượng khả năng tối đa của tham số π.
  2. Tìm một công thức tiệm cận cho lỗi tiêu chuẩn của công cụ ước tính này.
  3. Giải thích tại sao phân phối nhị thức âm sẽ xấp xỉ bình thường nếu tham số k đủ lớn. Các thông số của xấp xỉ bình thường này là gì?

Công việc của tôi là như sau:
1. Tôi cảm thấy đây là điều mong muốn nhưng tôi không chắc liệu tôi có chính xác ở đây không hoặc liệu tôi có thể đưa thông tin này được cung cấp thêm không?

p(x)=(x1k1)πk(1π)xkL(π)=Πinp(xn|π)(π)=Σinln(p(xn|π))(π)=Σinkπ(xk)(1π)
  1. Tôi nghĩ rằng sau đây là những gì được yêu cầu. Đối với phần cuối cùng, tôi cảm thấy mình cần thay thế π^ bằng kx

    (π^)=kπ^2+x(1π^)2se(π^)=1(π^)se(π^)=π^2k(1π^)2x
  2. Tôi không thực sự chắc chắn làm thế nào để chứng minh điều này và vẫn đang nghiên cứu nó. Bất kỳ gợi ý hoặc liên kết hữu ích sẽ được đánh giá rất cao. Tôi cảm thấy như nó có liên quan đến thực tế là phân phối nhị thức âm có thể được xem như là một tập hợp các phân phối hình học hoặc nghịch đảo của phân phối nhị thức nhưng không biết cách tiếp cận nó.

Bất kỳ trợ giúp sẽ được đánh giá rất cao


(1) Để tìm ước tính khả năng tối đa bạn cần tìm nơi hàm khả năng đăng nhập đạt đến mức tối đa. Tính điểm số (đạo hàm đầu tiên của hàm khả năng đăng nhập liên quan đến ) là một sự khởi đầu - giá trị này sẽ lấy giá trị nào ở mức tối đa? (Và hãy nhớ rằng bạn không cần phải ước tính .)π^πk
Scortchi - Tái lập Monica

Tôi đã quên thêm vào đạo hàm của khả năng đăng nhập = 0 với mục đích tìm ra mức tối đa. Nếu tôi đã tìm ra điều này một cách chính xác (vẫn đang làm việc trên nó kể từ khi đăng), thì cái tôi có làkπΣi=0n(xik)(1π)=0
Syzorr

Cẩn thận:Cũng lưu ý rằng bắt đầu lúc 1.i=1nkπi=1n(xik)(1π)= ?i
Scortchi - Tái lập Monica

Trong (2), hiếm khi xảy ra trường hợp đối ứng của sự khác biệt là sự khác biệt của các đối ứng. Lỗi này ảnh hưởng lớn đến công thức cuối cùng của bạn cho . se(π^)
whuber

Câu trả lời:


6

1.

p(x)=(xi1k1)πk(1π)xik

L(π;xi)=i=1n(xi1k1)πk(1π)xik

(π;xi)=i=1n[log(xi1k1)+klog(π)+(xik)log(1π)]d(π;xi)dπ=i=1n[kπ(xik)(1π)]

Đặt cái này thành 0,

nkπ=i=1nxink1π

π^=nki=1nx

    2.

Đối với phần thứ hai, bạn cần sử dụng định lý , là thông tin câu cá ở đây. Do đó, độ lệch chuẩn của sẽ là . Hoặc bạn gọi đó là lỗi tiêu chuẩn kể từ khi bạn sử dụng CLT tại đây.n(θ^θ)DN(0,1I(θ))I(θ)θ^[nI(θ)]1/2

Vì vậy, chúng ta cần tính toán thông tin Fisher cho phân phối nhị thức âm.

2log(P(x;π))π2=kπ2xk(1π)2

I(θ)=E(kπ2xk(1π)2)=kπ2+k(1π)(1π)2π

Lưu ý: cho nhị phân nhị phân âmE(x)=kπ

Do đó, lỗi tiêu chuẩn cho làπ^[n(kπ2+k(1π)(1π)2π)]1/2

Đơn giản hóa chúng ta có được chúng ta sẽ cóse(π)=π2(π1)kn

    3.

Phân phối hình học là trường hợp đặc biệt của phân phối nhị thức âm khi k = 1. Lưu ý là phân phối hình họcπ(1π)x1

Do đó, biến nhị thức âm có thể được viết dưới dạng tổng của k biến ngẫu nhiên, phân phối giống hệt (hình học).

Vì vậy, phân phối nhị thức âm CLT sẽ xấp xỉ bình thường nếu tham số k đủ lớn


1
Xin vui lòng đọc những chủ đề tôi có thể hỏi về đây? về các câu hỏi tự học: thay vì làm bài tập về nhà cho mọi người, chúng tôi cố gắng giúp họ tự làm.
Scortchi - Tái lập Monica

2
Bạn làm cần phải xem xét kích thước mẫu khi tính MLE. Bạn có thể nhầm lẫn một tài khoản của quan sát độc lập, mỗi quan sát không. các thử nghiệm cần thiết để đạt được thất bại ( ) với một tài khoản về một quan sát duy nhất về không. các thử nghiệm cần thiết để đạt được thất bại ( ). Cái trước cho khả năng ; cái sau, . nnkx1,x2,,xnkni=1nπ(1π)xikπk(1π)nk
Scortchi - Phục hồi Monica

1
Bạn nói đúng, tôi luôn bối rối về phần này. Cảm ơn rât nhiều. Tôi cũng hỏi rất nhiều câu hỏi trên diễn đàn này, nhưng tôi thực sự hy vọng mọi người có thể cho tôi câu trả lời rất chi tiết, sau đó tôi có thể tự mình nghiên cứu nó.
Sâu Bắc

Vâng. Tôi hiểu tại sao quy tắc chống lại việc cung cấp quá nhiều chi tiết nhưng câu trả lời này kết hợp với ghi chú của riêng tôi từ bài giảng đã cho phép tôi buộc rất nhiều kết thúc lỏng lẻo lại với nhau. Tôi dự định sẽ đi và nói chuyện với giảng viên của tôi ngày hôm nay về điều này để tôi có thể làm rõ từ anh ấy. Bây giờ là thứ sáu ở đây. Chuyển nhượng do thứ hai như đã nêu ở trên. Chúng tôi đã học được điều này vào thứ Tư và chỉ có một ví dụ duy nhất sử dụng phân phối nhị thức. Cảm ơn rất nhiều cho các chi tiết.
Syzorr

Có một số lỗi trong công việc của bạn ở đó vì tôi () = E [] không -E [] (điều này làm tôi bối rối cho đến khi tôi tìm kiếm các phương trình bạn đã sử dụng) Cuối cùng đã kết thúc vớise(π)=π2(π1)kn
Syzorr
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.