# The Mathematics of Duelling

Duelling with pistols. If you were the one issuing the challenge, your dilemma was that custom dictated that your adversary be allowed to shoot first. Only then, if you were still able to shoot, would you be permitted to seek “satisfaction”.

How much of an advantage does the first shooter really have? In this article, we build a simple probability model, and implement a numerical model in a few lines of R code.

Two gentleman face off in the snow. Convention dictates the challenged shoots first.

The Mathematics of Duelling & First Shooter Advantage

Duelling with pistols has a fascinating history. The dilemma facing a would-be challenger was that custom dictated that his adversary be allowed to shoot first. Once the challenged had fired his shot, then, assuming the challenger was still alive and able to shoot, he would be allowed to seek “satisfaction” by firing off a returning shot. Back and forth the volleys would go, turn by turn, though it was considered civilised to declare all satisfaction achieved and full honours restored if both parties remained unharmed after three rounds.

The game theoretic aspects of the first-shooter convention are interesting. For the one perceiving that a slight has been given, the convention presents an obvious deterrent. How likely are you to be seriously wounded before you have get the chance to win back your honour? Or in other words, how much of an advantage does the first shooter really have? For a professional agents provocateur, 1 the opposite question was key. How much does the skill of a marksman influence first-shooter advantage? Indeed, how long could a marksman-provocateur expect to live in such a career?

To model this problem, I’ll use a simple probability model to calculate the “first-shooter advantage” (FSA). It’s also a good opportunity to illustrate how one can numerically model such problems using just a few lines of code in R, an Open-Source Statistical Computing platform. This allows varying the assumptions very quickly and recalculating.

Assumptions
To model the problem, we’ll assume:

• that as per tradition [3], both shooters fire alternating shots
• that this continues back and forth indefinitely until “satisfaction” is obtained by one party.
• for simplicity, that both are equally talented shooters, i.e. they share the same single shot probability of success.
• that emotions (anger, fear, or ebulliance at having survived a round) don’t influence the single shot probabilities.

What do we find?

It turns out that the First Shooter Advantage (FSA) is not as high as one might expect. Even if both shooters have a 7 in 10 chance of success on any given shot, victory for the first shooter is by no means assured, standing at a rather uncertain 50-50 odds.

An interesting fact emerges. The most widely accepted duelling code, The Code Duello (1777, Ireland, [1]), appears to deliberately reduce overall mortality by requiring the use of “smooth-bore barrels as opposed to rifled barrels that cause the bullet to spin and give it greater accuracy and range”. How Duels Work, [4].

With the marksmanship of combatants deliberately kept low, first shooter advantage also decreases. Indeed, assuming two shooters share a low 1 in 4 (25%) likelihood of a strike, the first shooter’s advantage is now less than 15%, and declines even faster as shooting accuracy falls further.

There is a theory that the heightened conventions of politeness common in gentlemanly speech were reinforced by the dueling culture.
The Art of Manliness — Duelling [2]. But looking at the evidence, we can also say that the convention to decrease the accuracy of dueling weapons, by reducing first shooter advantage, kept the dueling culture viable, as there was less risk to being a challenger.

So let’s turn now to the technical details.

The Mathematics of Duelling & First Shooter Advantage

How might one go about calculating first shooter advantage?

(Warning: Spoiler ahead. Stop reading now if you want to work out the problem yourself first.)

Solution Outline

Each shooter engages in a sequence of Bernoulli trials that stops after a first success by either shooter. So for each shooter, we obtain a geometric distribution which gives the probability of a fixed number of trials to first success.

Then Shooter A wins in k shots if his first success is the kth shot and if Shooter B has missed on all k-1 previous shots. Shooter B wins on the kth shot if B is successful on shot k and if Shooter A has missed on all k previous shots. Notice the asymmetry here.

Victory on the kth shot is therefore:
$\displaystyle V_A(k) = (1-p_A)^{k-1}(p_A)(1-p_B)^{k-1}$

$\displaystyle V_B(k) = (1-p_B)^{k-1}(p_B)(1-pA)^{k}$

where $p_A, p_B$ are the individual shot sucess probabilities of Shooter A and B respectively.

If we assume that the duel continues to the death then k is unknown and unbounded, so each shooter’s victory probability must sum over the natural numbers. If instead, we respect the less severe maximum three rounds convention, then we can restrict the sum to the first three terms of the series.

