Saya memiliki dua variabel acak poisson independen, dan , dengan dan . Saya ingin menguji versus alternatif .
Saya sudah mendapatkan estimasi kemungkinan maksimum berdasarkan hipotesis nol dan alternatif (model), dan berdasarkan pada itu saya menghitung statistik uji rasio kemungkinan (LRT) (kode R diberikan di bawah).
Sekarang saya tertarik untuk menghitung kekuatan tes berdasarkan:
- Alfa tetap (kesalahan tipe 1) = 0,05.
- Dengan menggunakan ukuran sampel yang berbeda (n), katakanlah n = 5, 10, 20, 50, 100.
- Kombinasi berbeda dari dan , yang akan mengubah statistik LRT (dihitung sebagai berikut).
LRTstat
Ini kode R saya:
X1 = rpois(λ1); X2 = rpois(λ2)
Xbar = (X1+X2)/2
LLRNum = dpois(X1, X1) * dpois(X2, X2)
LLRDenom = dpois(X1, Xbar) * dpois(X2, Xbar)
LRTstat = 2*log(LLRNum/LLRDenom)
Dari sini, bagaimana saya bisa melanjutkan dengan perhitungan daya (lebih disukai di R)?