implement.SBFNAP_onet {NAP} | R Documentation |
Implement Sequential Bayes Factor using the NAP for one-sample t
-tests
Description
In a N(\mu,\sigma^2)
population with unknown variance \sigma^2
, consider the two-sided one-sample t
-test for testing the point null hypothesis H_0 : \mu = 0
against H_1 : \mu \neq 0
. For a sequentially observed data, this function implements the Sequential Bayes Factor design when a normal moment prior is assumed on the standardized effect size \mu/\sigma
under the alternative.
Usage
implement.SBFNAP_onet(obs, tau.NAP = 0.3/sqrt(2),
RejectH1.threshold = exp(-3), RejectH0.threshold = exp(3),
batch.size, return.plot = TRUE, until.decision.reached = TRUE)
Arguments
obs |
Numeric vector. The vector of sequentially observed data. |
tau.NAP |
Positive numeric. Parameter in the moment prior. Default: |
RejectH1.threshold |
Positive numeric. |
RejectH0.threshold |
Positive numeric. |
batch.size |
Integer vector. The vector of batch sizes at each sequential comparison. The first element (the first batch size) needs to be at least 2. Default: |
return.plot |
Logical. Whether a sequential comparison plot to be returned. Default: |
until.decision.reached |
Logical. Whether the sequential comparison is performed until a decision is reached or until the data is observed. Default: |
Value
A list with three components named N
, BF
, and decision
.
$N
contains the number of sample size used.
$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.
Author(s)
Sandipan Pramanik and Valen E. Johnson
References
Pramanik, S. and Johnson, V. (2022). Efficient Alternatives for Bayesian Hypothesis Tests in Psychology. Psychological Methods. Just accepted.
Johnson, V. and Rossell, R. (2010). On the use of non-local prior densities in Bayesian hypothesis tests. Journal of the Royal Statistical Society: Series B, 72:143-170. [Article]
Examples
out = implement.SBFNAP_onet(obs = rnorm(100))