powered by
Implements Gibbs sampling for estimating a two-way table of means under a order restriction.
ordergibbs(data,m)
matrix of simulated draws of the normal means where each row represents one simulated draw
data matrix with first two columns observed sample means and sample sizes
number of cycles of Gibbs sampling
Jim Albert
data(iowagpa) m=1000 s=ordergibbs(iowagpa,m)
Run the code above in your browser using DataLab