import pandas as pd aq=pd.read_csv('data/anscombesQuartet_hier.csv',header=[0,1],index_col=[0]) aq aq.mean() aq.var() tmp=pd.melt(aq) tmp['index']=list(range(int(len(tmp)/8)))*8 tmp[8:15] df=tmp.set_index(['group','index', 'var']).unstack() df.columns = [col[1].strip() for col in df.columns.values] df[7:15] df.reset_index(inplace=True) df[7:15] from ggplot import * ggplot(aes(x='x', y='y'), data=df)+facet_wrap('group',scales='fixed')+stat_smooth(method='lm',se=False) ggplot(aes(x='x', y='y'), data=df)+facet_wrap('group',scales='fixed')+stat_smooth(method='lm') ggplot(aes(x='x', y='y'), data=df)+geom_point()+facet_wrap('group')