Rを使用してビン化されていない最尤推定を行うチュートリアルの例を数多く見てきました。ただし、現時点では、ビン化された最尤推定値の例を R でまだ見つけることができていません。確かに別のソフトウェアを試すことはできますが、 Rで簡単に実装できるかどうかを確認することにまだ興味があります。
例として、いくつかの放射性同位元素の崩壊を測定しているとします。250 時間ビンにわたる事象の数が提供され、それぞれの持続時間は $\Delta t$ です。各ビンのイベント数は、各ビンの予想数が $$ \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を使用して比較的簡単に実行できるかどうかです。ありがとう！