| JAGSrun {bayesmix} | R Documentation |
Calls jags for MCMC sampling.
JAGSrun(y, prefix = yname, model = BMMmodel(k = 2),
control = JAGScontrol(variables = c("mu", "tau", "eta")), tmp = TRUE,
cleanup = TRUE, jags = getOption("jags.exe"), ...)
y |
a numeric vector. |
prefix |
character prefix of files. |
model |
object of class JAGSmodel or output from BMMmodel. |
control |
specification of control by a JAGScontrol object. |
tmp |
logical: shall the files be written in a temporary directory. |
cleanup |
logical: shall the created files be removed. |
jags |
string indicating location of jags executable. |
yname |
a character string with the actual y argument name. |
... |
further parameters handed over to BMMmodel where it
is used for the function specifiying the initial values, e.g.,
initsFS. |
If an error occurs when runing jags, the created files are not
removed. This function is a wrapper calling JAGSsetup,
JAGScall and JAGSread.
Returns a jags object with components
call |
the matched call. |
results |
results read in from ``jags.out'' if run was successful or from ``jags.dump'' if an error occurred. |
model |
a JAGSmodel object. |
variables |
vector containing the names of the monitored variables. |
data |
a numeric vector. |
Bettina Gruen
JAGSsetup, JAGScall,
JAGSread, BMMmodel, initsFS
data(fish)
prefix <- "fish"
variables <- c("mu","tau","eta")
k <- 3
modelFish <- BMMmodel(k = k, priors = list(kind = "independence",
parameter = "priorsFish", hierarchical = "tau"))
controlFish <- JAGScontrol(variables = c(variables, "S"), draw = 100, seed = 1)
## Installation of JAGS necessary for applying these functions
if (haveJAGS()) {
z1 <- JAGSrun(fish, prefix, model = modelFish, initialValues = list(S0 = 2),
control = controlFish, cleanup = TRUE, tmp = FALSE)
zSort <- Sort(z1, "mu")
BMMposteriori(zSort)
}
data(darwin)
prefix <- "darwin"
k <- 2
modelDarwin <- BMMmodel(k = k, priors = list(kind = "independence",
parameter = "priorsUncertain"), aprioriWeights = c(1, 15),
no.empty.classes = TRUE, restrict = "tau")
## Installation of JAGS necessary for applying these functions
if (haveJAGS()) {
z2 <- JAGSrun(darwin, prefix, model = modelDarwin, control =
JAGScontrol(variables = variables, draw = 3000, burnIn = 1000,
seed = 1), cleanup = TRUE, tmp = FALSE)
plot(z2, variables = "mu")
}