我正在研究一个名为expression
的基因表达数据框。我的样本属于不同的子组,并在名称中指出(即,所有在其名称中包含“ adk”的样本都属于同一子组)
adk1 adk2 bas1 bas2 bas3 non1 ...
gene1 1.1 1.3 2.2 2.3 2.8 1.6
gene2 2.5 2.3 4.1 4.6 4.2 1.9
gene3 1.6 1.8 0.5 0.4 0.9 2.2
...
我已经定义了子集使用
adk <- expression[grepl('adk', names(expression))]
然后我使用以下方法对此数据集进行了PCA>
pca = prcomp (t(expression), center = F, scale= F)
我现在想在PCA双线图中绘制从PCA获得的PC相互绘制的图。为此,我希望所有属于同一子组的样本具有相同的颜色(例如,所有“ adk”样本应为绿色,所有“ bas”样本应为红色,所有“ non”样本应为蓝色)。我尝试使用ggfortify中
color
函数的autoplot
参数,但我无法使其使用我定义的子集。
如果有人可以帮助我,我会很高兴!谢谢:)
我正在研究一种称为表达的基因表达数据框架。我的样本属于不同的子组,在名称中表示(即,所有在其名称中包含“ adk”的样本都属于...