<p><div><div><div><div data-show-faces="false" data-height="21" data-width="100" data-layout="button_count" data-send="true" data-href="https://www.r-bloggers.com/quidditch-is-it-all-about-the-snitch/"></div></div><div></div></div><div></div></div>
<div>
(This article was first published on <strong><a href="http://mattansb.blogspot.com/2018/11/quidditch-is-it-all-about-snitch.html"> Mattan S. Ben-Shachar</a></strong>, and kindly contributed to <a rel="nofollow" href="https://www.r-bloggers.com/">R-bloggers)</a> </div>
<p>Much has been said about how the game of Quidditch is ruined by the scoring system – specifically how it makes no sense that the snitch is worth 150 points <i>and</i> that catching it ends the game [<a rel="nofollow" href="https://www.reddit.com/r/HPMOR/comments/2cqq7n/why_does_the_snitch_ruin_quidditch/" target="_blank">1</a>, <a rel="nofollow" href="https://www.youtube.com/watch?v=Zfg1I5f2HZY" target="_blank">2</a>, <a rel="nofollow" href="https://www.reddit.com/r/FanTheories/comments/2nk118/in_defense_of_quidditch/" target="_blank">3</a>]. Most of these arguments seem to revolve around the claim that it is nearly impossible to win a match of Quidditch without catching the snitch. Is this true? Let’s try and answer this question formally using statistics and R simulations:</p>
<blockquote><p>What is the <i>probability </i>of winning a Quidditch match without catching the snitch?</p></blockquote>
<table cellspacing="0" align="center" cellpadding="0">
<tbody>
<tr>
<td><a rel="nofollow" href="https://i2.wp.com/3.bp.blogspot.com/-vdR0-S17fuE/W_lfELTf0wI/AAAAAAAA3dA/EBjvzpDC11Ac4mmbuBDrkkmERribllBFgCLcBGAs/s1600/freegifmaker.me_2dxgP.gif?ssl=1" target="_blank"><img border="0" src="https://i1.wp.com/3.bp.blogspot.com/-vdR0-S17fuE/W_lfELTf0wI/AAAAAAAA3dA/EBjvzpDC11Ac4mmbuBDrkkmERribllBFgCLcBGAs/s640/freegifmaker.me_2dxgP.gif?resize=444%2C250&ssl=1" data-original-width="444" width="444" data-original-height="250" data-recalc-dims="1" height="250"></a></td>
</tr>
<tr>
<td><i>Victor Krum caught the snitch for Bulgaria, but Ireland still won (1994 World cup)</i></td>
</tr>
</tbody>
</table>
<p><a name="more"></a></p>
<h2> The Game, the Rules</h2>
<p>Quidditch is a fictional team sport played on broomsticks, devised by author J. K. Rowling, that is played by wizards in the fantasy world of Harry Potter. Like any sport, the game of Quidditch has many rules, but I will summarize here only the ones pertinent to our discussion:</p>
<ol>
<li>10 points are gained by throwing the Quaffle through one of the opponent team’s three hoops.</li>
<li>150 points are gained by catching the Golden Snitch, which can only be caught by the Seeker.</li>
<li>The game ends when the Golden Snitch is caught.</li>
<li>The winning team is the one has the most points at the end of the game.</li>
</ol>
<div>Given these rules and game objectives, in order to estimate the probability of winning a match we must first estimate (1) how long does it take a seeker to catch the Snitch on average, which determines how long a game lasts and thus how much time is allotted for the teams to score goals,  and (2) what is the rate of goal scoring.</div>
<h3>Estimating the Average Game Length</h3>
<div>Unfortunately, not many pro-Quidditch matches are described in the books, so this is quite hard to estimate. On the one hand, the only game that is described in its entirety is the 1994 Quidditch World Cup final, and it was extremely short: only 15 minutes or so [<a rel="nofollow" href="http://harrypotter.wikia.com/wiki/1994_Quidditch_World_Cup" target="_blank">4</a>]; On the other hand, the longest game ever lasted 6 months [<a rel="nofollow" href="https://www.quora.com/What-was-the-longest-Quidditch-match-ever-played-and-how-did-it-end" target="_blank">5</a>]. So the average is somewhere between 15 minutes, and 260,000 minutes. .. great.</div>
<div>For the sake of our simulations, then, I will assume the average game lasts 90 minutes (assuming Rowling has soccer or rugby in mind). I will also assume that the length of a Quidditch game is exponentially distributed, since the probability of catching the Snitch remains constant (as opposed to it getting easier as the game progresses) – the Snitch’s ability to escape capture does not change (magic) and players are substituted throughout the game, so let’s assume an unchanging level of vigilance.</div>
<div>Together, the distribution of game lengths looks something like this:</div>
<div>
<div><a rel="nofollow" href="https://i0.wp.com/1.bp.blogspot.com/-hTjWMZaW2K4/W_lQgXAZuXI/AAAAAAAA3cw/bLGsPGjKt6stP2ACLyoiDvLqL32qpjMRwCLcBGAs/s1600/fig1.png?ssl=1" target="_blank"><img border="0" src="https://i0.wp.com/1.bp.blogspot.com/-hTjWMZaW2K4/W_lQgXAZuXI/AAAAAAAA3cw/bLGsPGjKt6stP2ACLyoiDvLqL32qpjMRwCLcBGAs/s400/fig1.png?resize=450%2C1062&ssl=1" data-original-width="450" width="450" data-original-height="1062" data-recalc-dims="1"></a></div>
</div>
<div>This makes the 6 month long game <i>very</i> unlikely ( $Pr(time>260,000)&lt;10^{-10,000,000}$ ), but the 15 minute game rather probable ( $Pr(time&lt;15)=0.15$ ).</div>
<h3>Estimating the Rate of Goals Scored</h3>
<div>Throughout the books, the final scores of several Quidditch games are mentioned, but since we don’t know their lengths, it is hard to estimate the rate of goals-per-minute (GPM). Going by the 15 minute world cup, in which the final score was 170-160, or 170-10 if only counting goal points, we can see that there were 18 goals in 15 minutes, or an average of 9 goals per 15 minutes per team, or a GPM of 0.6 (or 0.07 for Bulgaria, and 1.13 for Ireland). Though this is only one game, and a very short one at that, these numbers seem close to the averages seen in basketball, where GPM are around 40/48 = 0.83 [<a rel="nofollow" href="https://www.teamrankings.com/nba/stat/field-goals-made-per-game" target="_blank">6</a>]. This is also in line with the game’s depiction in the films as a fast sport, with many goals scored quickly and in succession.</div>
<div>I will also assume that this rate is constant and unchanging throughout the game, regardless of the score (or fatigue of the players, see above), and so will assume GPM has a poisson distribution, with an average of 0.8, making the GPM distribution look something like this:</div>
<div>
<div><a rel="nofollow" href="https://i1.wp.com/4.bp.blogspot.com/-7-dbMQTNP_0/W_lQkxyAhJI/AAAAAAAA3c0/03ZuEzs2NR8e-6_ZZMEaoCYPekvcOOBiQCLcBGAs/s1600/fig2.png?ssl=1" target="_blank"><img border="0" src="https://i1.wp.com/4.bp.blogspot.com/-7-dbMQTNP_0/W_lQkxyAhJI/AAAAAAAA3c0/03ZuEzs2NR8e-6_ZZMEaoCYPekvcOOBiQCLcBGAs/s400/fig2.png?resize=450%2C1062&ssl=1" data-original-width="450" width="450" data-original-height="1062" data-recalc-dims="1"></a></div>
</div>
<h2>Simulating the Games</h2>
<div>(If you’re not <i>that</i> into reading raw code, you can just skip to the conclusions at the end.)</div>
<h3>R Setup</h3>
<div></div>
<h3>Equal Teams</h3>
<div>Let’s start by simulating 1 million games between two teams with equal GPMs and equally good seekers.
<div></div>
<p>Under these conditions, the probability of winning a game without catching the snitch is:</p>
<div></div>
<p>3.75%! Hmm, that does seem quite unfair… Though it does put Fred and George’s gambling odds at around 1:20!</p>
<p>How many games were determined by the snitch – i.e., had the snitch been caught at the same time by the other team, they would have won?</p>
<div></div>
<p>94.9%! Not the best of odds…</div>
<h3>Unequal (but matched) Teams</h3>
<div>But just like in real sports, different teams have different strengths… Let’s see what happens when we create two teams who have an equal probability of winning, but differ in terms of their seekers and GPMs.</div>
<div></div>
<div>After some trial and error, I found that if we have a team with a seeker who catches the snitch within 67 minutes on average and a GPM of 0.71, versus a team with a seeker who catches the snitch in within 113 minutes on average but with a GPM of 0.90, both teams have a 50% chance of winning.
<p> So, what is the chance of winning <i>without catching the snitch</i>?</p>
</div>
<div></div>
<div>12.2%! Still pretty low… very disappointing. Let’s explore each team separately:
</div>
<div></div>
<div>It seems that 99.9% of team A’s wins were due to their better seeker…</div>
<div>…but that only 24.3% of team B’s wins were achieved by their weaker seeker <img alt="™" src="https://s.w.org/images/core/emoji/2.4/72x72/2122.png" class="wp-smiley">!</div>
<div></div>
<div>How many games were determined by the snitch?
</div>
<div></div>
<div>90.1%!</div>
<h2>Conclusion</h2>
<div>
<ul>
<li>If you’ve caught the snitch, you’ve probably won.</li>
<li>If you’ve won, you probably caught the snitch.</li>
<li>If you’ve lost, catching the snitch would have probably changed that.</li>
</ul>
</div>
<div>But… the majority of points are non-snitch related. And it seems to me unlikely that you would have a seeker that was twice as good as the opposing team’s, and so compensating for a weaker seeker seems very likely given this simulation. SO it’s probably best to invest in a better chaser who can bring up your GPM and a better keeper who can lower your opponents’ GPM, than in a seeker who can shave 5 minutes off of his snitch catching time…</div>
<div></div>
<div>Also, maybe someone should change the rules? Perhaps the snitch should <a rel="nofollow" href="https://en.wikipedia.org/wiki/Quidditch_(sport)" target="_blank">only be worth 30 points</a>?</div>
	<p>
