1

我正在使用 ggpubr 库中的 ggscatter 函数制作散点图。我的数据框看起来像这样

1   a   b   chr17   +   0.003   0.005   0,2 282232  4,0 253259  non_sig
10  a   b   chr22   -   0.733   0.6855  16,17   3,3 24,45   11,4    non_sig
12  a   b   chr13   +   0.7625  0.7965  22,14   1,7 7,18    1,4 non_sig
14  a   b   chr13   +   0.4555  0.369   20,16   19,12   4,23    17,11   non_sig
15  a   b   chr13   +   0.488   0.384   27,15   19,12   7,18    17,11   non_sig
16  a   b   chr16   -   0.9715  0.978   200141  3,2 260280  3,3 non_sig
21  a   b   chr1    +   0.9365  0.933   149118  1,12    133175  11,5    non_sig
22  a   b   chrX    +   0.6475  0.7265  129,57  58,35   104,78  37,29   non_sig
26  a   b   chr3    +   0.05    0.0475  54,32   721503  46,27   519617  non_sig
27  a   b   chr3    +   0.0475  0.045   57,34   721503  47,30   519617  non_sig

这是我正在使用的命令

library("ggpubr")
df <- read.table("test.txt",header =F,sep="\t")
ggscatter(df,x= "V6",y= "V7",color = "V12", shape = 21, size = 1,add = "reg.line",cor.coef = TRUE, cor.method = "pearson",conf.int = TRUE,title="A3SS(4561)",xlab="Ψ2",ylab = "Ψ1",
                  palette = c("black", "red"))

我想使用第12列为点non_sig着色signon_sigblacksigred

当我使用上面的代码时，它会做我想要的，但是我如何才能在这里专门为

sig=>red
non_sig=>black

谢谢您的帮助！！

4

2 回答 2

5

我假设“使用 12the column 为点着色”是指使用基于 column 的颜色填充点V12。

请注意，您的示例数据仅包含V12 = "non_sig"条目，因此我手动将一个条目更改为"sig"

library(ggpubr)
ggscatter(
    df,
    x= "V6", y= "V7",
    fill = "V12",
    shape = 21,
    size = 5,
    add = "reg.line",
    cor.coef = TRUE,
    cor.method = "pearson",
    conf.int = TRUE,
    title="A3SS(4561)",
    xlab="Ψ2",
    ylab = "Ψ1",
    palette = c("black", "red"))

在此处输入图像描述

样本数据

df <- read.table(text =
    "1   a   b   chr17   +   0.003   0.005   0,2 282232  4,0 253259  non_sig
10  a   b   chr22   -   0.733   0.6855  16,17   3,3 24,45   11,4    non_sig
12  a   b   chr13   +   0.7625  0.7965  22,14   1,7 7,18    1,4 non_sig
14  a   b   chr13   +   0.4555  0.369   20,16   19,12   4,23    17,11   non_sig
15  a   b   chr13   +   0.488   0.384   27,15   19,12   7,18    17,11   sig
16  a   b   chr16   -   0.9715  0.978   200141  3,2 260280  3,3 non_sig
21  a   b   chr1    +   0.9365  0.933   149118  1,12    133175  11,5    non_sig
22  a   b   chrX    +   0.6475  0.7265  129,57  58,35   104,78  37,29   non_sig
26  a   b   chr3    +   0.05    0.0475  54,32   721503  46,27   519617  non_sig
27  a   b   chr3    +   0.0475  0.045   57,34   721503  47,30   519617  non_sig", header = F)

更新

作为对您的评论的回应，您可以为您的论点使用命名向量；palette例如

df <- read.table(text =
    "1   a   b   chr17   +   0.003   0.005   0,2 282232  4,0 253259  non_sig
10  a   b   chr22   -   0.733   0.6855  16,17   3,3 24,45   11,4    non_sig
12  a   b   chr13   +   0.7625  0.7965  22,14   1,7 7,18    1,4 non_sig
14  a   b   chr13   +   0.4555  0.369   20,16   19,12   4,23    17,11   non_sig
15  a   b   chr13   +   0.488   0.384   27,15   19,12   7,18    17,11   sig
16  a   b   chr16   -   0.9715  0.978   200141  3,2 260280  3,3 non_sig
21  a   b   chr1    +   0.9365  0.933   149118  1,12    133175  11,5    non_sig
22  a   b   chrX    +   0.6475  0.7265  129,57  58,35   104,78  37,29   test
26  a   b   chr3    +   0.05    0.0475  54,32   721503  46,27   519617  non_sig
27  a   b   chr3    +   0.0475  0.045   57,34   721503  47,30   519617  non_sig", header = F)


ggscatter(
    df,
    x= "V6", y= "V7",
    fill = "V12",
    shape = 21,
    size = 5,
    palette = c(test = "black", sig = "red", non_sig = "orange"))

在此处输入图像描述

于 2018-07-23T23:36:55.527 回答
1

假设第 12 列中的变量是因子，它们的默认排序是按字母顺序排列的。因此，在您的示例中，第一个调色板颜色（“黑色”）进入第一个因子级别（“non_sig”）；第二种颜色（“红色”）属于第二个因素（“sig”）。

如果要以不同方式分配颜色，则需要重新排序调色板中的因子级别或颜色名称。例如，要将“black”、“red”和“green”分配给因子“sig”、“non_sig”和“new_var”，您可以执行以下操作：

df$V12 <- factor(df$V12, levels = c("sig", "non_sig", "new_var"))

然后在情节中：

palette = c("black", "red", "green")
于 2018-07-23T23:40:44.777 回答