R code for drawing PCA plots as in Belmonte et al.,2013