library(nclbayes)

#Too cold
posterior1 = metropolis(N=10000, initial=0, a=0.6)
mcmcAnalysis(posterior, rows=1)

#Just right
posterior2 = metropolis(N=10000, initial=0, a=6)
mcmcAnalysis(posterior2, rows=1)

#Too hot
posterior3 = metropolis(N=10000, initial=0, a=60)
mcmcAnalysis(posterior3, rows=1)

#Bad starting value:
posterior5 = metropolis(N=10000, initial=20, a=6)
mcmcAnalysis(posterior5, rows=1)

#Bad starting value and too much dependence!
posterior6 = metropolis(N=10000, initial=30, a=0.1)
mcmcAnalysis(posterior6, rows=1)

#Run for longer, then process:
posterior6a = metropolis(N=100000, initial=30, a=0.1)
mcmcAnalysis(posterior6a, rows=1)

##Thinning needed and removal of burn-in necessary
posterior6b = mcmcProcess(input = posterior6a, burnin = 5000, thin = 100)
mcmcAnalysis(posterior6b, rows=1)

#Go back to "just right":
posterior2 = metropolis(N=10000, initial=0, a=6)
mcmcAnalysis(posterior2, rows=1)

##Run for longer, thinby around 20?
posterior2a = metropolis(N=100000, initial=0, a=6)
posterior2b = mcmcProcess(input = posterior2a, thin = 20)
mcmcAnalysis(posterior2b, rows=1)