1

如何翻译这个使用来自GNU 科学图书馆的超几何分布函数的 C++ 程序

http://www.gnu.org/software/gsl/manual/html_node/The-Hypergeometric-Distribution.html

转换为使用Boost 库中类似函数的 C++ 程序？

#include <cstdlib>
#include <iostream>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_cdf.h>

int main(int argc,char *argv[]) {
  unsigned int n1 = 256;
  unsigned int n2 = 1583;
  unsigned int t = 300;
  unsigned int k = 40;
  std::cout << gsl_ran_hypergeometric_pdf(k, n1, n2, t)
            << std::endl
            << gsl_cdf_hypergeometric_P(k, n1, n2, t)
            << std::endl;
  return EXIT_SUCCESS;
}
4

1 回答 1

1

这是对 Boost 函数的翻译：

#include <cstdlib>
#include <iostream>
#include <boost/math/distributions/hypergeometric.hpp>
#include <boost/math/policies/policy.hpp>

int main(int argc, char *argv[]) {
  unsigned int n1 = 256;
  unsigned int n2 = 1583;
  unsigned int t = 300;
  unsigned int k = 40;
  boost::math::hypergeometric_distribution<double> hg_dist(n1, t, n1 + n2);
  std::cout << boost::math::pdf<double>(hg_dist, k)
            << std::endl
            << boost::math::cdf<double>(hg_dist, k)
            << std::endl;
  return EXIT_SUCCESS;
}

有关更多信息，请参阅 http://www.boost.org/doc/libs/1_52_0/libs/math/doc/sf_and_dist/html/math_toolkit/dist/dist_ref/dists/hypergeometric_dist.html

于 2012-11-07T17:12:58.253 回答