about · email me · subscribe
Spurious correlation #4,761 · View random

A linear line chart with years as the X-axis and two variables on the Y-axis. The first variable is Krispy Kreme Doughnuts store count in the US and the second variable is Total likes of SmarterEveryDay YouTube videos.  The chart goes from 2007 to 2021, and the two variables track closely in value over that time. Small Image
Download png
, svg

AI explanation

As the aroma of freshly glazed doughnuts wafts through the air, viewers find themselves drawn to the nearest computer or mobile device, where they can simultaneously indulge in both delicious content and delectable treats. It's a well-known fact that the gravitational pull of sugary confections is directly linked to an uptick in intellectual curiosity and appreciation for science. The combination of sugar rushes and mind-bending experiments is the ultimate recipe for boosting SmarterEveryDay's online popularity. So, whether it's the mesmerizing sight of doughnuts being made or the subconscious connection between circular shapes and the laws of physics, one thing is certain – when doughnuts delight, educational videos ignite!

Model: dalle-3
Prompt: Generate an image of a delicious spread of freshly baked Krispy Kreme doughnuts on a table, with a large screen in the background displaying captivating science-filled videos from the SmarterEveryDay YouTube channel. The foreground should be filled with people interacting with the doughnuts, while the background shows the increasing number of Krispy Kreme stores in the US, symbolizing the rising popularity of both Krispy Kreme and SmarterEveryDay. Capture the blend of indulgence and curiosity as viewers are drawn to both the mouthwatering treats and the intriguing educational content. The depiction should include only scenes, people, and objects, it should not include text or speech bubbles.
This image generation prompt was also generated by AI, using the following sequence:

Model: gpt-3.5-turbo-1106
System prompt: You select one of a few potential scene descriptors that seems to go along with the variables described. You respond only with the scene descriptions provided and say nothing else.
User prompt: Here are two variables: Krispy Kreme Doughnuts store count in the US and Total likes of SmarterEveryDay YouTube videos. Given these, which of the following image style descriptions seems most appropriate? Pick the one that seems to be related to one of the variables at hand. Please respond with only one of the following potential descriptions: Steampunk a funny image macro Cybernetic art An image that looks like it is a scene from 'Star Wars' Minimalism
ChatGPT: a funny image macro

Model: gpt-3.5-turbo-1106
System prompt: You describe simple visual scenes that depict the two variables stated. You combine both variables into one scene.
User prompt: Please describe a visual scene that prominently depicts both Krispy Kreme Doughnuts store count in the US and Total likes of SmarterEveryDay YouTube videos, in the following situation: As the aroma of freshly glazed doughnuts wafts through the air, viewers find themselves drawn to the nearest computer or mobile device, where they can simultaneously indulge in both delicious content and delectable treats. It's a well-known fact that the gravitational pull of sugary confections is directly linked to an uptick in intellectual curiosity and appreciation for science. The combination of sugar rushes and mind-bending experiments is the ultimate recipe for boosting SmarterEveryDay's online popularity. So, whether it's the mesmerizing sight of doughnuts being made or the subconscious connection between circular shapes and the laws of physics, one thing is certain – when doughnuts delight, educational videos ignite!. Keep the description simple.
ChatGPT: Amidst the enticing scent of freshly baked Krispy Kreme doughnuts, a screen showcases the mesmerizing science-filled videos of SmarterEveryDay. In the background, the rising number of Krispy Kreme stores in the US fuels the surge in likes for SmarterEveryDay's YouTube channel, as the two elements blend into a delightful mix of sugary indulgence and educational curiosity.

Model: gpt-3.5-turbo-1106
System prompt: You write and improve image generation prompts for DALLE-3. You maintain the same subjects, but enhance the visuals.
User prompt: Please respond with an image generation prompt that would be likely to create an interesting visual using DALLE-3, using the following style: [ a funny image macro ] and depicting this scene: Amidst the enticing scent of freshly baked Krispy Kreme doughnuts, a screen showcases the mesmerizing science-filled videos of SmarterEveryDay. In the background, the rising number of Krispy Kreme stores in the US fuels the surge in likes for SmarterEveryDay's YouTube channel, as the two elements blend into a delightful mix of sugary indulgence and educational curiosity. Keep the visual depiction simple. Do not include any words or speech bubbles.
ChatGPT: Generate an image of a delicious spread of freshly baked Krispy Kreme doughnuts on a table, with a large screen in the background displaying captivating science-filled videos from the SmarterEveryDay YouTube channel. The foreground should be filled with people interacting with the doughnuts, while the background shows the increasing number of Krispy Kreme stores in the US, symbolizing the rising popularity of both Krispy Kreme and SmarterEveryDay. Capture the blend of indulgence and curiosity as viewers are drawn to both the mouthwatering treats and the intriguing educational content.

