docs: protected static double calcPosterior(double f, double alpha, double mu, double sigma, double lambda) {
      final double pi = calcP_i(f, mu, sigma);
      final double qi = calcQ_i(f, lambda);
      return (alpha * pi) / (alpha * pi + (1.0 - alpha) * qi);
    }