Rev Language Reference

srMaxTime - Maximum time stopping rule

Cause an MCMC run to terminate once the specified time has elapsed.

Usage

srMaxTime(RealPos maxTime, String unit)

Arguments

maxTime : RealPos (pass by value)
The maximum time to run.
unit : String (pass by value)
The unit in which we measure the maximum time.
Default : seconds
Options : seconds|minutes|hours

Example

Binomial example: estimate success probability given 7 successes out of 20 trials

r ~ dnExp(10) p := Probability(ifelse(r < 1, r, 1)) n <- 20 k ~ dnBinomial(n, p) k.clamp(7) mymodel = model(k)

moves = VectorMoves() moves.append( mvSlide(r, delta=0.1, weight=1) )

paramFile = “parameters.log”

monitors = VectorMonitors() monitors.append( mnModel(filename=paramFile, printgen=100, p) )

Stop when the five seconds have elapsed

stopping_rules[1] = srMaxTime(5, “seconds”)

Create the MCMC object

mymcmc = mcmc(mymodel, monitors, moves)

Begin the MCMC run

mymcmc.run(rules = stopping_rules)

Edit this page