Я ищу реализацию Inverse Incomplete Beta Function, возможно, уже написанную на C ++ или простую для реализации самостоятельно.Однако мне нужно, чтобы это было БЫСТРО!Например, я собираюсь выполнить это во внутреннем цикле оптимизатора, так что, надеюсь, это займет менее пары сотен тактов.
Здесь уже есть пара потоков, но в этом случае я готов отбросить много точности для скорости.Кроме того, домен несколько ограничен, так как я использую только целочисленные значения для a и b.
Дополнительные сведения об этой проблеме: я даю целое число испытаний n и целое число k <= nиз этих испытаний, которые были успешными.Я предполагаю, что фоновое распределение для основной вероятности успешного испытания является равномерным в [0,1], поэтому, учитывая, что я видел некоторое количество испытаний и успехов, мое последующее распределение должно быть бета-распределением.В байесовской модели я, по сути, пытаюсь найти p-й процентиль вероятных основных вероятностей. </p>
Опять же, мне не нужно, чтобы это было чрезвычайно точно, просто быстро.Я могу справиться с погрешностью до +/- 1%.Тем не менее, это не может быть очень неточным для небольших чисел: мои входные данные варьируются от нуля до десятков тысяч.
Заранее спасибо!Если понадобится какое-либо разъяснение, дайте мне знать.