factor potential example

jsalvatier · jsalvatier · commit 5f3e90649b66 · 2013-04-26T00:51:44.000-07:00
diff --git a/examples/factor_potential.py b/examples/factor_potential.py
@@ -0,0 +1,10 @@
+from pymc import * 
+
+with Model() as model: 
+    x = Normal('x', 1,1)
+    model.AddPotential(-x**2)
+
+    
+    h = find_hessian()
+    step = Metropolis(h)
+    trace = sample(3000, step)
