将箱形图数据分组，同时将其各自的X轴标签保留在R

Question

我正在根据样本类型（x轴）绘制Gene1计数（y轴）的等级。我希望根据其来源组织（乳房，结肠直肠，肺）和颜色代码对样本类型进行分组，分别来自癌症或红色和绿色的正常组织。

我制作了图表1）BOXPLOTS WITH FACETS（请参见下文），这与我的愿景很接近，但显示了一些主要问题。我在图表上有一些问题需要改进：

[IMG] http://i57.tinypic.com/10yfmmw.png[/IMG]

1）每个方面最终都有9个车道（列），其中许多车道没有被一个箱子占据。如何删除每个方面中未被框占据的通道（列）？

2）我是否可以在不使用构面的情况下绘制此图形，同时仍保留图中所示的分组？

3）是否可以创建两个方面标签？即我想贴标签＆＃34; Gene1＆＃34;在现有的facet标签之上和之上。这将使我能够为Gene2生成如下所示的相同图形，因此我可以使用＆＃34; master＆＃34;每个图表顶部的构面标签。

我希望这是有道理的。感谢大家的建议和想法。

请参阅以下代码，该代码可让您下载我的数据并重现图表：

测试文件导入

fileURL <- "https://dl.dropboxusercontent.com/u/4098921/testfile.csv"
test <- read.csv(fileURL,header=T)
head(test)


> head(test)
  Subset Tissue        Type id Gene1 Gene2
1 Normal Breast GTEx_Breast  1  5027 12597
2 Normal Breast GTEx_Breast  2  5287 12338
3 Normal Breast GTEx_Breast  3  2385 12543
4 Normal Breast GTEx_Breast  4  3174 12266
5 Normal Breast GTEx_Breast  5  6593 11350
6 Normal Breast GTEx_Breast  6  4648 10932

1）带有FACETS的BOXPLOT

library(ggplot2)
ggplot(test,aes(x=Type, y=Gene1, fill=Subset))+
geom_boxplot(notch=T, notchwidth=0.5,outlier.shape=1,outlier.size=2, coef=1.5)+
theme(axis.text=element_text(color="black"))+
theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4))+
theme(panel.grid.minor=element_blank())+
labs(size= "Type",x = "",y = "Rank of Gene1 count", title = "1) BOXPLOT WITH FACETS")+
scale_fill_manual(values=c("red","lawngreen"),name="Subset",
labels=c("Cancer (TCGA)", "Normal (GTEx)"))+
facet_grid(~Tissue)

Answer 1

您可以通过将scales = "free_x"添加到facet_grid来删除x轴上未使用的标签。如果您还添加space = "free"，您将获得相同大小的箱线图。用：

ggplot(test,aes(x=Type, y=Gene1, fill=Subset))+
  geom_boxplot(notch=T, notchwidth=0.5,outlier.shape=1,outlier.size=2, coef=1.5)+
  theme(axis.text=element_text(color="black"))+
  theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4))+
  theme(panel.grid.minor=element_blank())+
  labs(size= "Type",x = "",y = "Rank of Gene1 count", title = "1) BOXPLOT WITH FACETS")+
  scale_fill_manual(values=c("red","lawngreen"),name="Subset",
                    labels=c("Cancer (TCGA)", "Normal (GTEx)"))+
  facet_grid(~Tissue, scales = "free_x", space = "free")

你会得到以下情节：

如果您不想使用facet但保留分组，则可能需要创建一个新的变量来考虑分组。你可以用interaction：

来做到这一点

# create the new variable
test$newType <- factor(interaction(test$Tissue,test$Type))
# set the correct order of the new variable
test$newType <- factor(test$newType,
                       levels=levels(test$newType)[order(levels(test$newType))],
                       labels=levels(test$Type)[order(levels(test$newType))])

然后你可以用：

创建一个新的情节

ggplot(test,aes(x=newType, y=Gene1, fill=Subset))+
  geom_boxplot(notch=T, notchwidth=0.5,outlier.shape=1,outlier.size=2, coef=1.5)+
  theme(axis.text=element_text(color="black"))+
  theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4))+
  theme(panel.grid.minor=element_blank())+
  labs(size= "Type",x = "",y = "Rank of Gene1 count", title = "1) BOXPLOT WITH FACETS")+
  scale_fill_manual(values=c("red","lawngreen"),name="Subset",
                    labels=c("Cancer (TCGA)", "Normal (GTEx)"))

产生以下图：

如果你想在你的情节中同时包含Gene1和Gene2，那么最好的做法是首先将数据重塑为长格式：

library(tidyr)
test2 <- test %>% gather(gene,value,5:6)

你可以制作一个情节：

ggplot(test2,aes(x=Type, y=value, fill=Subset))+
  geom_boxplot(notch=T, notchwidth=0.5,outlier.shape=1,outlier.size=2, coef=1.5)+
  theme(axis.text=element_text(color="black"))+
  theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4))+
  theme(panel.grid.minor=element_blank())+
  labs(size= "Type",x = "",y = "Rank of Gene1 count", title = "1) BOXPLOT WITH FACETS")+
  scale_fill_manual(values=c("red","lawngreen"),name="Subset",
                    labels=c("Cancer (TCGA)", "Normal (GTEx)"))+
  facet_grid(gene~Tissue, scales = "free_x", space = "free")

上面的代码将为您提供以下图表：

在不使用其他构面图层的情况下包含不同基因的选项，您可以使用：

ggplot(test2,aes(x=Type, y=value, fill=Subset, alpha=gene))+
  geom_boxplot(notch=T, notchwidth=0.5,outlier.shape=1,outlier.size=2, coef=1.5)+
  theme(axis.text=element_text(color="black"))+
  theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4))+
  theme(panel.grid.minor=element_blank())+
  labs(size= "Type",x = "",y = "Rank of Gene count", title = "BOXPLOT WITH FACETS")+
  scale_fill_manual(values=c("red","lawngreen"),name="Subset",
                    labels=c("Cancer (TCGA)", "Normal (GTEx)"))+
  facet_grid(.~Tissue, scales = "free_x", space = "free") +
  theme_bw() +
  theme(axis.text.x=element_text(angle=90, vjust=0.5, hjust=1))

结果：