Skip to content

Instantly share code, notes, and snippets.

View euxoa's full-sized avatar

Janne Sinkkonen euxoa

  • Helsinki, Finland
View GitHub Profile
data {
int<lower=0> N; // rows
int<lower=1> D; // delays, first is 0, last D-1
int<lower=0> T; // days, first is 1
int<lower=0> counts[N];
int<lower=1> t[N];
int<lower=0, upper=D-1> dt0[N];
int<lower=0, upper=D-1> dt1[N];
real<lower=0> dir_prior;
@euxoa
euxoa / shootings.stan
Created April 21, 2019 16:40
A simple Poisson state-space model with seasonality and overdispersion
data {
int<lower=1> n; // n = 2632
int y[n]; // number of shootings per day
}
transformed data {
int<lower=1,upper=365> yday[n];
for (i in 1:n) yday[i] = i % 365 + 1;
}
parameters {
vector[n] mu_innovations;
library(dplyr)
# install.packages("rgeos", type="source")
# install.packages("rgbif", type="source")
# library(rgbif)
# d <- occ_search(orderKey=797, datasetKey="a8d08280-1def-11de-be11-b8a03c50a862", year="2013", limit=2000)
library(RJSONIO)