我制作了一个PCA图,其中我根据各种基因的表达绘制了许多细胞。在这个图中,我想用单独的颜色对一些点进行着色。我试图通过创建“组”来实现这一点,我根据它们的表达或缺乏“gene1”的表达来分类细胞。R:使用自动绘图时基于PCA的PCA中的颜色点
这里就是我的数据帧的外表(基因1,基因2和cell_1,cell_2等都是colnames和rownames):
gene1 gene2 gene3 gene4 gene5
cell_1 0.0000 0.279204 25.995400 46.171700 94.234100
cell_2 0.0000 23.456000 77.339800 194.241000 301.234000
cell_3 2.0000 13.100000 45.309200 0.776565 0.000000
cell_4 0.0000 10.500000 107.508000 3.032500 0.000000
cell_5 3.0000 0.000000 0.266139 0.762981 123.371000
下面是我用它来尝试实现这一目标的代码:
library(ggplot2)
library(ggfortify)
# Group cells based on expression of a certain gene (to use for color labels in the next step)
groups <- factor(ifelse(df$gene1 > 0, "Positive", "Others"))
#Calculate PCs and plot PCA
autoplot(prcomp(log(df[]+1)), colour="Positive")
当我运行此代码时,出现以下错误:
Error in grDevices::col2rgb(colour, TRUE) : invalid color name 'Positive'
您应该使用组变量,而不是一个特定的标签,从组变量。 – Pieter
@Pieter谢谢。我应该提到 - 我尝试过与“团体”和“其他人”也有同样的错误。 (错误grDevices :: col2rgb(color,TRUE):无效的颜色名称'groups') – Galaffer