假设输出为n个总数,其中k是从速率为a的指数随机变量中采样的,n-k是从速率为b的指数随机变量中采样的,约束条件是0<a≤b和0≤k≤n,但a、b和k都是未知的。我的目标是估计a或b(不关心k,我不需要同时估计a和b:两种方法中的一种就可以了)。我想在理想的情况下用Python来做,但是我对任何自由软件都是开放的。我的第一种方法是使用sklearn.optimize优化一个似然函数,其中,对于数据集中的每一个数,我计算P(X=x)的指数率为a,对于指数的速率为b,只需选择两者中较大的一个:
from