Added first plot and associated R function