ref.synth {MBNMAdose} | R Documentation |
Synthesise single arm dose = 0 / placebo studies to estimate E0
Description
Synthesises single arm studies to estimate E0. Used in predicting responses from a dose-response MBNMA.
Usage
ref.synth(
data.ab,
mbnma,
synth = "fixed",
n.iter = mbnma$BUGSoutput$n.iter,
n.burnin = mbnma$BUGSoutput$n.burnin,
n.thin = mbnma$BUGSoutput$n.thin,
n.chains = mbnma$BUGSoutput$n.chains,
...
)
Arguments
data.ab |
A data frame of arm-level data in "long" format containing the columns:
|
mbnma |
An S3 object of class |
synth |
A character object that can take the value |
n.iter |
number of total iterations per chain (including burn in; default: 2000) |
n.burnin |
length of burn in, i.e. number of iterations to
discard at the beginning. Default is |
n.thin |
thinning rate. Must be a positive integer. Set
|
n.chains |
number of Markov chains (default: 3) |
... |
Arguments to be sent to |
Details
data.ab
can be a collection of studies that closely resemble the
population of interest intended for the prediction, which could be
different to those used to estimate the MBNMA model, and could include
single arms of RCTs or observational studies. If other data is not
available, the data used to estimate the MBNMA model can be used by
selecting only the studies and arms that investigate dose = 0 (placebo).
Defaults for n.iter
, n.burnin
, n.thin
and n.chains
are those used to estimate
mbnma
.
Value
A list of named elements corresponding to E0 and the between-study standard deviation for
E0 if synth="random"
. Each element contains the full MCMC results from the synthesis.
Examples
# Using the triptans data
network <- mbnma.network(triptans)
# Run an Emax dose-response MBNMA
emax <- mbnma.run(network, fun=demax(), method="random")
# Data frame for synthesis can be taken from placebo arms
ref.df <- triptans[triptans$agent=="placebo",]
# Meta-analyse placebo studies using fixed treatment effects
E0 <- ref.synth(ref.df, emax, synth="fixed")
names(E0)
# Meta-analyse placebo studies using random treatment effects
E0 <- ref.synth(ref.df, emax, synth="random")
names(E0)