model fans

{ for (i in 1:60)
      {t[i]~dweib(nu,lambda[i])
       log(lambda[i])<-mu + alpha[type[i]] + beta[speed[i]] + gamma[type[i],speed[i]]
       }

 nu~dgamma(1.5,1.5)
 
 mu~dnorm(-2,0.25)

 alpha[1]~dnorm(0,10)
 alpha[2]<-1-alpha[1]

 beta[1]<- - delta[1]+delta[2]
 beta[2]<- - 2*delta[2]
 beta[3]<-   delta[1]+delta[2]
 delta[1]~dnorm(0,10)
 delta[2]~dnorm(0,25)

 gamma[1,1]<- -eps[1]+eps[2]
 gamma[1,2]<-      -2*eps[2]
 gamma[1,3]<-  eps[1]+eps[2]
 gamma[2,1]<-  eps[1]-eps[2]
 gamma[2,2]<-       2*eps[2]
 gamma[2,3]<- -eps[1]-eps[2]
 eps[1]~dnorm(0,25)
 eps[2]~dnorm(0,100)

 }