In probability theory and statistics, the Hardy distribution is a discrete probability distribution that expresses the probability of the hole score for a given golf player. It is based on Hardy's (Hardy, 1945) basic assumption that there are three types of shots: 
 good , 
 bad  and 
 ordinary , 
where the probability of a good hit equals , the probability of a bad hit equals  and the probability of an ordinary hit equals . Hardy further assigned 

 a value of 2 to a good stroke, 
 a value of 0 to a bad stroke and 
 a value of 1 to a regular or ordinary stroke. 

Once the sum of the values is greater than or equal to the value of the par of the hole, the number of strokes in question is equal to the score achieved on that hole. A birdie on a par three could then have come about in three ways: ,  and , respectively, with probabilities ,  and .

Definitions

Probability mass function

A discrete random variable  is said to have a Hardy distribution, with parameters ,  and   if it has a probability mass function given by:

 if m is odd

and

 if m is even

with

and

where

  is the par of the hole ()

  is the golf hole score () if  is even

  is the golf hole score () if  is odd

  is the probability of a good shot ()

  is the probability of a bad shot () and ()

The moment generating function is given by:

  if m is odd

and

 if m is even

with

and

Each raw moment and each central moment can be easily determined with the moment generating function, but the formulas involved are too large to present here.

Hardy distribution for a par three, four and five 
For a par three:

For a par four:

Note the resemblance with . For a par five:

Note the resemblance with the formulas for  and .

History 
When trying to make a probability distribution in golf that describes the frequency distribution of the number of strokes on a hole, the simplest setup is to assume that there are only two types of strokes: 

 A good stroke with a probability of   
 A bad stroke with a probability of . 
 while
 a good shot then gets the value 1 and
 a bad shot gets the value 0. 

Once the sum of the shot values equals the par of the hole, that is the number of strokes needed for the hole. 
It is clear that with this setup, a birdie is not possible. After all, the smallest number of strokes one can get is the par of the hole. Hardy (1945) probably realized that too and then came up with the idea not to assume that there were just two types of strokes: good  and bad , but three types:

 good  with probability     
 bad  with probability     
 ordinary  with probability .   

In fact, Hardy called a good shot a supershot and a bad shot a subshot.  Minton later called Hardy's supershot an excellent shot  and Hardy's subshot a bad shot . In this article, Minton's excellent shot is called a good shot . Hardy came up with the idea of three types of shots in 1945, but the actual derivation of the probability distribution of the hole score was not given until 2012 by van der Ven.

Hardy assumed that the probability of a good stroke was equal to the probability of a bad stroke, namely . This was confirmed by Kang:

In retrospect, Hardy might well have been right, as the data in Table 2 in van der Ven (2013) show. This table shows the estimated - and -values for holes 1-18 for rounds 1 and 2 of the 2012 British Open Championship. The mean values were equal to 0.0633 and 0.0697, respectively. Later Cohen (2002) introduced the idea that  and  should be different. Kang says about this:

For the Hardy distribution the values of  and  may be different.

Goodness of fit 
The Hardy distribution gives the probability distribution of a single player's hole score. It takes several observations to perform a goodness-of-fit test (see Goodness of fit test) to check whether the Hardy distribution applies or not. This can be done with a single individual by having the individual play the same hole multiple times. Goodness-of-fit tests assume pure replications (see Replication (statistics)). This means that there should be no change in the player's golfing ability during repeated play of the hole. For example, there should not be an ongoing learning process (see Learning). Such effects cannot really be ruled out. One way around this problem is to use multiple players who can be assumed to have approximately the same golf proficiency. Such players are the participants in professional golf tournaments (see PGA Tour). Before using a goodness-of-fit test, it should first be checked that the participants indeed have approximately the same golf proficiency. This can be done separately for each hole by using, for example, the Pearson correlation coefficient between the hole score on the first day and the second day of a tournament. If there are no systematic differences (see Classical test theory) between players, the correlation (see Correlation) between the score achieved on Day 1 on a hole and the score achieved on Day 2 on that hole will not differ significantly (see Statistical significance) from zero. This can be easily tested statistically. In a study by van der Ven, the results of a goodness-of-fit test of the Hardy distribution were reported using the hole-by-hole scores from the 2012 Open Championship played at the St Andrews Golf Club. The distribution has been tested separately for each hole. Pearson's chi-squared test was used to determine whether the observed sample frequencies of the hole scores differed significantly from the expected frequencies according to the Hardy distribution. The fit between observed and expected frequencies was generally very satisfactory.

References 
Notes

Probability distributions