calculate.moller.hetero: calculate.moller.hetero
In georgeshirreff/multiNe: Multiple Ne estimators

Approximates the posterior distribution of Ne from a single genealogy at a regular grid of points using INLA package

1 2	calculate.moller.hetero(coal.factor, s, event, lengthout, prec_alpha = 0.01, prec_beta = 0.01, E.log.zero = -100, alpha = NULL, beta = NULL)

`coal.factor`	is a vector with coalescent times in increasing order
`s`	sampling times and coalescent times in increasing order
`event`	and indicator vector with 1 for coalescent times and 0 for sampling times that correspond to the s vector
`lengthout`	number of grid points
`prec_alpha`	alpha gamma hyperparameter of precision parameter for GP prior
`prec_beta`	beta gamma hyperparameter of precision parameter for GP prior
`E.log.zero`	internal
`alpha`	TODO
`beta`	TODO