aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorjeanpouget-abadie <jean.pougetabadie@gmail.com>2015-11-10 16:42:12 -0500
committerjeanpouget-abadie <jean.pougetabadie@gmail.com>2015-11-10 16:42:12 -0500
commit82b14ab18ccfd032f817979ec0ce74f274b69c3d (patch)
tree6ab8aee6e6c189c7957ee2b05b062c833ca959a5
parent4db9cdc26c52ecf596e552f569729278bc5dfd8e (diff)
downloadcascades-82b14ab18ccfd032f817979ec0ce74f274b69c3d.tar.gz
changing non-defined prior and source specification
-rw-r--r--simulation/bayes.py12
-rw-r--r--simulation/main.py6
2 files changed, 9 insertions, 9 deletions
diff --git a/simulation/bayes.py b/simulation/bayes.py
index 7d8567b..e2be6da 100644
--- a/simulation/bayes.py
+++ b/simulation/bayes.py
@@ -45,13 +45,13 @@ def mc_graph_setup(infected, susceptible, prior=None, *args, **kwargs):
# Container class for graph parameters
n_nodes = len(infected[0][0])
theta = np.empty((n_nodes,n_nodes), dtype=object)
- for i in xrange(n_nodes):
- for j in xrange(n_nodes):
- if prior is None:
+ if prior is None:
+ for i in xrange(n_nodes):
+ for j in xrange(n_nodes):
theta[i, j] = pymc.Beta('theta_{}{}'.format(i, j), alpha=1,
- beta=1)
- else:
- theta[i, j] = prior('theta_{}{}'.format(i, j), *args, **kwargs)
+ beta=1)
+ else:
+ theta = prior(theta=theta, *args, **kwargs)
# Container class for cascade realization
x = {}
diff --git a/simulation/main.py b/simulation/main.py
index e4b30f2..3e3de4b 100644
--- a/simulation/main.py
+++ b/simulation/main.py
@@ -108,15 +108,15 @@ def simulate_cascade(x, graph):
yield x, susc
-def uniform_source(graph):
+def uniform_source(graph, *args, **kwargs):
x0 = np.zeros(graph.shape[0], dtype=bool)
x0[nr.randint(0, graph.shape[0])] = True
return x0
def simulate_cascades(n, graph, source=uniform_source):
- for _ in xrange(n):
- x0 = source(graph)
+ for t in xrange(n):
+ x0 = source(graph, t)
yield simulate_cascade(x0, graph)