samplePowerPoisson: Compute the power of the Poisson rate test for a given sample...

In dkahle/bayesRates: Two-Sample Tests and Sample Size Determination from a Bayesian Perspective

Description

Compute the power of the Poisson rate test for a given sample size (t) with specified gamma priors

Usage

samplePowerPoisson(t, a1, b1, a2, b2, a = a1, b = b1, pi0 = 0.5,
  pi1 = 1 - pi0, c = 1)

Arguments

t	the sample size
a1	alpha, the hyperparameter of the gamma distribution of the first poisson rate
b1	beta, the hyperparameter of the gamma distribution of the first poisson rate
a2	alpha, the hyperparameter of the gamma distribution of the second poisson rate
b2	beta, the hyperparameter of the gamma distribution of the second poisson rate
a	alpha, the hyperparameter of the gamma distribution under the null
b	beta, the hyperparameter of the gamma distribution under the null
pi0	the prior probability of the null hypothesis
pi1	the prior probability of the alternative hypothesis
c	relative loss constant (loss due to type II error divided by loss due to type I error)

See Also

samplePower, samplePowerEst, findSize

Examples

## Not run: 

a1 <- 2; b1 <- 6
a2 <- 3; b2 <- 3
plot_gamma(c(a1, a2), c(b1, b2))
samplePowerPoisson(t = 10, a1, b1, a2, b2)
samplePowerPoisson(t = 11, a1, b1, a2, b2)
samplePowerPoisson(t = 10:11, a1, b1, a2, b2)
samplePowerPoisson(t = 250, a1, b1, a2, b2) # 56s
samplePowerEst(t = c(250, 500, 1000), a1, b1, a2, b2) # 5s, off by .001

a1 <- 4; b1 <- 2
a2 <- 6; b2 <- 2
plot_gamma(c(a1, a2), c(b1, b2))
samplePowerPoisson(t = 10, a1, b1, a2, b2)



## End(Not run)

dkahle/bayesRates documentation built on May 15, 2019, 9:07 a.m.