1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/nWayAOV.R
\name{nWayAOV}
\alias{nWayAOV}
\title{Use ANOVA design matrix to compute Bayes factors or sample posterior}
\usage{
nWayAOV(
y,
X,
gMap,
rscale,
iterations = 10000,
progress = getOption("BFprogress", interactive()),
callback = function(...) as.integer(0),
gibbs = NULL,
posterior = FALSE,
ignoreCols = NULL,
thin = 1,
method = "auto",
continuous = FALSE,
noSample = FALSE
)
}
\arguments{
\item{y}{vector of observations}
\item{X}{design matrix whose number of rows match \code{length(y)}.}
\item{gMap}{vector grouping the columns of \code{X} (see Details).}
\item{rscale}{a vector of prior scale(s) of appropriate length (see
Details).}
\item{iterations}{Number of Monte Carlo samples used to estimate Bayes
factor or posterior}
\item{progress}{if \code{TRUE}, show progress with a text progress bar}
\item{callback}{callback function for third-party interfaces}
\item{gibbs}{will be deprecated. See \code{posterior}}
\item{posterior}{if \code{TRUE}, return samples from the posterior using
Gibbs sampling, instead of the Bayes factor}
\item{ignoreCols}{if \code{NULL} and \code{posterior=TRUE}, all parameter
estimates are returned in the MCMC object. If not \code{NULL}, a vector of
length P-1 (where P is number of columns in the design matrix) giving which
effect estimates to ignore in output}
\item{thin}{MCMC chain to every \code{thin} iterations. Default of 1 means
no thinning. Only used if \code{posterior=TRUE}}
\item{method}{the integration method (only valid if \code{posterior=FALSE}); one
of "simple", "importance", "laplace", or "auto"}
\item{continuous}{either FALSE if no continuous covariates are included, or
a logical vector of length equal to number of columns of X indicating
which columns of the design matrix represent continuous covariates}
\item{noSample}{if \code{TRUE}, do not sample, instead returning NA. This is
intended to be used with functions generating and testing many models at one time,
such as \code{\link{anovaBF}}}
}
\value{
If \code{posterior} is \code{FALSE}, a vector of length 2 containing
the computed log(e) Bayes factor (against the intercept-only null), along
with a proportional error estimate on the Bayes factor. Otherwise, an
object of class \code{mcmc}, containing MCMC samples from the posterior is
returned.
}
\description{
Computes a single Bayes factor, or samples from the posterior, for an ANOVA
model defined by a design matrix
}
\details{
This function is not meant to be called by end-users, although
technically-minded users can call this function for flexibility beyond what
the other functions in this package provide. See \code{\link{lmBF}} for a
user-friendly front-end to this function. Details about the priors can be
found in the help for \code{\link{anovaBF}} and the references therein.
Argument \code{gMap} provides a way of grouping columns of
the design matrix as a factor; the effects in each group will share a common
\eqn{g} parameter. \code{gMap} should be a vector of the same length as the number of
nonconstant rows in \code{X}. It will contain all integers from 0 to
\eqn{N_g-1}{Ng-1}, where \eqn{N_g}{Ng} is the total number of \eqn{g}
parameters. Each element of \code{gMap} specifies the group to which that
column belongs.
If all columns belonging to a group are adjacent, \code{struc} can instead
be used to compactly represent the groupings. \code{struc} is a vector of
length \eqn{N_g}{Ng}. Each element specifies the number columns in the
group.
The vector \code{rscale} should be of length \eqn{N_g}{Ng}, and contain the
prior scales of the standardized effects. See Rouder et al. (2012) for more
details and the help for \code{\link{anovaBF}} for some typical values.
The method used to estimate the Bayes factor depends on the \code{method}
argument. "simple" is most accurate for small to moderate sample sizes, and
uses the Monte Carlo sampling method described in Rouder et al. (2012).
"importance" uses an importance sampling algorithm with an importance
distribution that is multivariate normal on log(g). "laplace" does not
sample, but uses a Laplace approximation to the integral. It is expected to
be more accurate for large sample sizes, where MC sampling is slow. If
\code{method="auto"}, then an initial run with both samplers is done, and
the sampling method that yields the least-variable samples is chosen. The
number of initial test iterations is determined by
\code{options(BFpretestIterations)}.
If posterior samples are requested, the posterior is sampled with a Gibbs
sampler.
}
\note{
Argument \code{struc} has been deprecated. Use \code{gMap}, which is the \code{\link{inverse.rle}} of \code{struc},
minus 1.
}
\examples{
## Classical example, taken from t.test() example
## Student's sleep data
data(sleep)
plot(extra ~ group, data = sleep)
## traditional ANOVA gives a p value of 0.00283
summary(aov(extra ~ group + Error(ID/group), data = sleep))
## Build design matrix
group.column <- rep(1/c(-sqrt(2),sqrt(2)),each=10)
subject.matrix <- model.matrix(~sleep$ID - 1,data=sleep$ID)
## Note that we include no constant column
X <- cbind(group.column, subject.matrix)
## (log) Bayes factor of full model against grand-mean only model
bf.full <- nWayAOV(y = sleep$extra, X = X, gMap = c(0,rep(1,10)), rscale=c(.5,1))
exp(bf.full[['bf']])
## Compare with lmBF result (should be about the same, give or take 1\%)
bf.full2 <- lmBF(extra ~ group + ID, data = sleep, whichRandom = "ID")
bf.full2
}
\references{
Rouder, J. N., Morey, R. D., Speckman, P. L., Province, J. M.,
(2012) Default Bayes Factors for ANOVA Designs. Journal of Mathematical
Psychology. 56. p. 356-374.
}
\seealso{
See \code{\link{lmBF}} for the user-friendly front end to this
function; see \code{\link{regressionBF}} and \code{anovaBF} for testing
many regression or ANOVA models simultaneously.
}
\author{
Richard D. Morey (\email{richarddmorey@gmail.com}), Jeffery N.
Rouder (\email{rouderj@missouri.edu})
}
\keyword{htest}
|
