post<-read.table('sum.txt',header=TRUE)
 cor1<-read.table('cor1.txt')
 cor2<-read.table('cor2.txt')
 
 c2<-1000
 
 #c2<-exp(700)
 
 cor<-matrix(cor1[,3],nrow=3)
 
 cor<-cbind( cor, cor2[,3])
 cor<-rbind( cor, c(cor2[,3],1))

m<-post[c(1:3,7),2]
sd<-post[c(1:3,7),3]


-0.5*4*log(2*3.14)+sum(log(sd))+log(det(cor))



beta <- post[c(1:3),2]
tau  <- post[7,2]
s    <- 1/sqrt(post[7,2])

eta<- x4%*%beta


0.5*4*log(2*3.14)+sum(log(sd))+0.5*log(det(cor))+sum( dnorm( y, eta, s, log=TRUE ) )+ sum( dnorm( beta, 0, sqrt(c2)*s, log=TRUE ) )+ dgamma( tau, 0.001, 0.001, log=TRUE ) 


0.5*4*log(2*3.14)+sum(log(sd))+0.5*log(det(cor))
sum( dnorm( y, eta, s, log=TRUE ) )
sum( dnorm( beta, 0, sqrt(1000)*s, log=TRUE ) )
dgamma( tau, 0.001, 0.001, log=TRUE )