Re: [R] draw circle on PCA plot

2009-11-26 Thread Wai Kuan Yip
Phoebe, Hi, this was how I solved that problem using plot(). I found using symbols to be hard to adjust... #plot the 1st vs 2nd component plot(x,y,ylab="2nd component", xlab="1st component", xlim=range(x), ylim=range(y), main="1st vs 2nd component analysis") #show the arrows, length defined by

[R] draw circle on PCA plot

2009-11-25 Thread phoebe kong
Hi, I have a hard time in drawing circle on PCA. I have 60 samples. The corresponding PC1 scores and PC2 scores were stored as "mergedata". Here are the summary of PCs scores. > range(mergedata[,"PC1"]) [1] 0.0085 0.0100 > range(mergedata[,"PC2"]) [1] 0.0032 0.0075 > mean(mergedata[,"PC1"]) [1] 0