#Script du modèle - 1 covariable


model{
# Y= teneur en protéines du blé
# covar_cr = variable explicative (SPAD ou INN) centrée réduite
# alpha=paramètres de la régression;   prec=paramètre de précision = 1/variance 
# m, s= moyenne et écart-type empirique de covar
#Modèle d’observations
for (i in 1:N) {         
 		Y[i] ~ dnorm(mu[i],prec)       
		#mu[i]<-alpha_nonc[1]+alpha_nonc[2]*covar[i] => Problème de convergence! 
		 mu[i]<-alpha[1]+alpha[2]*covar_cr[i] 
		#Calcul de la log_vraisemblance de l'observation i
  		loglik[i] <- -0.5*log(2*pi)+0.5*log(prec)-0.5*prec*pow(Y[i]-mu[i],2) 
}  
alpha_nonc[1]<- alpha[1]-alpha[2]*(m/s)
alpha_nonc[2]<- alpha[2]/s

 #Lois a priori
for (j in 1:2) {          
	alpha[j]~dnorm(0.0,1.0E-6)   
}
sd  ~ dunif(0,100)   
prec<- 1/pow(sd,2)
}


#Données - covariable INN + INN_pred
list(Y=c(1.19000E+01, 1.03000E+01, 1.24000E+01, 1.21000E+01, 1.24000E+01, 1.27000E+01, 1.15000E+01, 1.36000E+01, 9.90000E+00, 1.11000E+01, 1.06000E+01, 1.20000E+01, 1.28000E+01, 1.16000E+01, 1.12000E+01, 1.03000E+01, 1.29000E+01, 1.03000E+01, 1.15000E+01, 1.41000E+01, 1.26000E+01, 1.27000E+01, 1.19000E+01, 9.20000E+00, 1.25000E+01, 1.19000E+01, 1.13000E+01, 1.10000E+01, 1.00000E+01, 1.59000E+01, 1.06000E+01, 1.16000E+01, 1.14000E+01, 1.09000E+01, 1.46000E+01, 1.07000E+01, 9.40000E+00, 1.30000E+01, 1.25000E+01, 1.29000E+01, 1.11000E+01, 1.29000E+01, 1.32000E+01,     NA), covar_cr=c(1.26709E+00, -2.61237E-01, -2.61237E-01, 1.26709E+00, -1.02540E+00, 5.02926E-01, 1.26709E+00, 5.02926E-01, -1.78956E+00, -1.02540E+00, -1.02540E+00, -1.02540E+00, 5.02926E-01, -2.61237E-01, -2.61237E-01, -1.78956E+00, 5.02926E-01, -1.10182E+00, -8.72567E-01, 2.03125E+00, 1.03784E+00, 1.11426E+00, 2.73677E-01, -1.63673E+00, -7.96151E-01, 1.95483E+00, -5.66902E-01, 1.97261E-01, -6.43318E-01, 1.49634E+00, -6.43318E-01, 1.97261E-01, -1.08404E-01, -4.90486E-01, 3.50093E-01, 4.44281E-02, -8.72567E-01, -1.84821E-01, -7.19735E-01, 8.08591E-01, -4.14069E-01, 6.55758E-01, 1.80200E+00, 8.85007E-01), N=4.40000E+01, m=8.34186E-01, s=1.30862E-01, pi=3.14159E+00)



#Initialisation des chaînes
#Chaîne 1
list(
Y = c(
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA, 1.49),
alpha = c(
0.0,0.0),
sd = 1.0)


#Chaîne 2
list(
Y = c(
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA, 0.06),
alpha = c(
5.0,5.0),
sd = 0.1)


#Chaîne 3
list(
Y = c(
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,            NA,            NA,
            NA,            NA,            NA,	10.0),
alpha = c(
-5.0,-5.0),
sd = 10.0)