ordergibbs {LearnBayes} | R Documentation |
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)