*Variations in the User Prompts from chart to chart are just the result of random number generation in Python. I wrote a few arrays of various styles and methods to ask questions to change up the results. Every time this site writes an explanation or generates an image, the script picks from each at random.

I sequence the requests into multiple prompts because I find GPT 3.5 to perform much better with short, well-managed contexts. Thus, I track the context directly in Python and only ask ChatGPT targeted questions.

Model: gpt-3.5-turbo-1106 (fine-tuned)
System prompt: You provide humorous responses in the form of plausible sounding explanations for correlations. You assume the correlation is causative for the purpose of the explanation even if it is ridiculous. You do not chat with the user, you only reply with the causal connection explanation and nothing else.
User prompt: Please make up a funny explanation for how an increase in Krispy Kreme Doughnuts store count in the US caused Total likes of SmarterEveryDay YouTube videos to increase.


Random correlation

Discover a new correlation

View all correlations

View all research papers

Report an error


Data details

Krispy Kreme Doughnuts store count in the US
Source: Statista
See what else correlates with Krispy Kreme Doughnuts store count in the US

Total likes of SmarterEveryDay YouTube videos
Detailed data title: Total likes of SmarterEveryDay YouTube videos.
Source: SmarterEveryDay on YouTube
See what else correlates with Total likes of SmarterEveryDay YouTube videos

Correlation r = 0.5114681 (Pearson correlation coefficient)
Correlation is a measure of how much the variables move together. If it is 0.99, when one goes up the other goes up. If it is 0.02, the connection is very weak or non-existent. If it is -0.99, then when one goes up the other goes down. If it is 1.00, you probably messed up your correlation function.

r2 = 0.2615996 (Coefficient of determination)
This means 26.2% of the change in the one variable (i.e., Total likes of SmarterEveryDay YouTube videos) is predictable based on the change in the other (i.e., Krispy Kreme Doughnuts store count in the US) over the 15 years from 2007 through 2021.

p < 0.10(Null hypothesis significance test)
The p-value is 0.051. 0.0513230935498341850000000000
The p-value is a measure of how probable it is that we would randomly find a result this extreme. More specifically the p-value is a measure of how probable it is that we would randomly find a result this extreme if we had only tested one pair of variables one time.

But I am a p-villain. I absolutely did not test only one pair of variables one time. I correlated hundreds of millions of pairs of variables. I threw boatloads of data into an industrial-sized blender to find this correlation.

Who is going to stop me? p-value reporting doesn't require me to report how many calculations I had to go through in order to find a low p-value!
On average, you will find a correaltion as strong as 0.51 in 5.1% of random cases. Said differently, if you correlated 19 random variables Which I absolutely did.
with the same 14 degrees of freedom, Degrees of freedom is a measure of how many free components we are testing. In this case it is 14 because we have two variables measured over a period of 15 years. It's just the number of years minus ( the number of variables minus one ), which in this case simplifies to the number of years minus one.
you would randomly expect to find a correlation as strong as this one.

[ -0, 0.81 ] 95% correlation confidence interval (using the Fisher z-transformation)
The confidence interval is an estimate the range of the value of the correlation coefficient, using the correlation itself as an input. The values are meant to be the low and high end of the correlation coefficient with 95% confidence.

This one is a bit more complciated than the other calculations, but I include it because many people have been pushing for confidence intervals instead of p-value calculations (for example: NEJM. However, if you are dredging data, you can reliably find yourself in the 5%. That's my goal!


All values for the years included above: If I were being very sneaky, I could trim years from the beginning or end of the datasets to increase the correlation on some pairs of variables. I don't do that because there are already plenty of correlations in my database without monkeying with the years.

Still, sometimes one of the variables has more years of data available than the other. This page only shows the overlapping years. To see all the years, click on "See what else correlates with..." link above.
200720082009201020112012201320142015201620172018201920202021
Krispy Kreme Doughnuts store count in the US (Number of Stores)2512252222292302392530294307316355364365373
Total likes of SmarterEveryDay YouTube videos (Total likes)853565779502742614327570691345135320018129002132170174874018425202571370312148036405502646970




Why this works

  1. Data dredging: I have 25,237 variables in my database. I compare all these variables against each other to find ones that randomly match up. That's 636,906,169 correlation calculations! This is called “data dredging.” Instead of starting with a hypothesis and testing it, I instead abused the data to see what correlations shake out. It’s a dangerous way to go about analysis, because any sufficiently large dataset will yield strong correlations completely at random.
  2. Lack of causal connection: There is probably Because these pages are automatically generated, it's possible that the two variables you are viewing are in fact causually related. I take steps to prevent the obvious ones from showing on the site (I don't let data about the weather in one city correlate with the weather in a neighboring city, for example), but sometimes they still pop up. If they are related, cool! You found a loophole.
    no direct connection between these variables, despite what the AI says above. This is exacerbated by the fact that I used "Years" as the base variable. Lots of things happen in a year that are not related to each other! Most studies would use something like "one person" in stead of "one year" to be the "thing" studied.
  3. Observations not independent: For many variables, sequential years are not independent of each other. If a population of people is continuously doing something every day, there is no reason to think they would suddenly change how they are doing that thing on January 1. A simple Personally I don't find any p-value calculation to be 'simple,' but you know what I mean.
    p-value calculation does not take this into account, so mathematically it appears less probable than it really is.




