Add unit testing with testthat and initial unit tests for negbin and poisson 2d simulations