Gibbs sampling for a hierarchical regression model
Implements Gibbs sampling for estimating a two-way table of means under a order restriction.
ordergibbs(data,m)
data |
data matrix with first two columns observed sample means and sample sizes |
m |
number of cycles of Gibbs sampling |
matrix of simulated draws of the normal means where each row represents one simulated draw
Jim Albert
data(iowagpa) m=1000 s=ordergibbs(iowagpa,m)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.