r - 来自 ggplot2 geom_raster() 的双色矩形

Question

我的问题类似于此链接。在此示例中，1 表示突变，0 表示野生型，NA 表示不可用。我的数据框设置相同，但是每个基因可能包含两种或多种类型的突变。我想生成一个类似的图形，除非一个基因有两种类型的突变，我想把正方形切成两半，两种类型的突变都涂上颜色，类似于这个例子。目前，如果受试者中的基因有两个突变，则第二个突变会覆盖第一个突变。提前感谢您抽出时间提供帮助。

dat <- expand.grid(gene=1:10, subj=1:50)
dat$mut <- as.factor(sample(c(rep(0,300),rep(1,200)),500))
dat$mut[sample(500,300)] <- NA
dat[501,] = c(10,50,1) #included from comment below
ggplot(dat, aes(x=subj, y=gene, fill=mut)) +
  geom_raster() +
  scale_fill_manual(values = c("#8D1E0B","#323D8D"), na.value="#FFFFFF") +
  scale_x_discrete("Subject") +
  scale_y_continuous(breaks=1:10,
    labels=c("D0","D1","D2","D3","D4","D5","D6","D7","D8","D9")) +
  guides(fill=FALSE) +
  theme(
    axis.ticks.x=element_blank(), axis.ticks.y=element_blank(),
    axis.text.x=element_blank(), axis.text.y=element_text(colour="#000000"), 
    axis.title.x=element_text(face="bold"), axis.title.y=element_blank(),
    panel.grid.major.x=element_blank(), panel.grid.major.y=element_blank(),
    panel.grid.minor.x=element_blank(), panel.grid.minor.y=element_blank(), 
    panel.background=element_rect(fill="#ffffff")
  )

Answer 1

我看不到您的数据曾经有一个主题和基因的多个条目？什么都没有被覆盖，因为没有什么可以覆盖的。

我添加了最后一行的重复，但将 mut 更改为 1 以显示。我还从栅格更改为平铺，并更改了不透明度，以便具有多个值的平铺具有不同的颜色。

如果您想要链接到的绘图之类的东西，则还需要创建移位和高度向量，就像该帖子显示的那样，因此每个图块再次被细分。

    dat[501,] = c(10,50,1)

    ggplot(dat, aes(x=subj, y=gene)) +
      geom_tile(alpha=.5,aes(fill=mut), show.legend = F) +
      scale_fill_manual(values = c("#8D1E0B","#323D8D"), na.value="transparent") +
      scale_x_discrete("Subject") +
      scale_y_continuous(breaks=1:10,
                         labels=c("D0","D1","D2","D3","D4","D5","D6","D7","D8","D9")) +
      theme(
        axis.ticks.x=element_blank(), axis.ticks.y=element_blank(),
        axis.text.x=element_blank(), axis.text.y=element_text(colour="#000000"), 
        axis.title.x=element_text(face="bold"), axis.title.y=element_blank(),
        panel.grid.major.x=element_blank(), panel.grid.major.y=element_blank(),
        panel.grid.minor.x=element_blank(), panel.grid.minor.y=element_blank(), 
        panel.background=element_rect(fill="#ffffff")
      )