# observation
  for (i in 1:N.mouse) { # age = 1
    for (age in 1:Age.max) {
      Y[i, age] ~ dnorm(y[i, age], Tau.obs)
      log(y[i, age]) <- log.y[i, age]
    }
  }
  # state space
  for (i in 1:N.mouse) {
    log.y[i, 1] ~ dnorm(0.0, tau[1]) # age = 1
    r[i, 1] ~ dnorm(0.0, Tau.obs)
    for (age in 2:Age.max) {
      log.y[i, age] ~ dnorm(m[i, age], Tau.obs)
      r[i, age] ~ dnorm(0.0, tau[2]) # individual.rnd
      m[i, age] <- (
        log.y[i, age - 1] 
        + age.effect[age] 
        + r[i, age]
      )
    }
  }
  # parameter
  age.effect[1] ~ dnorm(0.0, Tau.obs)
  for (age in 2:Age.max) {
    age.effect[age] ~ dnorm(0.0, Tau.Noninformative)
  }
  for (k in 1:N.tau) {
    tau[k] <- 1.0 / (s[k] * s[k])
    s[k] ~ dunif(0, 1.0E+4)
  }
}