c++ - 如何在 C++ 中计算概率

Question

我正在转换一些 SAS code进入 C++，但目前停留在如何计算 probit 上功能。

我相信我已经找到了 probit 公式 ( http://en.wikipedia.org/wiki/Probit )，但我的统计知识有点生疏。所以基本上我需要有人将其转换为 C++，或者告诉我是否有其他名称下的相同方法。

Answer 1

使用 boost 数学库和您在维基百科文章中找到的 probit 的定义，我得出了以下结论。这里没有保证:-)

#include <boost/math/special_functions/erf.hpp>

namespace bm = boost::math;

template<typename T>
T probit(T p)
{
    T root_2 = sqrt(2);

    return root_2 * bm::erf_inv(2*p-1);
}

int main()
{
    double val = 0.9;
    double res = probit(val);
}