Can i sample without any observation but only potential?

Questions
4

I’ve written code for what I think you want to do in PyMC3. I apologize that I can’t answer with PyMC2 because I don’t know much about how to use it.

import pymc3 as pm
import numpy as np

gs = np.random.uniform(size=(4, 3))

with pm.Model() as model:
    seed  = pm.Normal('seed', 0, 1, shape=[gs.shape[1], 1])
    label = pm.Deterministic('label',pm.math.gt(pm.math.dot(gs,seed),0))
    value = pm.Potential('value',pm.math.sum(pm.math.dot(label.T,gs)))
    step  = pm.Metropolis()
    trace = pm.sample(step=step)
2 Likes