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
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/contingency.R
\name{contingencyTableBF}
\alias{contingencyTableBF}
\title{Function for Bayesian analysis of one- and two-sample designs}
\usage{
contingencyTableBF(
x,
sampleType,
fixedMargin = NULL,
priorConcentration = 1,
posterior = FALSE,
callback = function(...) as.integer(0),
...
)
}
\arguments{
\item{x}{an m by n matrix of counts (integers m,n > 1)}
\item{sampleType}{the sampling plan (see details)}
\item{fixedMargin}{for the independent multinomial sampling plan, which margin is fixed ("rows" or "cols")}
\item{priorConcentration}{prior concentration parameter, set to 1 by default (see details)}
\item{posterior}{if \code{TRUE}, return samples from the posterior instead
of Bayes factor}
\item{callback}{callback function for third-party interfaces}
\item{...}{further arguments to be passed to or from methods.}
}
\value{
If \code{posterior} is \code{FALSE}, an object of class
\code{BFBayesFactor} containing the computed model comparisons is
returned.
If \code{posterior} is \code{TRUE}, an object of class \code{BFmcmc},
containing MCMC samples from the posterior is returned.
}
\description{
This function computes Bayes factors for contingency tables.
}
\details{
The Bayes factor provided by \code{contingencyTableBF} tests the independence assumption in
contingency tables under various sampling plans, each of which is described below.
See Gunel and Dickey (1974) for more details.
For \code{sampleType="poisson"}, the sampling plan is assumed to be
one in which observations occur as a poisson process with an overall
rate, and then assignment to particular factor levels occurs with
fixed probability. Under the null hypothesis, the assignments to the
two factors are independent. Importantly, the total N is not fixed.
For \code{sampleType="jointMulti"} (joint multinomial), the sampling
plan is assumed to be one in which the total N is fixed, and observations
are assigned to cells with fixed probability. Under the null hypothesis, the
assignments to the two factors are independent.
For \code{sampleType="indepMulti"} (independent multinomial), the
sampling plan is assumed to be one in which row or column totals are fixed,
and each row or column is assumed to be multinomially distributed.
Under the null hypothesis, each row or column is assumed to have the
same multinomial probabilities. The fixed margin must be given by
the \code{fixedMargin} argument.
For \code{sampleType="hypergeom"} (hypergeometric), the sampling
plan is assumed to be one in which both the row and column totals are fixed.
Under the null hypothesis, the cell counts are assumed to be governed by the
hypergeometric distribution.
For all models, the argument \code{priorConcentration} indexes
the expected deviation from the null hypothesis under the alternative,
and corresponds to Gunel and Dickey's (1974) "a" parameter.
}
\note{
Posterior sampling for the hypergeometric model under the alternative
has not yet been implemented.
}
\examples{
## Hraba and Grant (1970) doll race data
data(raceDolls)
## Compute Bayes factor for independent binomial design, with
## columns as the fixed margin
bf = contingencyTableBF(raceDolls, sampleType = "indepMulti", fixedMargin = "cols")
bf
## Posterior distribution of difference in probabilities, under alternative
chains = posterior(bf, iterations = 10000)
sameRaceGivenWhite = chains[,"pi[1,1]"] / chains[,"pi[*,1]"]
sameRaceGivenBlack = chains[,"pi[1,2]"] / chains[,"pi[*,2]"]
hist(sameRaceGivenWhite - sameRaceGivenBlack, xlab = "Probability increase",
main = "Increase in probability of child picking\nsame race doll (white - black)",
freq=FALSE, yaxt='n')
box()
}
\references{
Gunel, E. and Dickey, J., (1974)
Bayes Factors for Independence in Contingency Tables. Biometrika, 61, 545-557
}
\author{
Richard D. Morey (\email{richarddmorey@gmail.com})
Tahira Jamil (\email{tahjamil@gmail.com})
}
\keyword{htest}
|
