我已经看到许多使用R完成未合并最大似然估计的教程示例。但是，到目前为止，我还没有在 R 中找到任何合并最大似然估计的示例。虽然我当然可以尝试不同的软件，但我仍然有兴趣看看它是否可以在R中轻松实现。
例如，假设我正在测量一些放射性同位素的衰变，并且我已经获得了超过 250 个时间段的事件数量，每个时间段的持续时间为 $\Delta t$。每个 bin 中的事件数量可以通过泊松统计建模，每个 bin 的预期数量为 $$ \bar{N}_{bin} = \Delta t N_0\lambda e^{-\lambda t_{bin}} $ $ 所以分箱最大似然函数的对数由 $$ \sum_{bin} \left( N_{bin} \ln(\bar{N}_{bin}) - \bar{N}_{bin} -\ln(N_{bin}!) \right) $$ 我要估计的参数是$\lambda$和$N_0$。
所以我的问题是这是否可以使用R相对简单地完成。谢谢！