comparepost: Line plots of selected parameters from multiple models.

In mbtyers/jagsplot: Diagnostic plotting tools for R2jags output

Description

Line plots of selected parameters from multiple models.

Usage

comparepost(postlist, years1 = NA, years2 = NA, pnames = NULL)

Arguments

postlist	A single, or list of, R2jags output objects or mcmc.list objects. The legend will be created from the names of the list objects.
years1	An optional vector of years to use for plotting labels, possibly useful for random effects
years2	An optional vector of years to use for plotting labels, possibly useful for random effects
pnames	A character list of which parameters to plot

Author(s)

Matt Tyers

Examples

# ------------ two competing R2jags models ----------- #

## Not run: 
library(R2jags)
N <- 1000
x <- 1:N
yreff <- rnorm(10,0,10)
yr <- rep(1:20, each=50)
epsilon <- rnorm(N, 0, 1)
y <- x + yreff[yr] + epsilon

cat('model {
 for (i in 1:N){
   y[i] ~ dnorm(mu[i], tau)
   mu[i] <- b0 + b1*x[i] + a.yr[yr[i]]
}
b0 ~ dnorm(0, .0001)
b1 ~ dnorm(0, .0001)
tau <- pow(sigma, -2)
sigma ~ dunif(0, 100)
for(j in 1:Nyr) {
  a.yr[j] ~ dnorm(0, tau.a)
}
tau.a <- pow(sigma.a, -2)
sigma.a ~ dunif(0, 100)
}', file="test.jags")

test.data <- list(x=x,y=y,N=N,yr=yr,Nyr=max(yr))
test.jags.out1 <- jags(model.file="test.jags", data=test.data, 
     parameters.to.save=c("b0","b1","a.yr","sigma","sigma.a"), n.chains=3, n.iter=10000)

cat('model {
 for (i in 1:N){
   y[i] ~ dnorm(mu[i], tau)
   mu[i] <- b0 + b1*x[i] + a.yr[yr[i]]
}
b0 ~ dnorm(0, .0001)
b1 ~ dnorm(0, .0001)
tau <- pow(sigma, -2)
sigma ~ dunif(0, 100)
for(j in 1:Nyr) {
  a.yr[j] ~ dnorm(0, 0.0001)
}
}', file="test.jags")

test.data <- list(x=x,y=y,N=N,yr=yr,Nyr=max(yr))
test.jags.out2 <- jags(model.file="test.jags", data=test.data, 
     parameters.to.save=c("b0","b1","a.yr","sigma"), n.chains=3, n.iter=10000)

# ------------- calling comparepost() ------------ #
par(mfrow=c(3,2))
comparepost(postlist=list(hierarchical=test.jags.out1,non_hier=test.jags.out2),years1=1998:2017)

### or just one posterior...
comparepost(test.jags.out1,years1=1998:2017)

## End(Not run)

mbtyers/jagsplot documentation built on May 22, 2019, 12:58 p.m.