# How to model football (soccer) scores

I'm trying to model football scores using the poisson distribution. I'm calculating a lambda for each team and then multiplying the probabilities of two poisson calculations to get the probability of a specific score.

e.g. poison(lambdaHome, 1)*poisson(lambdaAway, 2) to calculate the probability of a 1-2 result at the end of the match.

The problem I'm having is that the poisson distribution assumes that the second goal is just as likely as the third goal is just as likely as a fourth goal. But in reality this shouldn't be true. If a match has a score of 2-0 then intuitively it is less likely for the team to score another goal is it would probably spend more time on defence to keep the score than it would do attacking to score more goals.

Is there a probability distribution similar to poisson, but rather than having a constant probability for all goals, makes each goal less and less likely?

-