$\displaystyle V_A = \sum_{k=1}^{\infty} (1-p_A)^{k-1}(p_A)(1-p_B)^{k-1}$ $\displaystyle V_B = \sum_{k=1}^{\infty} (1-p_B)^{k-1}(p_B)(1-p_A)^{k}$

What is the First Shooter Advantage? It is the difference between the victory probabilities of shooters A and B, i.e.

$\displaystyle FSA = V_A - V_B$

The asymmetry arises from the fact that shots are alternated, starting with Shooter A, and because of the possibility of immediate death.

Numerical Model

A little code in R calculates the results quickly, and allows exploration of possibilities by changing the parameters of the model.

First calculate the victory probabilities for A or B winning on shot k, given their individual shot success probabilities pA and pB respectively, as follows:


prAfirst <- function(k,pA,pB) { (1-pA)^(k-1) * (pA) * (1-pB)^(k-1) }
prBfirst <- function(k,pA,pB) { (1-pB)^(k-1) * (pB) * (1-pA)^(k) }


Then calculate a partial sum that, in the limit as N gets arbitrarily large, converges to the First Shooter Advantage (FSA) that shooter A enjoys over shooter B even when both shooters are equally skilled with success probability p=pA=pB.


FSA <- function(p,N) {
sum(prAfirst(seq(1,N),p,p)) - sum(prBfirst(seq(1,100),p,p))
}


Now use sapply, one of R’s powerful functional programming paradigms, to iterate FSA over each element of an array of success probabilities. Each invocation of FSA calculates the first N=100 terms of the partial series, which is adequate in this case to reach acceptable convergence. The results are gathered into a table fsa_tab, showing for each shot probability p, the advantage that shooter A has over shooter B.


options(digits=3)
fsa_tab <- cbind(seq(0.1,1,0.1), sapply(seq(0.1,1,0.1), FSA, 100))
colnames(fsa_tab) <- c('p', 'FSA')
fsa_tab   # show table


First Shooter Advantage given both shooters have shot success probability p

Finally, the results can be plotted. Observe that for a crack marksman approaching 100% success probability p, survival as an agents provocateur would be more or less assured. But given the strict convention of accuracy suppressing equipment (smooth-bore barrels instead of spiral-grooved barrels), success probabilities even for an excellent marksman are likely to have been significantly less than 100%. In this case, the plot shows that first shooter advantage comes down geometrically (roughly quadratically) with decreasing shot-probability.


plot( seq(0.1,1,0.1), fsa_tab, type="b",
main="The Mathematics of Duelling",
xlab="Probability p of Shot Success",
ylab="First Shooter Advantage (A over B)"
)


Plot showing First Shooter Advantage

Extensions

Introduce into the model above the convention of maximum three rounds to a duel.

Problem How likely is it that both contestants come away unharmed?

I leave this as an exercise to the reader. Feel free to reply in the comments.

Take-aways

So, does this change how we view duelling with pistols? Or the decisions we would make were we in the last few decades of the 19th century?

Apart from generally avoiding conflict, you’d certainly want to hold your peace if your antagonist was a known first-class marksman. And you’d definitely want your second to check that his weapon was of the less accurate kind specified in the Duelling Code.

How likely would it be that you were singled out by an agent provocateur? Most mathematicians know the story of the brilliant young mathematician Evariste Galois who was cut down before his prime in a duel with pistols amidst some suspicious circumstances. These accounts of Galois’ death attribute the cause of the duel to agent provocateurs. More recent scholarship by Tony Rothman has shown most convincingly that such a set-up could not have been further from the truth.

Marksmen aside, in all other cases the first shooter advantage declines quite rapidly with decreasing accuracy. The convention of using less accurate duelling pistols meant that for your average gentleman the custom of deciding controversies by arms was as much dependent on luck as skill.

Ah, the ways of honour!

[5] Genius and Biographers: The Fictionalization of Evariste Galois, American Mathematical Monthly, 89, 84, 1982

[7] All of Statistics: A Concise Course in Statistical Inference, by Larry Wasserman, Springer, 2004

There are many excellent references for using R. I have found the following by Hadley Wickham to be one of the best in flattening out the learning curve in going from basic to advanced R programming

[9] Advanced R by Hadley Wickham, Chapman & Hall, forthcoming June 2014

Picture Credits
The picture of two men facing off in the snow is a painting by Ilya Repin, entitled Eugene Onegin and Vladimir Lensky’s duel. Courtesy of Wikipedia.

Footnotes

1. An agent provocateur, literally ‘inciting agent’, is one hired to give sufficient offence to provoke a challenge but without raising suspicions of a set-up.