Joyplot

This example is taken from the “Common Plots II/ Ridge, or ‘joy’, plots” sub-section of the book Coding for Economists (1st edition).

import pandas as pd

from lets_plot import *

LetsPlot.setup_html()

df = pd.read_pickle("https://github.com/aeturrell/coding-for-economists/raw/main/data/berkeley_data.pkl")
df.head()

	Date Number	Year	Month	Day	Day of Year	Anomaly
0	1880.001	1880	1	1	1	-0.786
1	1880.004	1880	1	2	2	-0.695
2	1880.007	1880	1	3	3	-0.783
3	1880.01	1880	1	4	4	-0.725
4	1880.012	1880	1	5	5	-0.802

final_year = df["Year"].max()
first_year = df["Year"].min()

breaks = [y for y in list(df.Year.unique()) if y % 10 == 0]
(
    ggplot(df, aes("Anomaly", "Year", fill="Year"))
    + geom_area_ridges(scale=20, alpha=1, size=0.2, trim=True, show_legend=False)
    + scale_y_continuous(breaks=breaks, trans="reverse", format='d')
    + scale_fill_viridis(option="inferno")
    + ggtitle(
        "Global daily temperature anomaly {0}-{1} \n(°C above 1951-80 average)".format(
            first_year, final_year
        )
    )
)