У меня есть некоторые интенсивные вычисления, связанные с распределением гипергеомотри c, и я хотел бы знать, существует ли реализация этой функции, которая поддерживает «реальные» вычисления векторизации? (Q1) Под «реальным» я подразумеваю: не использовать np.vectorize, который, кажется, реализует for l oop.
Я знаю, R реализация закона гипергеометрии c позволяет для такого рода вычислений, что делает код намного более эффективным.
Если это поможет: я пытаюсь продублировать алгоритм, разработанный в R, до Python, чтобы вычислить точные доверительные интервалы для гипергеометрии c параметры. Концепция была разработана в следующей статье: http://www.wright.edu/~weizhen.wang/paper/37-2015jasa_wang.pdf
Решение включает в себя итерации по всей возможной комбинации паремтеров для выбора подмножества, удовлетворяющего указанному c условию.
(Q2) В качестве альтернативы, можно ли было бы создать свой собственный udf, используя определение закона и метод np.choose, но примененный к векторам?
(Q3) в качестве окончательного варианта Можно ли вызвать мой скрипт в R из Python?
Спасибо!