One option is to use the y data from the graphs - probably the most useful for emissions (pilots)
_, bp = pd.DataFrame.boxplot(df, return_type='both') outliers = [flier.get_ydata() for flier in bp["fliers"]] boxes = [box.get_ydata() for box in bp["boxes"]] medians = [median.get_ydata() for median in bp["medians"]] whiskers = [whiskers.get_ydata() for whiskers in bp["whiskers"]]
But it might be easier to get other values ββ(including IQR) with
quantiles = df.quantile([0.01, 0.25, 0.5, 0.75, 0.99])
or as suggested by WoodChopper
stats = df.describe()
source share