compare | R Documentation |
This method is used primarily in the backend, and will only rarely be called by the end user. But see the examples below for a demonstration.
compare(numerator, denominator, data, ...)
numerator |
first model |
denominator |
second model (if omitted, compare to predefined null) |
data |
data for the comparison |
... |
arguments passed to and from related methods |
The compare function will return a model comparison object, typically a Bayes factor
## Sample from the posteriors for two models
data(puzzles)
## Main effects model; result is a BFmcmc object, inheriting
## mcmc from the coda package
mod1 = lmBF(RT ~ shape + color + ID, data = puzzles, whichRandom = "ID",
progress = FALSE, posterior = TRUE, iterations = 1000)
plot(mod1)
## Full model
mod2 = lmBF(RT ~ shape*color + ID, data = puzzles, whichRandom = "ID",
progress = FALSE, posterior = TRUE, iterations = 1000)
## Each BFmcmc object contains the model used to generate it, so we
## can compare them (data is not needed, it is contained in the objects):
compare(mod1, mod2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.