This function tries Newton's method. If not working, then use bisection.
Usage
vga_pois_solver(
init_val,
x,
s,
beta,
sigma2,
maxiter = 1000,
tol = 1e-05,
method = "newton"
)
Arguments
- init_val
initial value for posterior mean
- x, s
data and scale factor
- beta, sigma2
prior mean and variance. Their length should be equal to n=length(x)