mudiff.acc.equalvar: Bayesian sample size determination for differences in normal...

mudiff.acc.equalvarR Documentation

Bayesian sample size determination for differences in normal means when variances are equal using the Average Coverage Criterion

Description

The function mudiff.acc.equalvar returns the required sample sizes to reach a given coverage probability on average for a posterior credible interval of fixed length for the difference between two normal means, when variances are equal.

Usage

mudiff.acc.equalvar(len, alpha, beta, n01, n02, level = 0.95, equal = TRUE)

Arguments

len

The desired fixed length of the posterior credible interval for the difference between the two unknown means

alpha

First prior parameter of the Gamma density for the common precision (reciprocal of the variance)

beta

Second prior parameter of the Gamma density for the common precision (reciprocal of the variance)

n01

Prior sample size equivalent for the mean for the first population

n02

Prior sample size equivalent for the mean for the second population

level

The desired average coverage probability of the posterior credible interval (e.g., 0.95)

equal

logical. Whether or not the final group sizes (n1, n2) are forced to be equal:

when equal = TRUE, final sample sizes n1 = n2;
when equal = FALSE, final sample sizes (n1, n2) minimize the expected posterior variance given a total of n1+n2 observations

Details

Assume that a sample from each of two populations will be collected in order to estimate the difference between two independent normal means. Assume that the precisions of the two normal sampling distributions are unknown but equal, with prior information in the form of a Gamma(alpha, beta) density. Assume that the means are unknown, but have prior information equivalent to (n01, n02) previous observations, respectively. The function mudiff.acc.equalvar returns the required sample sizes to attain the desired average coverage probability level for the posterior credible interval of fixed length len for the difference between the two unknown means.

This function uses a fully Bayesian approach to sample size determination. Therefore, the desired coverages and lengths are only realized if the prior distributions input to the function are used for final inferences. Researchers preferring to use the data only for final inferences are encouraged to use the Mixed Bayesian/Likelihood version of the function.

Value

The required sample sizes (n1, n2) for each group given the inputs to the function.

Note

The sample sizes returned by this function are exact.

Author(s)

Lawrence Joseph lawrence.joseph@mcgill.ca and Patrick Bélisle

References

Joseph L, Bélisle P.
Bayesian sample size determination for Normal means and differences between Normal means
The Statistician 1997;46(2):209-226.

See Also

mudiff.alc.equalvar, mudiff.modwoc.equalvar, mudiff.acc, mudiff.alc, mudiff.modwoc, mudiff.varknown, mudiff.mblacc.equalvar, mudiff.mblalc.equalvar, mudiff.mblmodwoc.equalvar, mudiff.mblacc, mudiff.mblalc, mudiff.mblmodwoc, mudiff.mbl.varknown, mudiff.freq, mu.acc, mu.alc, mu.modwoc, mu.varknown, mu.mblacc, mu.mblalc, mu.mblmodwoc, mu.mbl.varknown, mu.freq

Examples

mudiff.acc.equalvar(len=0.2, alpha=2, beta=2, n01=10, n02=50)

SampleSizeMeans documentation built on Aug. 23, 2023, 1:09 a.m.