サッカーのスコアのモデリング

Dixon、Coles（1997）では、彼らは（4.3）の2つの修正された独立したポアソンモデルの最尤推定を使用して、サッカーのスコアをモデル化しました。

パッケージを使用せずにアルファとベータ、およびホームエフェクトパラメーター（pg。274、表4）を "再現"するためにRを使用しようとしています（通常の独立したポアソンモデルを使用しても問題ありません）。bivpoisパッケージを使用してみましたが、パラメーターを変更する方法がわかりません。

誰かがRコードを使ってデータをモデル化できるように手助けしていただければ幸いです。英語プレミアリーグのシーズン2012/13のホームチームとアウェーチームのスコア。

あなたが読んでいる論文は、Maher（1982）によって記述されているように、攻撃と防御のパラメーターを参照するために暗黙的にとを使用しています。$\alpha_i$$\beta_i$

主な違いは、Maherがチームごとに4つのパラメーター（ホームアタック、ホームディフェンス、アウェイアタック、アウェイディフェンス）を使用するのに対し、ディクソンとコールズはアタックおよびディフェンスパラメーターと別のパラメーターを使用してホームアドバンテージを表すことです。

ポアソン分布のMLEは単純です： $\lambda_{MLE}= \frac{1}{n} \sum_{i=1}^n k_i$

..ポアソン分布への変更を再現する限り（簡単に見ると、時間依存性と二変量の両方になっていることがわかります）、誰もあなたのためにそうすることはないでしょう。実際に意味のあるツールを使用する方がずっと良いです。

二変量ポアソンは必要ありません。独自の関数を定義してから、optimのような一般的な最適化スクリプトを使用できます。

bivpoisRパッケージを自由に参照してください。私の以前のプロジェクトは対角二変量毒を適用しました。https://github.com/scibrokes/odds-modelling-and-testing-inefficiency-of-sports-bookmakersを参照してください。