我想根据某些行索引在成对图中着色点。这里是我用来绘制一个变量与另一个变量的代码。成对图中的着色点
cases<-which(rownames(data_no_na) %in% colnames(tumor_data))
controls<-which(rownames(data_no_na) %in% colnames(control_data))
plot(y=range(pca[,1]),x=range(pca[,2]),type='n',xlab="Principle Component 2",ylab="Principle Component 1", main="Iterative Thresholding Sparse PCA")
points(y=pca[cases,1], x=pca[cases,2], col = 'red')
points(y=pca[controls,1], x=pca[controls,2], col = 'blue');
一个简单的对情节是这样的:
pairs(pca[,1:3])
编辑:例:
cases<-1:10
controls<-11:20
pca<-matrix(c(rnorm(3*10,0,1),rnorm(3*10,5,1)),nrow=20,ncol=3)
你应该给一个重复的例子。 – agstudy 2013-03-24 14:37:36
这意味着我们需要查看您的真实数据或可以说明上述情节的一些虚拟数据的输出。如果你的数据不包含许多列,试着粘贴'dput(head(cases))'&'dput(head(controls))'的输出。 – 2013-03-24 14:47:27
你好!请通过看看[**如何使一个伟大的可重复的例子**],使您的文章可重现(http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-例如)为我们提供帮助。谢谢。 – Arun 2013-03-24 14:58:29