powered by
logp(x)
mcmc.normal
NR <- 100 NC <- 2 b <- mcmc.normal(rnorm(NC),mu=0,tau=0.0001) tau.y <- mcmc.gamma(runif(1),alpha=0.1,beta=0.1) logp(b) logp(tau.y)
Run the code above in your browser using DataLab