code: def draw_jointplot(figname, x, y, data=None, kind="reg", color=None,
                     xlim=None, ylim=None, format="pdf"):
      
      import seaborn as sns
      sns.set_context('talk')
      plt.clf()
  
      register = {"MeanCoverage": "Sample Mean Coverage",
                  "HD.FDP": "Depth of full spanning reads",
                  "HD.PDP": "Depth of partial spanning reads",
                  "HD.PEDP": "Depth of paired-end reads",
                  "HD.2": "Repeat size of the longer allele"}
  
      g = sns.jointplot(x, y, data=data, kind=kind, color=color,
                        xlim=xlim, ylim=ylim)
      g.ax_joint.set_xlabel(register.get(x, x))
      g.ax_joint.set_ylabel(register.get(y, y))
      savefig(figname + "." + format, cleanup=False)
docs: Wraps around sns.jointplot