2014-05-18 52 views
1

我可以通过一个名为“日”,这从1-10变化的变量进行分区大型数据矩阵一个传奇,我想用这个上色我的PCA图着色在PCA曲线

该作品很好,当我使用下面的代码:

plot(pca$x,col = day,pch=20) 

但是当我尝试添加一个传奇表现哪天每种颜色是指:

legend("topleft", pch = 20, col=day, legend = c("Day 1","Day 2","Day 3","Day 4","Day 5","Day 6","Day 7","Day 8","Day 9","Day 10"), bty='n', cex=.75) 

我得到这个。如何让传说的情节??显示各个颜色: enter image description here

+1

这可能是因为'day'包含尽可能多的颜色数据点,但传说你只需要十位。也许'唯一(一天)'就够了,但这取决于你的数据是如何组织的。 –

回答

2

您可以使用

legend("topleft", col=1:10, legend = paste("Day", 1:10), 
    pch = 20, bty='n', cex=.75) 

要清除了一点东西。图例只应指定您在图中使用的不同颜色。

1

对于mtcars这会工作:

with(mtcars, plot(mpg, wt, col=cyl)) 
with(mtcars, legend("topright", pch = 20, col=unique(cyl), legend = unique(cyl))) 

也许你可以用你的数据做了同样的方式。

或者,如果你使用GGPLOT2:

qplot(wt, mpg, data=mtcars, colour=factor(cyl))