Initial commit. Poisson code works, Custom code (and R script) incomplete/untested