SaveOutput {Rtwalk}R Documentation

Save the twalk MCMC output to a file.

Description

Save the twalk MCMC output to a file.

Usage

SaveOutput( info, file, pars=1:(info$dim), from=1, to=info$Tr,
	row.names=FALSE, col.names=paste("X", pars), ...)

Arguments

info

as returned from Runtwalk.

file

name of file to write results to (of not "").

from

iteration number to start saving (from=0 begings at initialization point).

to

last iteration to save.

pars

parameters to save (defaults to all).

row.names, col.names, ...

parameters passed to write.table.

Value

write.table object.

Author(s)

J Andres Christen (CIMAT, Guanajuato, MEXICO).

References

Christen JA and Fox C (2010). A general purpose sampling algorithm for continuous distributions (the t-walk)., Bayesian Analysis, 5 (2), 263-282. URL: http://ba.stat.cmu.edu/journal/2010/vol05/issue02/christen.pdf

See Also

Runtwalk for running the twalk.

Examples

  

#### We first load the twalk package:
library(Rtwalk)



#### A ver simple inline example, 4 independent normals N(0,1):
######        dimension,  num of it, -log of objective function besides a const, support,
info <- Runtwalk( dim=4,  Tr=1000,  Obj=function(x) { sum(x^2)/2 }, Supp=function(x) { TRUE },
	x0=runif(4, min=20, max=21), xp0=runif(4, min=20, max=21)) 
####  and two (intentionally bad) initial points


### Save the twalk MCMC output as columns in a table
SaveOutput( info, file="Tsttwalk.dat")


[Package Rtwalk version 1.8.0 Index]