mardi 3 novembre 2020

Compare the power of the Bonferoni and needle in the haystack

I have the sample of 100 from Poisson distribution. I have 1000 of similar hypothesis and the meta problem of testing the global hypothesis

⋂H0i: E(Xi)=5 vs. ⋂H1i: E(Xi)>5

And I want to use these simulations to compare the power of the Bonferoni and Fisher tests for: Needle in the haystack E(X1) = 7, E(Xj ) = 5 for j ∈ {2, . . . , 1000}

I tried calculating the power for Bonferroni and the pv values like this, but I have no idea how to compare it to the Needle in the haystack. Thank you for your help!

alfa=0.05
sample <- rpois(100,5)
pv <- 1-ppois(5,5)
power <- 1 - exp(-alfa)

Aucun commentaire:

Enregistrer un commentaire