Because of the discreteness issue in one-sample proportion tests, it is not always 'effective' (the sense is explained in the Details section) to use just any value as a maximum sample size for designing a MSPRT. Suppose, we have a desired value for the maximum sample size. Given this, the function finds a value, defined as the effective maximum sample size, which should be used as the maximum sample size for designing the MSPRT (as N.max
in design.MSPRT()
).
effective.N(N, side = "right", type1 = 0.005, null = 0.5, plot.it = T)
a numeric; a desired value of the maximum sample size for designing the MSPRT.
a character; direction of the alternative hypothesis H1.
Has to be one of "right"
or "left"
.
Default: "right"
.
a numeric in (0,1)
; the probability at which we want to control the Type 1 error of the MSPRT.
Default: 0.005
.
a numeric; the hypothesized value of proportion under the null hypothesis.
Default: 0.5.
a logical; if TRUE, returns a plot; otherwise it doesn't.
Default: TRUE.
Returns a numeric, the effective maximum sample size.
Suppose we are provided with a simple null hypothesis, and Type 1 & 2 error probabilities. Because of the discreteness issue in one-sample proportion tests in a fixed design, the fixed design alternative will not always decrease even if we are increasing the sample size. So, we first shortlist only those values from 1 to N
which results in strictly decreasing UMPBT point alternatives (as is originally defined in Johnson (2013)). The effective maximum sample size is then chosen to be the maximum among those shorlisted value. So, obviously, the effective maximum sample size is either N
or smaller than that.
Main article and supplemental file of MSPRT
Johnson, Valen E., Uniformly most powerful Bayesian tests., Ann. of Stat., 41, (4), 2013, pp. 1716-1741
# NOT RUN {
effective.N(N = 30, null = .2)
# }
Run the code above in your browser using DataLab