# Statistical significance: Is my new bot really better than the previous one?

A test statistic (generically named *T*) indirectly measures the unlikeliness that some observed result can be explained by random deviation from a simple reference case where nothing special happens (examples: The coin is fair, the dice is fair, the two bots are as good as each other).
In the use scope of this paper, the test statistic we’ll calculate is named χ21${\chi}_{1}^{2}$, however we’ll use a signed version of it.

### Use scope

- The formulas below are for comparing two bots only: Games with more than two bots (1 vs. 1 vs. 1 for instance) are not covered here.
- The formulas allows for draws, and you don’t need to know how the game behaves for a match between two identical bots.

### Formula of statistical significance

with:

*W*the number of matches won by the new bot,*D*the number of drawn matches,*L*the number of matches lost by the new bot,*N*=*W*+*D*+*L*the number of matches played.

## Proof

A test statistic compares actual observation with what we would expect if nothing special happens: That expected outcome is the null hypothesis *H*₀, and *T* = 0 if the actual observation exactly fits *H*₀.
As the outcomes to compare are triplets (values being labeled respectively “wins”, “draws” and “losses” for the new bot) for which we increment one of the values at each match result, we’re comparing trinomial distributions.
We’ll thus use Pearson’s *χ*² test of goodness of fit, which is an asymptotic test statistic, not exact but precise enough for our need (calculating an exact test statistic would cost a lot of computing resources for a slight and useless gain in precision).

In our case the null hypothesis is that the two bots are as good as each other (if not identical), so the expected distribution is an equal number of wins and losses; The expected number of draws is undefined because it depends on the game rules and on the bot (especially if it does not play perfectly), this is a free parameter for *H*₀.
So, for a given number *N* of matches played, the actual observation is a (*W*, *D*, *L*) triplet while the expected outcome is a (*α*, *N*−2*α*, *α*) triplet with *α* ∈ [0, *N*/2] (otherwise the triplet would contain negative values).
We have 2 degrees of freedom for the actual observation (as *N* = *W* + *D* + *L* is fixed, we can choose freely only two of the three values) and 1 free parameter for *H*₀ (namely *α*), so the remaining number of degrees of freedom for the test statistic is 2 − 1 = 1: Pearson’s *χ*² test statistic is thus χ21${\chi}_{1}^{2}$ in our case.

Whatever its degrees of freedom, Pearson’s *χ*² test statistic is:

with Oᵢ actual observations and Eᵢ expected values. So, in our case:

The free parameter *α* of *H*₀ is estimated by minimizing χ21${\chi}_{1}^{2}$ in the [0, *N*/2] domain.
Its first derivative is:

We have an extremum when the first derivative vanishes:

As we want *α*₀ ∈ [0, *N*/2], *N*−2*α*₀ and *α*₀ are both positive, so we use the positive root:

As we have taken the positive root, *N*−2*α*₀ and *α*₀ are either both positive or both negative, but that latter case is impossible as *N* > 0, so *α*₀ ∈ [0, *N*/2] as expected.
χ21${\chi}_{1}^{2}$ has a unique extremum in the [0, *N*/2] domain, at *α*₀; We then need to check that this extremum is a minimum thanks to the second derivative:

so

χ21${\chi}_{1}^{2}$ is thus minimal at *α*₀, and its value is:

A *χ*² test statistic tells only how far we are from *H*₀, without a notion of direction: Its value is the same by exchanging *W* and *L*.
As we have only 1 degree of freedom, we can introduce a direction by adding a sign to the test statistic (Note: It has no meaning for more than 1 d.f.).
As we want to evaluate the new bot against the old one, we want the test statistic to be positive if *W* > *L* and negative if *W* < *L* (given *H*₀, it is already 0 if *W* = *L*):

Note: You don’t need to bother about sign(0), as the right-hand multiplier’s value is 0 if *W* = *L*.
So you can simply implement the sign function as `sign(W, L) := W > L ? 1 : -1`

.

### Numerical interpretation

Below is a table giving the relation between the absolute value of *T* and the invert of the asymptotic one-tailed *p*-value, i.e., of the probability *p*₁ that mere luck could have given the observed result or something more extreme, for some values of *T*.
Graphs above show that *p*₁ is not about “better than” but about “more extreme than”, i.e., “farther from 0”, hence the flip when *W* < *L*.
Values for 1/*p*₁ are rounded to 2 significant digits only: The test statistic is asymptotic so the probability is not exact, and the precise probability is not important, its order of magnitude is enough for a good feeling of what it represents.

|T| | 1/p₁ |
---|---|

0 | 2 |

10 | 1300 |

20 | 260ᴇ3 |

30 | 46ᴇ6 |

40 | 7.9ᴇ9 |

50 | 1.3ᴇ12 |

60 | 210ᴇ12 |

70 | 34ᴇ15 |

80 | 5.3ᴇ18 |

90 | 840ᴇ18 |

100 | 130ᴇ21 |

## Proof

|*T*| is χ21${\chi}_{1}^{2}$, Pearson’s *χ*² test statistic with 1 degree of freedom.
The *χ*² distribution with 1 degree of freedom has PDF:

Note: The symbol “χ21${\chi}_{1}^{2}$” is shared by the test statistic value and the distribution function, even though they are two different mathematical objects: They must not be confused.

By definition of Pearson’s *χ*² test statistic with 1 d.f. the asymptotic two-tailed *p*-value *p*₂ for |*T*| verifies:

with “erf” the error function. So

with “erfc” the complementary error function.
As *p*₁ = *p*₂/2, we finally have 1/*p*₁ = 2/erfc(√(|*T*|/2)).

For instance, as the world population will soon reach 8 billion people, we can expect that, at any given time, someone on the planet is reaching a statistical significance of |*T*| ≈ 40 for some fair random process: Probably someone else, but possibly you while measuring your new bot performance (very unlikely, but still reasonably possible).

### Formula of win rate

There are several ways to define a win rate taking drawn matches into account, the simplest one is to count draws as half wins and half losses:

### Practical use

Use a run length limit of 5000 to 10,000 matches (depending on the time it would take) and a threshold for |*T*| around 50:

- no less than 40, as 1/
*p*₁ would be below the world population, - no more than 60, as the statistical significance is high enough.

Once the length limit or the threshold is reached, launch a new **independent** run in order to verify that the result is similar.
Note: in case of win rate being precisely 50 % (identical bots) or extremely close to that, *T* will vary randomly in the [−10, 10] range (with possible temporary pikes outside that range), so the first run may end with *T* ≈ 10 when the length limit is reached (after a lot of movements in the range) and the second run with *T* ≈ −10, they will nevertheless be “similar” in the sense that |*T*| does not increase visibly on average.