model biodisperse;

#Modeling dose-response relation with a logistic regresion model
#using the overdispersed version of the model (page 441, GCSR)

var
  dose[N],                    # dosis recibida en cada grupo
  nanimals[N],                # numero de animales en cada grupo
  deaths[N],                  # numero de muertes por grupo
  prob[N],		      # "verdadera" probabilidad de muerte
  alpha, beta,                # parametros de la regresion
  e[N],			      # random effects model...
  tau,			      # precision del modelo
  ld50;			      # parametro que define la dosis a la
			      # cual hay un 50% de tasa de supervivencia 

data in "bio.dat";
inits in "bio.2.in";
{
  alpha ~ dnorm(0.0,1.0E-06); # intercepto
  beta  ~ dnorm(0.0,1.0E-06); # pendiente
  tau ~ dgamma(1.0E-1,1.0E-1);
  for(i in 1:N) {                          
	e[i] ~ dnorm(0.0,tau);
	logit(prob[i]) <- alpha + beta* dose[i] + e[i];
        deaths[i] ~ dbin(prob[i],nanimals[i]);
  }
  ld50<-(-alpha/beta);
}