Nothing
## ----fig.height=3.5, fig.width=5----------------------------------------------
# load package
library("metaBMA")
# load data set
data(towels)
# Half-normal (truncated to > 0)
p1 <- prior("norm", c(mean=0, sd=.3), lower = 0)
p1
p1(1:3)
plot(p1)
# custom prior
p1 <- prior("custom", function(x) x^3-2*x+3, lower = 0, upper = 1)
plot(p1, -.5, 1.5)
## ----fig.height=3.5, fig.width=5----------------------------------------------
# Fixed-effects
progres <- capture.output( # suppress Stan progress for vignette
mf <- meta_fixed(logOR, SE, study, towels,
d = prior("norm", c(mean=0, sd=.3), lower=0))
)
mf
# plot posterior distribution
plot_posterior(mf)
## ----fig.height=3.5, fig.width=5----------------------------------------------
# Random-effects
progres <- capture.output( # suppress Stan progress for vignette
mr <- meta_random(logOR, SE, study, towels,
d = prior("norm", c(mean=0, sd=.3), lower=0),
tau = prior("t", c(location=0, scale=.3, nu=1), lower=0),
iter = 1500, logml_iter = 2000, rel.tol = .1)
)
mr
# plot posterior distribution
plot_posterior(mr, main = "Average effect size d")
plot_posterior(mr, "tau", main = "Heterogeneity tau")
## ----fig.height=4.5, fig.width=6----------------------------------------------
mb <- meta_bma(logOR, SE, study, towels,
d = prior("norm", c(mean=0, sd=.3), lower=0),
tau = prior("t", c(location=0, scale=.3, nu=1), lower=0),
iter = 1500, logml_iter = 2000, rel.tol = .1)
mb
plot_posterior(mb, "d", -.1, 1.4)
plot_forest(mb)
## ----eval = FALSE, fig.height=4.5, fig.width=6--------------------------------
# mp <- predicted_bf(mb, SE = .2, sample = 30)
# plot(mp)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.