R 足球得分建模（迪克森和科尔）

在Dixon，Coles（）中，他们使用了（4.3）中两个修正的独立泊松模型的最大似然估计来模拟足球中的得分

我试图使用R来“重现”α和β以及home effect参数（第274页，表4），而不使用任何包（使用通常的独立泊松模型也可以）。我尝试过使用bivpois包，但我不确定如何修改它的参数

如果有人能帮我用R代码来模拟2012/13赛季英超主客场球队的得分数据，我将不胜感激。所以，基本上我需要帮助，用R中的optim函数对方程4.3或4.5进行编码

泊松分布密度（仅1个独立泊松模型）

正态泊松似然函数

ll.poisson <- function(par, y) {
(alpha*beta*home) <- exp(par)
out <- sum(y * log((alpha*beta*home))) - length(y) * (alpha*beta*home)
return(out)
}

---

您好，根据一般经验，有关StackOverflow的问题应该显示出一些工作和努力。具体来说，对代码的请求应该明确地有一些尝试的指示，特别是您需要的帮助。请随意编辑您的问题，并指出您遇到的困难，并可作为您问题的参考。

ll.poisson <- function(par, y) {
(alpha*beta*home) <- exp(par)
out <- sum(y * log((alpha*beta*home))) - length(y) * (alpha*beta*home)
return(out)
}

opt <- optim(par = 2, fn = ll.poisson, method = "BFGS", control = list(fnscale = -1),
y = mydata$HS)$par
mle <- exp(opt)

    HS  AS           Home             Away 
1    2  1           Arsenal       Aston Villa
2    1  2           Arsenal           Chelsea
3    0  0           Arsenal           Everton 
                .
                .
                .
378  2  2    Wigan Athletic Tottenham Hotspur 
379  1  2    Wigan Athletic   West Bromwich …  
380  2  1    Wigan Athletic   West Ham United