Create plot_hmr_pmd.R