Playing Poker after taking Viagra. Analysing the effect from the math perspective

Introduction

If you watched a movie called Limitless or read a book called 'The Dark Fields' by Alan Glynn (Limitless is an adaptation), you probably contemplated the idea of taking a pill that would make you smarter. If you did not watch the movie or read the book, then I can describe it briefly: the main character is a loser (an underdog) who, after taking a pill, becomes extremely smart and changes his life significantly. Even though Viagra makes a different impact on men, there is still an interesting idea I want to test. And the book was just an inspiration. I want to understand whether conditions that people try to avoid, like fear, adrenaline, pain, unsatisfied sexual desire, etc., can actually be a way to channel our energy into better decision-making and improve cognitive performance.

Experiment Design

The experiment is going to be an AB test, where Group A (control) - is me playing without Viagra. Group B - me playing after taking Viagra.

Structure of the experiment:

Mathematical Structure of the Experiment:

So, I needed to play over 3k hands in each group, that would be over 6k hands total. And I considered it was time to learn more about Viagra - the dosage, side effects, durability, price, etc. Here is the list of the main things that I learned:

I also found interesting facts about Viagra while doing the research. Here are some of them with links:

So, it was time for me to calculate how many Viagra pills I will need. It was time to do some math again. I decided to play 3 tables for 4 hours. In 6-max NLHE, the number of hands you play per hour is about 90 (based on my statistics from HM3, I had 94 with unknown standard deviation; I did not want to overcomplicate anything). I concluded I would play 90 * 3 = 270 hands per hour. 4 hours of play would give me 4 * 270 = 1080 hands per day. So, that would take me 3 days of playing in each group, and 150mg of Viagra to finish the experiment.

The first problem was that the drug was not available. It was surprising considering how expensive that is, and the economy is in a recession. I did not buy anything online because I found an article that about 50% of the online Viagra sold is fake and dangerous. Was that the risk I was willing to take? And what do they mean "dangerous"? Pollution-dangerous? Sterility-dangerous? Death-dangerous? Not worth it. I decided to drive around the city and found in one pharmacy 2 packs of Viagra, 100mg in each. So, I had to buy both. In the end, I had 200mg of Viagra.

Playing under Viagra

Playing without Viagra was not anything special for me. The same ups and downs, the same players, the same default settings. I remember the times when it used to be so scary to bluff, and the insane amounts of adrenaline I got from bluffing. It is in our human nature to get used to everything that happens around us too often. For example, if you are perpetually exposed to heights, high-speed driving, or talking to girls, you get less and less rush every consecutive time. I noticed that I was still not very confident when I had to bluff. This was just something I had to do to play a balanced game.

The Viagra experience was, of course, diabolical. 30 minutes after I took it for the first time, I already felt my heart beating. And I could feel a pulse on my dick even without touching it. Honestly, the first hour was the most difficult to stay focused. I was thinking about too many things, like I thought I was going to die. But then I managed to channel my extensive energy into the game itself.

Do you remember I mentioned that Viagra requires sexual stimulation to work? It seems like it does not have to be sexual, because when I overcame the feeling that I would pass away with a hard-on and focused on the game, I noticed the difference. I wanted to make more risky decisions. Not risky in terms of just an intuition, but based on statistical factors. I was more aggressive with tight players, trying to find an optimal balance of getting as much as I could from that person by losing as little as possible. Of course, I also did it without Viagra, but I was still more conservative. It was like my mind was clear, like I was looking into the future with so much anticipation. Placebo? Maybe. Maybe the confidence thing. But not from the biological standpoint. Remember, I mentioned that an erect penis makes a man more attractive and increases confidence? You often need that confidence in poker.

Groups Analysis

Before doing the t-test, I had to make sure the groups happened to be balanced. In other words, I needed to know how many pairs I had in group A and group B, how many premium hands I had in group A and B, the share of hands by position in each group, what positions I played the most (though I was confident the positions would be balanced), etc.

Overall, I played 3247 hands in Group A (without Viagra), and 3338 hands in Group B (with Viagra). As you see: 3247!=33383247 != 3338

We must mathematically prove that we can work with these group volumes. I would not like to incorporate other techniques to fix that. We will use χ2\chi^2 test to test our null hypothesis:

H0: the number of hands in Group A == the number of hands in Group B

Here is the math formula:

χ2=i=1rj=1c(OijEij)2Eij\chi^2 = \sum_{i=1}^{r} \sum_{j=1}^{c} \frac{(O_{ij} - E_{ij})^2}{E_{ij}} \qquad

but I will use the Python scipy library to calculate the p-value:

from scipy.stats import chisquare

counts = df["group"].value_counts()
expected = [len(df) / 2, len(df) / 2]
chi2, p_value = chisquare(counts.values, f_exp=expected)

I got p-value:

np.float64(0.26211439424422495)

Thankfully our p>α(0.05)p > \alpha(0.05). It means we do not reject our H0H0 hypothesis, that states:

H0: the number of hands in Group A == the number of hands in Group B

Let us move to internal analysis.

Here are the results:

Overall Premium Table

GroupCategoryHandsPercent
AMedium662.03%
APremium421.29%
AStrong260.80%
AWeak311395.87%
BMedium682.04%
BPremium401.20%
BStrong220.66%
BWeak320896.11%

Overall Category Percentages

GroupMediumPremiumStrongWeak
A2.03%1.29%0.80%95.87%
B2.04%1.20%0.66%96.11%

Hands per Group × Position

