Рассчитать эксцесс для цикла для каждого именованного столбца в фрейме данных и графике

Question

Я хочу рассчитать эксцесс для каждого столбца моего фрейма данных и построить его с помощью цикла for, мои столбцы имеют имена a1, a2, a3, a4. после этого я хочу нанести их с заголовками "a1", "a2" и т. д. В настоящее время я делаю каждый столбец отдельно, но это не короткий путь. Как я могу зациклить это для кадра данных размером 4 столбца?

from scipy.stats import kurtosis
import matplotlib.pyplot as plt
import numpy as np 
import pandas as pd

temp=df['a1']
y=temp.kurtosis(axis=0, level=0)
plt.title('a1')
plt.show()

temp=df['a2']
y=temp.kurtosis(axis=0, level=0)
plt.title('a2')
plt.show()