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')
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
