Quadratic approximation to a posterior (eq. 28.5, p. 344)
Arguments
- dw
parameter offset from the posterior mode
- a
Hessian of the negative log posterior at the mode
Value
list(quadform, logratio, ratio, errorbars)
Examples
A <- rbind(c(2, 0.3), c(0.3, 1.5))
morie_postgapx(dw = c(0.2, -0.1), a = A)
#> $quadform
#> [1] 0.083
#>
#> $logratio
#> [1] -0.0415
#>
#> $ratio
#> [1] 0.9593493
#>
#> $errorbars
#> [1] 0.7179582 0.8290267
#>
