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

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

# model file model-012f.bug fits an intercept and gradient model with an
# additional squared terms make the model a sextic fit corresponding R file
# model-012.r

data{
  n <- length(weight)
  day2 <- day^2
  day3 <- day^3
  day4 <- day^4
  day5 <- day^5
  day6 <- day^6
}


model{
  
  for (i in 1:n) {
    weight[i] ~ dnorm(mu[i], tau)
    mu[i] <- b0 + b1 * day[i] + b2 * day2[i] + b3 * day3[i] + b4 * day4[i] + 
      b5 * day5[i] + b6 * day6[i]
  }
  
  tau ~ dgamma(0.05, 0.05)
  b0 ~ dnorm(0, 1e-04)
  b1 ~ dnorm(0, 1e-04)
  b2 ~ dnorm(0, 1e-04)
  b3 ~ dnorm(0, 1e-04)
  b4 ~ dnorm(0, 1e-04)
  b5 ~ dnorm(0, 1e-04)
  b6 ~ dnorm(0, 1e-04)
}