add DM457overDM456 plotting R code