d3443536f3d4e39d5e0b8f90e1f9d0e25af77713,pymc/step_methods/metropolis.py,BinaryMetropolis,astep,#BinaryMetropolis#Any#Any#,159

Before Change

        docstring for astep

        // Calculate probabilities of each value
        logp_true = logp([True]*len(q0))
        logp_false = logp([False]*len(q0))

        p_true = exp(logp_true)
        p_false = exp(logp_false)

        // Stochastically set value according to relative
        // probabilities of True and False
        return (random(len(q0)) < p_true / (p_true + p_false)).astype(int)

After Change

        switch_locs = where(rand_array<p_jump)
        q[switch_locs] = True - q[switch_locs]

        q_new = metrop_select(logp(q) - logp(q0), q, q0)

        return q_new

def subst(a, val, index):
    ap = a.copy()
    ap[index]=val

Instances