modelDisplay('window')
modelPrecision(6)
modelCheck('MhBetamod.txt')
modelData('MhBetadat.txt')
modelCompile(3)
modelInits('MhBetaint1.txt', 1)
modelInits('MhBetaint2.txt', 2)
modelInits('MhBetaint3.txt', 3)
modelGenInits()
samplesSet(N)
samplesSet(mu)
samplesSet(tau)
modelUpdate(300000, TRUE)
samplesStats("*")
samplesHistory("*")
samplesTrace("*")
samplesDensity("*")
modelSaveLog('MhBetaLog')
samplesCoda(*,MhBetacoda)
modelQuit()