inst/models/model-022f.bugs.R

In jaggR: Supporting Files and Functions for the Book Bayesian Modelling with 'JAGS'

# model-022f.bug - model for means of grouped observations - corner constrained
# parameterisation and paper weights and angle as covariates corresponding R file
# model-022.r


data{
  n <- length(distance)
  n.papers <- max(paper)
  n.angles <- max(angle)
  n.designs <- max(design)
}


model{
  
  for (i in 1:n) {
    distance[i] ~ dnorm(mu[i], tau)
    mu[i] <- b0 + b1[paper[i]] + b2[angle[i]]
  }
  
  for (j in 2:n.papers) {
    b1[j] ~ dnorm(0, 1e-07)
  }
  b1[1] <- 0  # corner constrain the paper weight
  
  for (k in 2:n.angles) {
    b2[k] ~ dnorm(0, 1e-07)
  }
  b2[1] <- 0  # corner constrain the launch angle
  
  # dummy iterations
  for (m in 2:n.designs) {
    b3[m] <- design[m]
  }
  
  tau ~ dgamma(0.001, 0.001)
  b0 ~ dunif(0, 10000)
}