add R script to perform comparison plot