Решено этим подходом
fig = plt.figure(figsize=[20,12])
fig, axs = plt.subplots(2, 2, sharex='col', sharey='row',gridspec_kw={'hspace': 0, 'wspace': 0})
#erase extra
fig.delaxes(axs[1][1])
#rb Li
axs[0,0].axis(xmin=autom['7Li'].min()-2*autom['er7Li'].max(), xmax=autom['7Li'].max()+2*autom['er7Li'].max())
axs[0,0].axis(ymin=autom['87Rb'].min()-2*autom['er87Rb'].max(), ymax=autom['87Rb'].max()+2*autom['er87Rb'].max())
for index in autom.index:
axs[0, 0].errorbar(autom.at[index,'7Li'],autom.at[index,'87Rb'], xerr=autom.at[index,'er7Li'], yerr=autom.at[index,'er87Rb'],
marker=autom.at[index,'symb'], mec='black', c=autom.at[index,'colorea'], label=autom.at[index,'nombre'])