Try it yourself

You can calculate the values on this page on your own! Try running the Python code to see the calculation results. Step 1: Download and install Python on your computer.

Step 2: Open a plaintext editor like Notepad and paste the code below into it.

Step 3: Save the file as "calculate_correlation.py" in a place you will remember, like your desktop. Copy the file location to your clipboard. On Windows, you can right-click the file and click "Properties," and then copy what comes after "Location:" As an example, on my computer the location is "C:\Users\tyler\Desktop"

Step 4: Open a command line window. For example, by pressing start and typing "cmd" and them pressing enter.

Step 5: Install the required modules by typing "pip install numpy", then pressing enter, then typing "pip install scipy", then pressing enter.

Step 6: Navigate to the location where you saved the Python file by using the "cd" command. For example, I would type "cd C:\Users\tyler\Desktop" and push enter.

Step 7: Run the Python script by typing "python calculate_correlation.py"

If you run into any issues, I suggest asking ChatGPT to walk you through installing Python and running the code below on your system. Try this question:

"Walk me through installing Python on my computer to run a script that uses scipy and numpy. Go step-by-step and ask me to confirm before moving on. Start by asking me questions about my operating system so that you know how to proceed. Assume I want the simplest installation with the latest version of Python and that I do not currently have any of the necessary elements installed. Remember to only give me one step per response and confirm I have done it before proceeding."


# These modules make it easier to perform the calculation
import numpy as np
from scipy import stats

# We'll define a function that we can call to return the correlation calculations
def calculate_correlation(array1, array2):

    # Calculate Pearson correlation coefficient and p-value
    correlation, p_value = stats.pearsonr(array1, array2)

    # Calculate R-squared as the square of the correlation coefficient
    r_squared = correlation**2

    return correlation, r_squared, p_value

# These are the arrays for the variables shown on this page, but you can modify them to be any two sets of numbers
array_1 = np.array([251,225,222,229,230,239,253,0,294,307,316,355,364,365,373,])
array_2 = np.array([8535,65779,5027,42614,327570,691345,1353200,1812900,2132170,1748740,1842520,2571370,3121480,3640550,2646970,])
array_1_name = "Krispy Kreme Doughnuts store count in the US"
array_2_name = "Total likes of SmarterEveryDay YouTube videos"

# Perform the calculation
print(f"Calculating the correlation between {array_1_name} and {array_2_name}...")
correlation, r_squared, p_value = calculate_correlation(array_1, array_2)

# Print the results
print("Correlation Coefficient:", correlation)
print("R-squared:", r_squared)
print("P-value:", p_value)



Reuseable content

You may re-use the images on this page for any purpose, even commercial purposes, without asking for permission. The only requirement is that you attribute Tyler Vigen. Attribution can take many different forms. If you leave the "tylervigen.com" link in the image, that satisfies it just fine. If you remove it and move it to a footnote, that's fine too. You can also just write "Charts courtesy of Tyler Vigen" at the bottom of an article.

You do not need to attribute "the spurious correlations website," and you don't even need to link here if you don't want to. I don't gain anything from pageviews. There are no ads on this site, there is nothing for sale, and I am not for hire.

For the record, I am just one person. Tyler Vigen, he/him/his. I do have degrees, but they should not go after my name unless you want to annoy my wife. If that is your goal, then go ahead and cite me as "Tyler Vigen, A.A. A.A.S. B.A. J.D." Otherwise it is just "Tyler Vigen."

When spoken, my last name is pronounced "vegan," like I don't eat meat.

Full license details.
For more on re-use permissions, or to get a signed release form, see tylervigen.com/permission.

Download images for these variables:


View another random correlation

How fun was this correlation?

Kudos for rating!


Correlation ID: 4761 · Black Variable ID: 416 · Red Variable ID: 25506
about · subscribe · emailme@tylervigen.com · twitter

CC BY 4.0