Sıfır hipotezi altında binom testleri simüle edilirken p değerlerinin eşit olmayan dağılımı


17

Sıfır hipotezi altında p-değeri dağılımının eşit olması gerektiğini duydum. Bununla birlikte, MATLAB'daki binom test simülasyonları, ortalama 0.5'ten (bu durumda 0.518) daha büyük olan homojen dağılımlardan çok farklı geri dönüşler sağlar: resim açıklamasını buraya girin

coin = [0 1];
success_vec = nan(20000,1);

for i = 1:20000
    success = 0;
    for j = 1:200
        success = success + coin(randperm(2,1));
    end
    success_vec(i) = success;
end

    p_vec = binocdf(success_vec,200,0.5);
    hist(p_vec);

Rasgele sayılar üretme şeklimi değiştirmeye çalışmak yardımcı olmadı. Burada herhangi bir açıklama gerçekten takdir ediyorum.


6
Dikkate alınması gereken bir nokta, binom testindeki p-değerlerinin sadece belirli ayrık değerleri alacağıdır (pay ayrık olduğu için): örnek olarak, deney başına sadece 20 deneme [coin çevirme] ile sadece 11 ayrı p- döndürülebilecek değerler. Bu olası p değeridir, bu nedenle deney başına n = 200 deneme, 101 ayrı p değeri vardır. n/2+1
James Stanley

Matlab'ın "binom testi" tam olarak ne yapar?
whuber

2
Görünüşe göre bu posterin binom testi, binocdfsadece binom uk.mathworks.com/help/stats/binocdf.html'nin CDF'si
konjugateprior

Yanıtlar:


20

pH0

James Stanley'nin yorumlarda belirttiği gibi, test istatistiğinin dağılımı ayrıktır, bu nedenle sonuç geçerli değildir. Kodunuzda hiç hata olmayabilir (bir histogramla ayrı bir dağıtım göstermese de, cdf veya pmf'yi veya daha iyisini her ikisini de göstermeye eğilimliydim).

F(x)=x

resim açıklamasını buraya girin

Bu dağılımı simüle etmek yerine tam olarak hesaplamak oldukça mümkündür - ancak ipucunu takip ettim ve bir simülasyon yaptım (senden daha büyük olsa da).

n

α


Glen ve @ JamesStanley teşekkürler! Tam olarak p-değeri dağılımının tekdüze olmadığı ve hipotez testi açısından sonuçları ne anlama geldiğini anlamaya çalışıyorum - ama bunun için wikipedia'ya
dalacağım

3
α

F(x)x

A.Donda, Glen_b - teşekkürler! Çok yardımcı oldun.
TanZor
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.