#The data
library(arm)
library(R2WinBUGS)

 ### Prob 11-77 Montgommery & Runger - Linear regression - full Bayesian Hierarchical
 test=matrix(scan("Prob-11-77-data.txt"),ncol=2,byrow=T)
N=dim(test)[1]
data=list(n=N,x=test[,1],y=test[,2])
inits = function(){list(alphamean=0,betamean=1,taue=1,alpha=0,beta=
1,taualpha=1,taubeta=1)}

#inits = function(){list(alpha=0,beta=1)}
parameters = c("alpha","beta","taue","alphamean","betamean","taualpha","taubeta")
 linreg.sim = bugs(data,inits=inits,parameters,"linear-reg-examp-prob11-77-full-bayes-bug.txt",n.chains=3,n.iter=10000)
 attach.bugs(linreg.sim)
 hist(alpha)
 hist(beta)
hist(taue)