I've got a set of football results and I want to make a probabilty model of football scores as described in Dixon, Coles (1997, http://www.math.ku.dk/~rolf/teaching/thesis/DixonColes.pdf). They estimate the parameters based on maximum likelihood and the model assumes the variables are independent Poisson.
Now I have 2 questions:
- How should I transform my data to input them in the model?
- What packages in R are best for maximum likelihood estimation of this nature?