model bio;

#Modeling dose-response relation with a logistic regresion model


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
  ld50;			      # parametro que define la dosis a la
			      # cual hay un 50% de tasa de supervivencia 

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