<div>
<div>To <strong>leave a comment</strong> for the author, please follow the link and comment on their blog: <strong><a href="http://mattansb.blogspot.com/2018/11/quidditch-is-it-all-about-snitch.html"> Mattan S. Ben-Shachar</a></strong>.</div>
<hr>
<a rel="nofollow" href="https://www.r-bloggers.com/">R-bloggers.com</a> offers <strong><a rel="nofollow" href="https://feedburner.google.com/fb/a/mailverify?uri=RBloggers">daily e-mail updates</a></strong> about <a rel="nofollow" href="https://www.r-project.org/" title="The R Project for Statistical Computing">R</a> news and <a rel="nofollow" href="https://www.r-bloggers.com/search/tutorial" title="R tutorials">tutorials</a> on topics such as: <a rel="nofollow" href="https://www.r-bloggers.com/search/data%20science" title="Data science">Data science</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/Big%20Data" title="Big Data">Big Data, <a rel="nofollow" href="https://www.r-users.com/" title="R jobs">R jobs</a>, visualization (<a rel="nofollow" href="https://www.r-bloggers.com/search/ggplot2" title="ggplot and ggplot2 tutorials">ggplot2</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/boxplot" title="Boxplots using lattice and ggplot2 tutorials">Boxplots</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/map" title="Maps and gis">maps</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/animation" title="Animation in R">animation</a>), programming (<a rel="nofollow" href="https://www.r-bloggers.com/search/RStudio" title="RStudio IDE for R">RStudio</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/sweave" title="Sweave and literate programming">Sweave</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/LaTeX" title="LaTeX in R">LaTeX</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/SQL" title="SQL and databases">SQL</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/eclipse" title="Eclipse IDE for R">Eclipse</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/git" title="git and github, Version Control System">git</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/hadoop" title="Large data in R using Hadoop">hadoop</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/Web+Scraping" title="Web Scraping of google, facebook, yahoo, twitter and more using R">Web Scraping</a>) statistics (<a rel="nofollow" href="https://www.r-bloggers.com/search/regression" title="Regressions and ANOVA analysis tutorials">regression</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/PCA" title="principal component analysis tutorial">PCA</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/time+series" title="Time series">time series</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/trading" title="finance trading">trading</a>) and more...
<img width="1" alt="" src="http://feeds.feedburner.com/~r/RBloggers/~4/WheTnbvMmec" height="1"></a></div>