1. Função adaptabilidade e estabilidade bayesiana

#Leitura de dados# 
#Esse arquivo contém as médias dos genótipos em cada ambiente e o índice ambiental.
dados<-read.table("dadosR.txt",header=TRUE)
attach(dados)

#####################################################################

##Genótipo 
#Regressão considerando informação a priori
post1c<-MCMCregress(X1~i, data=dados, burnin=10000, mcmc= 100000, thin=2, b0=c(1515.90,1.09),
B0= matrix(c(0.000015639,0,0,33.33),ncol=2,nrow=2,byrow=TRUE),c0=1.015324 , d0=2522.665218,
marginal.likelihood="Laplace")
summary(post1c)

#Regressão sem informação a priori
post1s<-MCMCregress(X1~i,data=dados, burnin=10000, mcmc= 100000, thin=2, b0=c(0,0),
B0= matrix(c(0.0000001,0,0,0.0000001),ncol=2,nrow=2,byrow=TRUE),c0=0.00001 , d0=5000, 
marginal.likelihood="Laplace")
summary(post1s)