coded Fig3 with ggplot2