model{

 for (t in 1:Tmax){
   logR[t]~dnorm(mu[t],tau[p[t]])
     mu[t]<-log(a[p[t]])+log(S[t])-(b[p[t]]*S[t])
density[t]<-dnorm(logR[t],mu[t],tau[p[t]])
 }#t


#Priors
for (p in 1:3){
a[p]~dunif(0.001,20)
b[p]~dunif(0,0.01)

tau[p]<-pow(sig.tau[p],-2)
sig.tau[p]~dunif(0,4)
sig[p]<-sqrt(1/tau[p])
}
}