In case of two independent populations \(N(\mu_1,\sigma_0^2)\) and \(N(\mu_2,\sigma_0^2)\) with known common variance \(\sigma_0^2\), consider the two-sample \(z\)-test for testing the point null hypothesis of difference in their means \(H_0 : \mu_2 - \mu_1 = 0\) against \(H_1 : \mu_2 - \mu_1 \neq 0\). For a sequentially observed data, this function implements the Sequential Bayes Factor design when a normal moment prior is assumed on the difference between standardized effect sizes \((\mu_2 - \mu_1)/\sigma_0\) under the alternative.
implement.SBFHajnal_twoz(obs1, obs2, es1 = 0.3, sigma0 = 1,
RejectH1.threshold = exp(-3), RejectH0.threshold = exp(3),
batch1.size, batch2.size, return.plot = TRUE,
until.decision.reached = TRUE)
Numeric vector. The vector of sequentially observed data from Group-1.
Numeric vector. The vector of sequentially observed data from Group-2.
Positive numeric. \(\delta\) as above. Default: \(0.3\). For this, the prior on \((\mu_2 - \mu_1)/\sigma_0\) takes values \(0.3\) and \(-0.3\) each with equal probability 1/2.
Positive numeric. Known standard deviation in the population. Default: 1.
Positive numeric. \(H_0\) is accepted if \(BF \le\)RejectH1.threshold
. Default: exp(-3)
.
Positive numeric. \(H_0\) is rejected if \(BF \ge\)RejectH0.threshold
. Default: exp(3)
.
Integer vector. The vector of batch sizes from Group-1 at each sequential comparison. Default: rep(1, length(obs1))
.
Integer vector. The vector of batch sizes from Group-2 at each sequential comparison. Default: rep(1, length(obs2))
.
Logical. Whether a sequential comparison plot to be returned. Default: TRUE
.
Logical. Whether the sequential comparison is performed until a decision is reached or until the data is observed. Default: TRUE
. This means the comparison is performed until a decision is reached.
A list with three components named N1
, N2
, BF
, and decision
.
$N1
and $N2
contains the number of sample size used from Group-1 and 2.
$BF
contains the Bayes factor values at each sequential comparison.
$decision
contains the decision reached. 'A'
indicates acceptance of \(H_0\), 'R'
indicates rejection of \(H_0\), and 'I'
indicates inconclusive.
Hajnal, J. (1961). A two-sample sequential t-test.Biometrika, 48:65-75, [Article].
Schnuerch, M. and Erdfelder, E. (2020). A two-sample sequential t-test.Biometrika, 48:65-75, [Article].
# NOT RUN {
out = implement.SBFHajnal_twoz(obs1 = rnorm(100), obs2 = rnorm(100))
# }
Run the code above in your browser using DataLab