'Tie Fighter' Errorbar Charts in Python With Plotnine and Altair
It’s December 2020 and it’s time to check out different python plotting packages. I primarily will use pandas dataframe plotting with either seaborn’s
sns.set() or setting plotly as the plotting backend during data exploration, but there’s nothing in python that matches R’s
ggplot2. I’ve known about plotnine as the python equivalent but last time I tried it it was still lacking in the feature department. I’m happy to say that it’s totally viable now.
One chart that I often used in R was a ’tie fighter’ (error bar) chart. It’s handy to compare accuracy for various predictive models. Here’s a motivating example from ’thecodeforest’ where he compares various time serious forecasting models in R. I’m also motivated to stick with ggplot because then I’d be able to quickly make plots in R, python, and in Julia (given how seamless
RCall is in Julia and native support of dataframes).
Here’s what the syntax for a quick tie fighter chart in
plotnine using the included
mtcars dataset as a toy example:
from plotnine import * from plotnine.data import mtcars import pandas as pd import numpy as np from scipy.stats import sem # Calc mean, std error of mpg by number of cylinders and whether they're automatic (0) or manual (1) df = mtcars.groupby(['cyl', 'am'])['mpg'].agg([np.mean, sem]).reset_index() p = (ggplot(df, aes('factor(cyl)', 'mean', color='factor(am)')) + geom_point() + geom_errorbar(aes(ymin = 'mean - sem', ymax = 'mean + sem')) + coord_flip() + labs(title="A Tie Fighter Chart", x = 'Cylinders', y = 'Mean MPG +/-1 SE') ) p
And what I love about
plotnine) is the consistent syntax and how easy it is to add or remove customization and complexity compared to the other plotting libraries in python:
# Now it's trivial to add or remove stuff to our plots, even in python! p + theme(xkcd)
One thing we’re still missing is the ability to do subtitles and captions, but maybe someday.
Altair is another opinionated plotting package based on Vega-lite that has a very consistent syntax that’s worth a look.
Here’s the same example as above, but in Altair:
import altair as alt from plotnine.data import mtcars points = alt.Chart(mtcars).mark_point().encode( x='average(mpg):Q', y='cyl:O', color='am:N', ).properties(height=200) error_bars = alt.Chart(mtcars).mark_errorbar(extent='stderr', ticks=True).encode( x=alt.X('average(mpg):Q', scale=alt.Scale(zero=False)), y=alt.Y('cyl:O'), color='am:N', ) points+error_bars
Altair has baked in aggregation and transformations so I didn’t have to do my
pandas. And similar to what I like in
plotnine there’s the ability to quickly add and remove complexity. But as you can see the result was so small on my macbook pro, and I wish there was an easy way to scale the whole thing up to a reasonable aspect ratio (looks like I’d have to scale up each component separately or save as SVG).