#Hierarchical stock-recruitment model of six Alaskan chinook stocks
#####################################################

model chinook;
{
######Create y variable, log(R/S)#####
for (i in 1:N) { y[i] <- log(R[i]/S[i]) }

######Model the y variable######
for (i in 1:N) { 
  mean.y[i] <- alpha[stock[i]] - beta[stock[i]]*S[i];
  y[i] ~ dnorm(mean.y[i],prec.y[stock[i]]);
}

######Priors#####
for (i in 1:6) {
  alpha[i] ~ dnorm(alpha.mu,alpha.prec);
  beta[i] ~ dnorm(0,0.000001);
  prec.y[i] ~ dgamma(0.001,0.001);
}
alpha[7] ~ dnorm(alpha.mu,alpha.prec);
alpha.mu ~ dnorm(0,0.000001);
alpha.sigma ~ dunif(0,100);
alpha.prec <- 1/(alpha.sigma*alpha.sigma);

###Silly trick to avoid data loading error###
for (i in 1:N) {  year[i] ~ dnorm(1980,0.001) }
}

###Data###
#list(N=97)
###Inits###
#list(alpha.mu=1,alpha.sigma=1,prec.y=c(1,1,1,1,1,1))