Rev Language Reference

srMaxTime - Maximum time stopping rule

Terminates an MCMC run 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)

See Also

Edit this page