GroupPositionHands
ABB542
ABTN541
ACO541
AHJ541
ASB541
AUTG541
BBB556
BBTN557
BCO556
BHJ556
BSB557
BUTG556

Category by Position (first 12 rows)

GroupPositionCategoryCountTotalPercent
ABBMedium195423.51%
ABBPremium65421.11%
ABBStrong35420.55%
ABBWeak51454294.83%
ABTNMedium65411.11%
ABTNPremium75411.29%
ABTNStrong35410.55%
ABTNWeak52554197.04%
ACOMedium105411.85%
ACOPremium95411.66%
ACOStrong55410.92%
ACOWeak51754195.56%

Pivoted Table — Category Percentages by Group & Position

GroupPositionMediumPremiumStrongWeak
ABB3.51%1.11%0.55%94.83%
BTN1.11%1.29%0.55%97.04%
CO1.85%1.66%0.92%95.56%
HJ1.48%1.29%1.11%96.12%
SB1.48%1.66%1.11%95.75%
UTG2.77%0.74%0.55%95.93%
BBB1.98%1.26%0.54%96.22%
BTN2.33%1.44%0.90%95.33%
CO2.52%1.08%1.08%95.32%
HJ1.98%1.80%0.18%96.04%
SB1.97%0.36%0.18%97.49%
UTG1.44%1.26%1.08%96.22%

And visualization:

Even though some of the combinations, depending on how you group them, are significantly different, in order to get very accurate and balanced distribution you need to play about 100.000 hands. The more hands you play, the closer the results are to their true values. It is our fundamental Law of Large Numbers [2]. You probably understand what problems I would face trying to get to 100.000. Erectile dysfunction for life, maybe even death.

I did not perform an analysis of the games themselves because I initially expected the games to be different. Because of the Viagra effect. So, I did not statistically analyze how many flops or rivers I played in each group because an imporvement or deterioration of my strategy must be clear after performing an AB test. Preflop hands, though, can show whether I had an upstrick/downstrick in a certain group, which I did not.

Win Rate by Group

When I looked at my win rate, I was not very happy, not gonna lie. I set up the test to catch about 5 bb/100 difference, but we got:

Group A - 5.6 BB/100

Group B = 7.1 BB/100

It is only 1.5 bb/100 difference. In order to catch that type of difference and confidently say that Viagra improves your poker skills, I had to play about 28k hands in each damn group. But we must let the math do the talking here before concluding anything.

Based on the results, let us do the t-test:

SE=σ12n1+σ22n2SE = \sqrt{ \frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2} }

t=xˉ2xˉ1SEt = \frac{\bar{x}_2 - \bar{x}_1}{SE}

df=(σ12n1+σ22n2)2σ14n12(n11)+σ24n22(n21)df = \frac{ \left( \frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2} \right)^2 } { \frac{\sigma_1^4}{n_1^2 (n_1 - 1)} + \frac{\sigma_2^4}{n_2^2 (n_2 - 1)} }

p=2(1Ft(t,df))p = 2 \cdot (1 - F_t(|t|, df))

sp=(nA1)sA2+(nB1)sB2nA+nB2s_p = \sqrt{ \frac{(n_A - 1)s_A^2 + (n_B - 1)s_B^2}{n_A + n_B - 2} }

import numpy as np
from scipy.stats import t

mean1 = 5.6     # bb/100
mean2 = 7.1     # bb/100
n1 = a / 100 # this is the number of blocks
n2 = b / 100

# Even though the number of hands is not that big I decided to calculate sigmas for each group
sigma1 = 86.4   # this is 86.4 bb/100 deviation, which is quite normal for NLHE
sigma2 = 91.3   # this is with Viagra

# calculating standard error
se = np.sqrt(sigma1**2 / n1 + sigma2**2 / n2)

# t-statistics
t_stat = (mean2 - mean1) / se

# degrees of freedom calculated by Welch–Satterthwaite equation
# I consider sigmas not to be equal, othervise we could calcualte it like 'n1 + n2 - 2'
# I decided not to do another test on sigmas
df_num = (sigma1**2 / n1 + sigma2**2 / n2)**2
df_den = (sigma1**4 / (n1**2 * (n1 - 1))) + (sigma2**4 / (n2**2 * (n2 - 1)))
df = df_num / df_den

p_value = 2 * (1 - t.cdf(abs(t_stat), df))

print(f"t-stat = {t_stat:.4f}")
print(f"df = {df:.2f}")
print(f"p-value = {p_value:.4f}")
t-stat = 0.0685
df = 63.80
p-value = 0.9456

So, the p-value is 0.9456. And 0.9456>0.050.9456 > 0.05, so we have to conclude that the result is not statistically significant (p-value=0.95). We cannot reject our null hypothesis, which was:

H0 = Viagra does not improve poker skills

BUT, because the observed effect was much smaller (1.5bb/100) than the expected 5bb/100, 3000 hands is too small a number to detect that small difference. I want you to understand that it does not mean Viagra does not improve poker skills (or decision-making); it means it might improve the skills, but probably by less than 5bb/100, which we were unable to detect. Or I am just coping.

Conclusion

Yes, I failed to prove such an important concept in science. But I swear I could feel the effects. The confidence boost and desire to make riskier (but calculated) decisions. I truly believe it increased my stress tolerance. I can still feel it. Which means the experiment was not a waste of my time and resources. I felt it expanded my comfort zone, which, in theory, must open more opportunities for me. The only thing left for me to do is to find them. And I still have 50mg of Viagra to seize at least one of those opportunities.