我想基于某些行索引在一对图中着色点.这是我用于绘制1个变量与另一个变量的代码.
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' );
Run Code Online (Sandbox Code Playgroud)
一个简单的对图是这样的:
pairs(pca[,1:3])
Run Code Online (Sandbox Code Playgroud)
编辑:示例:
cases<-1:10
controls<-11:20
pca<-matrix(c(rnorm(3*10,0,1),rnorm(3*10,5,1)),nrow=20,ncol=3)
Run Code Online (Sandbox Code Playgroud)
Rol*_*and 21
像这样的东西?
cols <- character(nrow(iris))
cols[] <- "black"
cols[iris$Species %in% c("setosa","versicolor")] <- "blue"
cols[iris$Species == "virginica"] <- "red"
pairs(iris,col=cols)
Run Code Online (Sandbox Code Playgroud)