Displaying 1 result from an estimated 1 matches for "rat_pca".
2010 Apr 18
2
plotting pca of samples in different colors
...,132,147,163,166,215,275,drop=F)])
>#specifying columns to be colored in brown
> d=colnames(rat[,-c(1,5,9,12,16,154,227,276,284,298,305,311,318,336,338,2,3,4,6,13,14,7,8,10,11,15,29,35,36,50,54,57,59,61,63,70,71,79,99,110,132,147,163,166,215,275,drop=F)])
>#calculating pca for columns
> rat_pca=prcomp(t(rat),scale=T)
>plot(rat_pca$x,pch=20,col=c("red","blue","green","brown")[c(a1,a2,a3,d)])
But after executing last command there is nothing plotted in the graph
and it also does not give any error. Please help!
regards
Amit