We're assuming zero-vig so the expectation on each bet is zero.
The variance of a bet of size x at decimal odds of d is given by (d-1)*x
2. So for example the variance of a 1-unit bet at +130 would be 1
2 * (-1) = 1.3, while the variance of a 2-unit bet at -105 would be 2
2 * (-1) ≈ 3.8095.
If we ignore compounding the variance of a linear combination of independent bets is the sum of the individual variances. In the provided data set the total variance works out to be about 272.86 units
2. The standard deviation is the square root of the variance, which in this case is (272.86 units
2)
½ ≈ 16.519 units.
The realized return is
20/
16.510 ≈ +1.2108 standard deviations better than expected, implying a p-value of =NORMSDIST(1.2108) ≈ 88.701% from the normal distribution and =1-TDIST(1.2108,93-1,1) ≈ 88.545% from the t-distribution with 92 (i.e., 93 bets - 1) degrees of freedom (where NORMSDIST and TDIST refer to MS Excel's cumulative standard normal and Student t-distribution functions, respectively).
One problem with the above analysis, however, is that it completely discounts the effect of compounding. Nevertheless, because we're dealing with a relatively small number of bets, each at a relatively small percentage of bankroll, this effect is fairly small.
Taking compounding into account, the only real procedural difference would be in the calculation of variance. Because the result of previous bets impact the results of future bets (insofar as if we win an earlier bet we'd be betting more on a later bet it's no longer correct that the total variance equals the sum of the variances.
Rather than go into some long derivation, I'll just state the result here. The total variance is given by the sum of the variances plus the sum of the product of the variances taken 2 at a time plus the sum of the product of the variances taken 3 at a time ... plus the product of all N variances.
This obviously represents a huge number of terms (2
N - 1 terms for N bets, 9,903,520,314,283,042,199,192,993,792 terms for 93 bets), but if bets are sufficiently small we can approximate by only looking at terms up to only the second or third order. Given the 93 bets of the original problem, the total variance (up to the second order) would then be:
Total Variance ≈ ΣNi=1[σ2i] + ΣNi=1[σ2i * ΣNj=i+1[σ2j ]]
which works out to a standard deviation of 16.629%, implying a t-distribution (with 92 degrees of freedom) p-value of 88.391%. Were we to go to the
third order we'd of course find a slightly higher σ (~16.630%) and hence a slightly lower p-value (~88.390%).
This is
faily close to the p-value of 88.215% that I obtained from my 40,000,000-trial Monte Carlo run. The difference stems from having approximated the distribution of results with the t-distribution.