rm(list=ls())

getwd()

#setwd("pathway")



reg_data<-read.table("wg_priors.txt",header=T)

# reg_data$
  
names(reg_data)

cuadro<-array(0,c(8,3))

colnames(cuadro)<-c("b0","b1","b2")

names<-c("DHt122","DHt155","Dh187","DHt222","DHt264","DHt332","DHt412","KK")

rownames(cuadro)<-names

par(mfrow=c(3,3))

co<-c(1,3,5,7,9,11,13,15)



for(i in 1:8) {
  
  plot(reg_data[,co[i]],reg_data[,co[i]+1],main=names[i],xlab="number",ylab="height")
  
  y<-reg_data[,co[i]+1]
  
  x<-reg_data[,co[i]]
  x <- scale(x)
  
  x2<-reg_data[,co[i]]^2
  x2 <- scale(x2)
  
  model<-lm(y~x+x2)
  
  cuadro[i,]<- as.numeric(model$coeff)	
  
}

cuadro
colSums(cuadro)/8
sqrt(var(cuadro[,1]))
sqrt(var(cuadro[,2]))
sqrt(var(cuadro[,3]))