mcmc.list
mcmclist2coda.Rd
This function writes a coda file from an object of class mcmc.list
.
Note that only first entry (i.e. one chain) will be processed.
mcmclist2coda(mcmclist, name, coda.digits=5)
mcmclist | An object of class |
---|---|
name | Name of the coda file to be written |
coda.digits | Number of digits after decimal in the coda file |
The coda file and a corresponding index file are written into the working directory.
if (FALSE) { ############################################################################# # EXAMPLE 1: MCMC estimation 2PNO dataset Reading ############################################################################# data(data.read) # estimate 2PNO with MCMC with 3000 iterations and 500 burn-in iterations mod <- sirt::mcmc.2pno( dat=data.read, iter=3000, burnin=500 ) # plot MCMC chains plot( mod$mcmcobj, ask=TRUE ) # write sampled chains into codafile mcmclist2coda( mod$mcmcobj, name="dataread_2pl" ) }