display('log')
check('Test/NegBin_model.txt')
data('Test/NegBin_data.txt')
compile(1)
inits(1,'Test/NegBin_inits.txt')
update(500)
set(theta)
set(alpha)
set(beta)
update(2000)
stats(*)
history(*)
density(*)
autoC(*)
save('NegBinLog')