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