[R-sig-Geo] Openbugs doesn't seen to want to thin.. ?

James Rooney ROONEYJ4 at tcd.ie
Mon Aug 26 12:22:16 CEST 2013


Hi all,

Doing some spatial smoothing using the BYM model from the ASDAR book.
I was using Winbugs but had lots of troubles so switched to Openbugs.
Unfortunately having a new problem running out of memory.
I have tried to fix it with thinning - however Openbugs does not seem to want to thin.
My Openscript is below (generated by R2Openbugs).
But no matter what I set the thin level too it still gives me the full 40,000 iterations in the CODA file.
This is a problem because Openbugs is 32bit and above 40k iterations my model gets bigger than 4GB.

Anyone know why it is not thinning ?
Thanks,
James

Script (FYI - I've clipped the pathnames out). This runs fine but the 3 in the modelUpdate() command has no effect - no matter what I change it to it saves all 40,000 samples.

modelDisplay('log')
modelCheck('model.txt')
modelData('data.txt')
modelCompile(2)
modelSetRN(1)
modelInits('inits1.txt',1)
modelInits('inits2.txt',2)
modelGenInits()
modelUpdate(1000,3,1000)
samplesSet(theta)
samplesSet(alpha)
samplesSet(tau)
samplesSet(precv)
samplesSet(deviance)
summarySet(theta)
summarySet(alpha)
summarySet(tau)
summarySet(precv)
summarySet(deviance)
dicSet()
modelUpdate(40000,3,40000)
samplesCoda('*', '//')
summaryStats('*')
dicStats()
samplesHistory('*')
modelSaveLog('log.odc')
modelSaveLog('log.txt')


More information about the R-sig-Geo mailing list