Tôi đang viết một lớp con scipy.stats._distn_infrastructure.rv_discrete
cho phân phối nhị thức beta có PMF là
Trong đó là hàm Beta. Việc triển khai CDF và SF hiện tại của tôi (chức năng sinh tồn, tương đương với 1 - CDF) là không chính xác; chiến lược tôi sử dụng sẽ tính toán giá trị kỳ vọng của cdf nhị thức đối với thành phần beta:
scipy.stats.beta.expect
phương thức không được vector hóa hoàn toàn (nó sẽ sụp đổ trên bất cứ thứ gì ngoài mảng float hoặc 0d).
PPF thậm chí còn tệ hơn - đó là một vòng lặp mạnh mẽ trên các số nguyên sao cho
Theo Wikipedia, chức năng sinh tồn cho phân phối nhị thức beta là
trong đó là hàm siêu bội tổng quát. Có cách nào hiệu quả để tính toán điều này trong Python, vì vậy tôi có thể xóa tham chiếu đến không? Ngoài ra, làm cách nào để đảo ngược hàm này để giải quyết cho cho ?beta.expect