Displaying 1 result from an estimated 1 matches for "codasplit".
2007 Apr 29
0
Going from Coda Files to R2WinBUGS
...along the lines of c.bugs<-coda2bugs
(n.burnin=1000)
coda2bugs<-function(codafile="jags.out", indexfile="jags.ind",
n.chains=1,
n.iter=NA, n.burnin=NA, n.thin=1, DIC=FALSE, file.rm=T, ...){
require(R2WinBUGS)
#first, split up the coda file for R2WinBUGS
codaSplit(codafile, indexfile)
#get the parameter names
index.table<-read.table(indexfile)
varNames<-as.vector(index.table[,1])
#determine the n.iter
if(is.na(n.iter)){n.iter<-index.table[1,3]}
#you will need to put the n.burnin in yourself
#for the cypraea example, it is 100...