model{

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


#Priors
for (p in 1:3){
a[p]~dunif(0.001,20)
b[p]~dunif(0,0.01)
}
tau<-pow(sig.tau,-2)
sig.tau~dunif(0,4)
sig<-sqrt(1/tau)

}