On 11/06/2009 03:21 AM, Joel Fürstenberg-Hägg wrote:
Hi all,
I'm making a PCA plot with eight variables (columns taken from a larger data frame
"fieldTrial0809[idx,c(39:46)]"). I then want the symbols in the plot to be colored as a
gradient from red to blue, depending on the value of another column in "fieldTrial0809[idx,
c(48)]" containg temperatures from -12.1 to -5.4.
I don't want to use the heat.colors(n) function, because then I just get colors
from red to yellow...
I've created a vector "colTemperature" with 12 colors and I've also made a new vector
"tempACC", based on the values of "fieldTrial0809[idx, c(48)]" with 12 levels.
What I don't know, is how to incorporate this vector in the text() function
below:
df=fieldTrial0809[idx,c(39:46)]
PC=prcomp (df[,1:8])
hcl=hclust(dist(df))
plot(PC$x[,1],PC$x[,2], col="white", main="8 variables")
text(PC$x[,1],PC$x[,2], labels=cutree(hcl,8), col=cutree(hcl,8)) # Change
this, but how???
library(plotrix)
text(...,col=color.scale(fieldTrial0809[idx,48],
extremes=c("red","blue"))
Finally, I've got another question, how do you decide how many components to use in
"text(... labels=cutree(hcl,8), ...)"? Do you try with different values until
i looks nice or what..?
As in so many other areas of life, beauty is in the eye of the beholder.
Jim
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.