Esempio n. 1
0
def test_chains(proposal_name, nprocs):
    """Runs the PTSampler ``test_chains`` test using the bounded normal
    proposal.
    """
    model = PoissonModel()
    proposal = _setup_proposal(proposal_name, model.params, model.prior_bounds)
    _test_chains(PoissonModel, nprocs, SWAP_INTERVAL, proposals=[proposal])
Esempio n. 2
0
def test_chains(name, nprocs):
    """Runs the PTSampler ``test_chains`` test using the adaptive proposal.
    """
    model = Model()
    # we'll just use the adaptive normal for one of the params, to test
    # that using mixed proposals works
    proposal = _setup_proposal(model, name,  params=[list(model.params)[0]])
    _test_chains(Model, nprocs, SWAP_INTERVAL, proposals=[proposal])
Esempio n. 3
0
def test_chains(nprocs, proposal_name):
    """Runs the PTSampler ``test_chains`` test using the adaptive normal
    proposal.
    """
    model = Model()
    # we'll just use the adaptive normal for one of the params, to test
    # that using mixed proposals works
    proposal = _setup_proposal(model, proposal_name,
                               params=[list(model.params)[0]])
    _test_chains(Model, nprocs, SWAP_INTERVAL, proposals=[proposal],
                 init_iters=STABILITY_DURATION-2)
Esempio n. 4
0
def test_chains(proposal_name, nprocs):
    """Runs the PTSampler ``test_chains`` test using the bounded normal
    proposal.
    """
    model = Model()
    # we'll just the bounded normal proposal for one of the parameters,
    # to test that using mixed proposals works
    param = list(model.params)[0]
    proposal = _setup_proposal(proposal_name, [param],
                               {param: model.prior_bounds[param]})
    _test_chains(Model, nprocs, SWAP_INTERVAL, proposals=[proposal])
Esempio n. 5
0
def test_chains(proposal_name, nprocs):
    """Runs the PTSampler ``test_chains`` test using the angular proposal.
    """
    model = AngularModel()
    proposal = _setup_proposal(proposal_name, model.params)
    _test_chains(AngularModel, nprocs, SWAP_INTERVAL, proposals=[proposal])