x=c(94,96,94,95,104,106,108,113,115,121,131)
y=c(0.47,0.75,0.83,0.98,1.18,1.29,1.4,1.6,1.75,1.9,2.23)
require(MCMCpack)
plot(x,y)
posterior.out=MCMCregress(y~x)
summary(posterior.out)
plot(posterior.out)
beta=rnorm(100)
mean(beta<0.25)
mean(beta>0.25)
posterior.out