diff options
| author | Thibaut Horel <thibaut.horel@gmail.com> | 2015-09-14 23:08:02 -0400 |
|---|---|---|
| committer | Thibaut Horel <thibaut.horel@gmail.com> | 2015-09-14 23:08:02 -0400 |
| commit | ab0b1f3cefedb35327a19ec1b6afd560bfdf802d (patch) | |
| tree | b777f3e2c0ac0e712d8c5faab5107b1d236e2c3a /experiments/new.py | |
| parent | 960676226862d2d68c7a9c04c56d4f8157803025 (diff) | |
| download | criminal_cascades-ab0b1f3cefedb35327a19ec1b6afd560bfdf802d.tar.gz | |
Import supplements and repo reorganization
Diffstat (limited to 'experiments/new.py')
| -rw-r--r-- | experiments/new.py | 23 |
1 files changed, 0 insertions, 23 deletions
diff --git a/experiments/new.py b/experiments/new.py deleted file mode 100644 index 4274dd7..0000000 --- a/experiments/new.py +++ /dev/null @@ -1,23 +0,0 @@ -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) |
