Editing Probability (section)

===Complex combinations/What are the odds of a specific number of things happening?===

Having saved the rogue, we now go to 40k. We have a squad of ten tactical marines firing bolters into 3 plague zombies that we really need cleared off a victory point. For those that stumbled in here without knowing 40k (welcome to /tg/!), each marine needs to roll several d6s to take out a zombie. First, they need 3+ on a d6 to hit, then 3+ on a d6 to wound, then the zombie needs to fail a 5+ Feel No Pain roll. What are the chances those tac marines will wipe out the zombies? Should we point a valuable squad of assault marines at them for support?

This is going to have an easy part and a hard one. The easy part is the chance of ''one'' marine killing ''one'' zombie. A d6 has a 4/6 chance of getting a 3+, and also a 4/6 chance of missing a 5+. This marine, then, has a 4/6 * 4/6 * 4/6 = 64/216 = 8/27 = '''29% chance''' to clip a zombie from the group. We've got ten marines, so our 'expected' number of dead zombies is 10*.29, or 2.9, just a little under the three we want. That suggests we should probably call up the assault marines, as there's a good chance we won't clear that whole group out.

Now, the tricky part. This is a really close game, and those assault marines are desperately needed elsewhere. We want the ''exact'' possibility the plague zombies won't survive. Open up [http://www.wolframalpha.com Wolfram Alpha], and buckle up. '''Skip to the TL;DR if you just want a highly useful formula for mathhammering.'''

We have ten marines and we're looking to put at least three unsaved wounds downrange, so we're looking for the number of ways to score precisely 3 successes in 10 tests. To make matters a little worse, we actually also need the number of ways to score 4 successes, or 5, or 6, or more, up to ten. Inflicting more wounds is not a problem for us. To work this out, we'll need the probability for a given scenario of hits, and the number of ways that scenario can occur.

Let's look at inflicting exactly one wound. This requires one marine to succeed and all the others to fail. If we said that marine were the first guy in the squad we could tell how likely that is by multiplying everything together. We know the first marine succeeds and '''all''' the others fail, so this is [[File:Probability Eq6.png|x30px]] one success and nine fails, in that order. P(failure) is just 1-P(success), so we've got .29 and .71.

The probability of any scenario like this, then, is just [[File:Probability Eq7.png|x30px]], or in this case, <nowiki>.29*.71*.71*.71*.71*.71*.71*.71*.71*.71 = .29*.71^9 = </nowiki>  .29*.0458 = '''.0133'''. Wow, it sure got tedious in here. So, this 1.3% chance is the chance of scoring exactly one kill in exactly this way.

We know there ten scenarios exactly like this one, differing only in which marine got the kill. We don't care about the difference, so those ten scenarios with 1.33% each give us 10*.0133 = .133, or a '''13.3% chance to kill exactly one zombie'''.

Now we can compute the possibility of any given scenario, but to inflict two wounds we would need to know how many ways you can organize two successes out of ten tests. Same thing for three, or four, or whatever else. For one success in ten tests, that was easy. For two or more, it's really annoying. Fortunately, we have this function on our calculators fittingly called a 'Combination':

Simply put, a Combination is the number of ways an event can occur in a set, ignoring order. This is important to us because we don't care what order these wounds are inflicted in and we don't want to spend all day multiplying numbers (only most of the day). A Combination combines two numbers, the number of tests and the number of successes.

I'm going to call the number of tests T and the number of successes S. The Combination of T and S is written 'T choose S', or TcS.

For the curious, the internal math on a Combination:
<div class="toccolours mw-collapsible mw-collapsed" style="width:800px">
Know what a factorial is? In math, the exclamation point pulls double duty both as a point of excitement and also as notation for 'multiply this number by everything below it'. 5!, for example, is 5*4*3*2*1. Everyone learns about them in middle school, then nobody ever uses them.  ...until now.

[[File:Probability Eq7a.png|x90px]]

So, 10C3 = 10! / (3! * 7!) = (10*9*8*7*6*5*4*3*2*1)/((3*2*1)*(7*6*5*4*3*2*1)). Cool thing is, the sequence 7*6*5*4*3*2*1 will always show up in the numerator and the denominator.  We can cross it right off.

10C3 = (10*9*8)/(3*2*1). Way better. Even by hand it's not too hard to get to 120. However, we didn't get to be the dominant species on the planet by doing things the hard way. Just make Wolfram Alpha do it: [https://www.wolframalpha.com/input/?i=combination(10,3)]

</div>

For the less curious, 10c3 = Make WolframAlpha do it: [https://www.wolframalpha.com/input/?i=combination(10,3) 120].


Alright, we've got our tool. Back on track. Ten marines, three zombies. The probability of any scenario that delivers '''S''' wounds from '''T''' marines with a per-marine kill rate of '''P''' is [[File:Probability Eq8.png|x30px]]. There are '''TcS''' ways for that scenario to occur. We're down to the good stuff:

Probability of inflicting exactly S wounds = [[File:Probability Eq9.png|x30px]].

Exactly three wounds, ten marines, three zombies? 10c3 * (.29^3)*(.71^(10-3)). [https://www.wolframalpha.com/input/?i=combination(10,3)*(.29%5E3)*(.71%5E(10-3)) Wolfram Alpha, away!] '''26.6%'''

Ever so close. That's the chance of exactly three unsaved wounds. We would be fine with four, or five, or any other number above three. So we need to either [https://www.wolframalpha.com/input/?i=sum+combination(10,j)*(.29%5Ej)*(.71%5E(10-j)),+j%3D3+to+10 add up everything above S=3], or [https://www.wolframalpha.com/input/?i=1-sum+combination(10,j)*(.29%5Ej)*(.71%5E(10-j)),+j%3D0+to+2 add the results for zero, one, and two wounds and subtract that from one]. '''59% chance.''' This game is going to be close, but the marines are more likely than not to win.

'''TL;DR''' The general form of all this is that if you have a number of tries T each with a success rate of P, the chance of getting exactly S successes is [[File:Probability Eq9.png|x30px]]. If you're looking for at least S successes, add the results of this formula for everything between S and T. Or, if you're a real human with things to do, [https://www.wolframalpha.com/input/?i=sum+(combination(T,j)*(P%5Ej)*((1-P)%5E(T-j))),+j%3DS+to+T+where+P%3D.29,+S%3D3,+T%3D10 make a machine do it].

This is pretty easy to modify. Just change the inputs a little bit to try things out. [https://www.wolframalpha.com/input/?i=sum+(combination(T,j)*(P%5Ej)*((1-P)%5E(T-j))),+j%3DS+to+T+where+P%3D(2%2F3*2%2F3*2%2F3),+S%3D3,+T%3D20 Marines are rapid firing]? [https://www.wolframalpha.com/input/?i=sum+(combination(T,j)*(P%5Ej)*((1-P)%5E(T-j))),+j%3DS+to+T+where+P%3D(2%2F3*2%2F3),+S%3D3,+T%3D10 Plague zombies are normal cultists]? [https://www.wolframalpha.com/input/?i=sum+(combination(T,j)*(P%5Ej)*((1-P)%5E(T-j))),+j%3DS+to+T+where+P%3D(2%2F3*5%2F6*2%2F3),+S%3D4,+T%3D8 Eight devastators firing S8 AP3 krak missiles at a T6 3+ Carnifex with 5+ cover]? Enjoy.