Pemodelan untuk skor sepakbola

Dalam Dixon, Coles ( 1997 ), mereka telah menggunakan estimasi kemungkinan maksimum untuk dua model Poisson independen yang dimodifikasi dalam (4.3) untuk memodelkan skor dalam sepak bola.

Saya mencoba menggunakan R untuk "mereproduksi" alfa dan beta serta parameter efek rumah (hal 274, Tabel 4) tanpa menggunakan paket apa pun (menggunakan model Poisson independen yang biasa juga baik-baik saja). Saya telah mencoba menggunakan bivpoispaket tetapi saya tidak yakin tentang cara memodifikasi parameternya.

Saya akan sangat menghargai jika ada yang bisa membantu saya dengan kode R untuk memodelkan data - Skor dari tim tuan rumah dan tandang untuk Musim 2012/13 di Liga Primer Inggris.

Makalah yang Anda baca secara implisit menggunakan dan β i untuk merujuk pada parameter serangan dan pertahanan seperti yang dijelaskan oleh Maher (1982) .$\alpha_i$$\beta_i$

Perbedaan utama adalah bahwa Maher menggunakan empat parameter untuk setiap tim (serangan rumah, pertahanan rumah, serangan tandang dan pertahanan tandang) sementara Dixon dan Coles menggunakan parameter serangan dan pertahanan serta parameter lain untuk mewakili keunggulan kandang.

$\lambda_{MLE}= \frac{1}{n} \sum_{i=1}^n k_i$

.. sejauh mereproduksi perubahan mereka pada distribusi Poisson (sekilas memberitahu saya itu telah menjadi tergantung waktu dan bivariat), saya ragu ada yang akan melakukan itu untuk Anda. Anda jauh lebih baik menggunakan alat yang benar-benar masuk akal.

Anda tidak perlu Poisson bivariat. Anda dapat mendefinisikan fungsi Anda sendiri, dan kemudian menggunakan skrip optimisasi umum seperti optim.

Anda merasa bebas untuk merujuk ke bivpoispaket R. Proyek saya sebelumnya menerapkan racun bivariat diagonal. karena Anda dapat merujuk ke https://github.com/scibrokes/odds-modelling-and-testing-inefficiency-of-sports-bookmakers .