Nilai yang diharapkan dari estimasi parameter koin kemungkinan maksimum


8

Misalkan saya memiliki percobaan lemparan koin di mana saya ingin menghitung perkiraan kemungkinan maksimum dari parameter koin saat melempar koin kali. Setelah menghitung turunan dari fungsi kemungkinan binomial L (p) = {n \ pilih x} p ^ x (1-p) ^ {nx} , saya mendapatkan nilai optimal untuk p menjadi p ^ {*} = \ frac {x} {n} , dengan x menjadi jumlah keberhasilan.pnL(p)=(nx)px(1p)nxpp=xnx

Pertanyaan saya sekarang adalah:

  • Bagaimana saya menghitung nilai / varian yang diharapkan dari estimasi kemungkinan maksimum ini untuk p ?
  • Apakah saya perlu menghitung nilai / varian yang diharapkan untuk L(p) ?
  • Jika ya, bagaimana saya melakukannya?

1
Saya kira ini semacam belajar mandiri (Anda harus menandainya seperti itu). Apa sebenarnya yang kamu inginkan? Apakah inferensi pada parameter Anda?
pkofod

Apa yang Anda maksud inferensi pada parameter? Saya tidak yakin bagaimana saya akan menghitung nilai / varian yang diharapkan untuk kuantitas . Maksud saya, saya tahu apa arti / varians dan bagaimana menghitungnya untuk contoh sederhana, tetapi tidak mengerti bagaimana menerapkannya ke . pp
Manu

Jawaban:


6

Pertama-tama, ini adalah pertanyaan belajar mandiri, jadi saya akan membahas terlalu banyak detail teknis, tetapi saya juga tidak akan sibuk dengan derivasi. Ada banyak cara untuk melakukan ini. Saya akan membantu Anda dengan menggunakan properti umum dari estimator kemungkinan maksimum.

Informasi latar belakang

Untuk menyelesaikan masalah Anda, saya pikir Anda perlu mempelajari kemungkinan maksimum sejak awal. Anda mungkin menggunakan beberapa jenis buku teks, dan jawabannya harus benar-benar ada di suatu tempat. Saya akan membantu Anda mencari tahu apa yang harus dicari.

Kemungkinan Maksimum adalah metode estimasi yang pada dasarnya adalah apa yang kita sebut sebagai estimator-M (anggap "M" sebagai "maksimalkan / perkecil"). Jika kondisi yang diperlukan untuk menggunakan metode ini dipenuhi, kami dapat menunjukkan bahwa estimasi parameter konsisten dan terdistribusi secara asimptotik, sehingga kami memiliki:

N(θ^θ0)dNormal(0,A01B0A01),

di mana dan adalah beberapa matriks. Ketika menggunakan kemungkinan maksimum, kita dapat menunjukkan bahwa , dan dengan demikian kita memiliki ungkapan sederhana: Kami memiliki mana menunjukkan hessian. Inilah yang perlu Anda perkirakan untuk mendapatkan varians Anda.A0B0A0=B0

N(θ^θ0)dNormal(0,A01).
A0E(H(θ0))H

Masalah spesifik Anda

jadi bagaimana kita melakukannya? Di sini mari kita sebut vektor parameter kami apa yang Anda lakukan: . Ini hanya skalar, jadi "skor" kami hanyalah turunan dan "hessian" hanyalah turunan urutan kedua. Fungsi kemungkinan kami dapat ditulis sebagai: yang ingin kami maksimalkan. Anda menggunakan turunan pertama dari ini atau kemungkinan log untuk menemukan . Alih-alih menetapkan turunan pertama sama dengan nol, kita dapat membedakan lagi, untuk menemukan turunan urutan kedua . Pertama kita mengambil log: Kemudian 'skor' kami adalah: dan 'hessian' kami: θp

l(p)=(p)x(1p)nx,
pH(p)
ll(p)log(l(p))=xlog(p)+(nx)log(1p)
ll(p)=xp+nx1p,
H(p)=ll(p)=xp2nx(1p)2.
Maka teori umum kami dari atas hanya memberitahu Anda untuk menemukan . Sekarang Anda hanya perlu mengambil ekspektasi (Petunjuk: gunakan ), kalikan dengan dan ambil kebalikannya. Maka Anda akan memiliki varian penduga Anda.(E(H(p)))1H(p)E(x/n)=p1

Apakah benar? Var(p)=p21n1np
Manu

@ Manu: Tidak cukup, tapi sepertinya Anda baru saja membuat kesalahan di suatu tempat. Bisakah Anda memposting beberapa langkah lagi?
pkofod

(E(H(p)))1=E(H(p))1]=(E(xp2)+E(nx(1p)2))1=(p2np+(1p)2(nnp))1 . dari sana saya disederhanakan dengan mengalikan dan mengambil kebalikannya.
Manu

Itu semua benar, sekarang hanya menyederhanakan. Di bagian pertama p batal, dan di bagian kedua Anda dapat mengambil n di luar tanda kurung.
pkofod

2
(n/p+n/[1p])1 adalah apa yang Anda miliki di atas. Hanya faktor keluar, pakai denominator umum dan kemudian ambil kebalikannya. n
ekvall

6

Untuk memulai Anda, mari kita lakukan nilai yang diharapkan:

Jika adalah jumlah keberhasilan dalam melempar, maka adalah proporsi keberhasilan dalam sampel Anda. Pertimbangkan ; untuk setiap lemparan, probabilitas keberhasilan adalah sesuai dengan asumsi, jadi ketika melempar koin satu kali "jumlah keberhasilan" yang diharapkan adalah , kan? Jadi, jika Anda melempar koin kali, Anda akan mengharapkan kesuksesan kali karena lemparannya independen. Kemudian, karena adalah jumlah keberhasilan yang diharapkan dalam lemparan, Anda mendapatkanxnx/nExpp×1+(1p)×0=pnnpnpn

Ep=En1x=n1Ex=n1×np=p

Jadi estimatornya tidak bias. Bisakah Anda mengetahui cara melakukan varians dari sini?

Sunting: Mari kita lakukan varians juga. Kami menggunakan . Istilah kedua yang kita miliki dari perhitungan nilai yang diharapkan, jadi mari kita lakukan yang pertama: Untuk menyederhanakan beberapa , kita dapat menyatakan jumlah keberhasilan dalam melempar sebagai berikut: mana mengambil nilai 1 jika melempar sukses dan 0 sebaliknya. Karenanya, dan dengan demikian menyatukan semuanya Anda tiba di .Var(p)=Ep2(Ep)2

Ep2=n2Ex2
n
x=1nχi,
χii
Ex2=E(1nχi)2=E[1nχi2+2i<jχiχj]=np+n(n1)p2,
Var(p)=p(1p)n

Jika Anda melempar kepala berturut-turut, . Namun berapa nilai tepatnya yang akan diambil Var ( )? n=3pMLE=1.0p
piccolo
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.