Now working main_Test.r routines. Yet to fully test with range of parameters and analyse output. Added PowerGLM test script comparing Poisson samples of different sizes.