From f21723ec23a6e5960d5a94660193298cd8f0d5d3 Mon Sep 17 00:00:00 2001
From: Thibaut Horel <thibaut.horel@gmail.com>
Date: Tue, 18 Aug 2015 22:15:45 -0700
Subject: WIP

---
 experiments/new.py | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)
 create mode 100644 experiments/new.py

(limited to 'experiments/new.py')

diff --git a/experiments/new.py b/experiments/new.py
new file mode 100644
index 0000000..4274dd7
--- /dev/null
+++ b/experiments/new.py
@@ -0,0 +1,23 @@
+from math import log, exp
+
+T = 100
+N = 100
+
+
+def kernel(t, mu):
+    return mu * exp(-mu * t)
+
+
+def base_rate(t, lamb):
+    return lamb
+
+
+def ll(crimes, weights, mu, lamb):
+    r = 0
+    for i, crime in enumerate(crimes):
+        t, v = crime
+        a = sum(weights[(u, v)] * kernel(t - s, mu) for s, u in crimes[:t])
+        r += log(base_rate(t, lamb) + a)
+    for j in range(N):
+        a = sum(weights[(u, v)] * kernel(T - s, mu) for s, u in crimes)
+        r -= log(a)
-- 
cgit v1.2.3-70-g09d2