This function approximates the Gibbs sample mode.
Value
The (approximated) mode.
Examples
samples <- oeli::rmixnorm(
n = 1000, mean = matrix(c(-2, 2), ncol = 2),
Sigma = matrix(c(1, 1), ncol = 2), proportions = c(0.7, 0.3)
)
hist(samples)
mean(samples) # expected: 0.7 * (-2) + 0.3 * 2 = -0.8
#> [1] -0.8136716
mode_approx(samples) # expected: -2
#> [1] -2.061629
