我已经看到很多使用 R 完成 unbinned 最大似然估计的教程示例。但是，到目前为止，我还没有能够在R中找到 binned 最大似然估计的任何示例。虽然我当然可以尝试不同的软件，但我仍然有兴趣看看它是否可以很容易在 R 。
中实施作为一个例子，我说我正在测量一些放射性同位素的衰变，并且我已经获得了超过250个时间段的事件数量，每个时间段为$ \ Delta t $。每个仓中的事件数量可以通过泊松统计建模，每个仓的预期数量 $$ \ bar {N} _ {bin} = \ Delta t N_0 \ lambda e ^ { - \ lambda t_ {bin}} $$ 所以binned最大似然函数的对数由下式给出 $$ \ sum_ {bin} \ left（N_ {bin} \ ln（\ bar {N} _ {bin}） - \ bar {N} _ {bin} - \ ln（N_ {bin}！）\ right） $$ 我想估计的参数是$ \ lambda $和$ N_0 $。
所以我的问题是，使用 R 是否可以相对简单地完成。谢谢！