pm.wbugs: Run Pella-Masuda model via WinBUGS

In bbolker/mixstock: Functions for Mixed Stock Analysis

Run Pella-Masuda model via WinBUGS

Description

Sets up the Pella-Masuda model and passes it to WinBUGS via R2WinBUGS

Usage

pm.wbugs(x,
n.iter=20000,  n.burnin=floor(n.iter/2),
n.chains=x$R,
n.thin=max(1, floor(n.chains * (n.iter - n.burnin)/1000)),
...)

Arguments

x	a mixstock data object
n.iter	Total length of each chain
n.burnin	Number of burn-in iterations
n.chains	Number of chains (default, number of sources)
n.thin	thinning rate. Must be a positive integer. Set 'n.thin' > 1 to save memory and computation time if 'n.iter' is large. Default is 'max(1, floor(n.chains * (n.iter-n.burnin) / 1000))' which will only thin if there are at least 2000 simulations.
...	other arguments to bugs

Value

a BUGS object (as returned from R2WinBUGS); can be converted to CODA format using as.mcmc.bugs

Note

tmcmc is in general be much more convenient and efficient than pm.wbugs: pm.wbugs is included for completeness and testing of WinBUGS methods.

Author(s)

Ben Bolker

See Also

tmcmc

Examples

data(bolten98)
## Not run: 
bolten98.wbugs <- pm.wbugs(bolten98,tot=1000,clearWD=TRUE)
bolten98.wbugs

## End(Not run)

bbolker/mixstock documentation built on July 23, 2024, 12:18 p.m.