R/EB_poisson_mean_routines.R
pois_mean_GP_opt_obj.Rd
calculate objective function
pois_mean_GP_opt_obj(theta, x, s, beta, sigma2, n)