model 

{for (i in 1:30)
     {loss[i]~dnorm(lossmean[i],tau)
      lossmean[i]<-alpha+beta[1]*(hard[i]-60)+beta[2]*(tens[i]-200)
      hard[i]~dnorm(muhard,tau.hard)
      tens[i]~dnorm(tensmean[i],tau.tens)
      tensmean[i]<-delta[1]+delta[2]*(hard[i]-60)
      }

 alpha~dnorm(150,0.000625)
 beta[1]~dnorm(0,0.0025)
 beta[2]~dnorm(0,0.0025)
 beta0<-alpha+60*beta[1]+200*beta[2]
 muhard~dnorm(60,0.0025)
 delta[1]~dnorm(200,0.0004)
 delta[2]~dnorm(0,1)

 tau.tens~dgamma(1,2000)
 tau.hard~dgamma(1,100)
 tau~dgamma(2,3200)

 }