2017-04-09 172 views
2

受多重因素影响小面我有一个data.frame,我试图在facetted方式与Rggplotgeom_boxplot绘制:颜色和ggplot

set.seed(1) 

vals <- rnorm(12) 
min.vals <- vals-0.5 
low.vals <- vals-0.25 
max.vals <- vals+0.5 
high.vals <- vals+0.25 


df <- data.frame(sample=c("c0.A_1","c0.A_2","c1.A_1","c1.A_2","c2.A_1","c2.A_2","c0.B_1","c0.B_2","c1.B_1","c1.B_2","c2.B_1","c2.B_2"), 
       replicate=rep(c(1,2),6),val=vals,min.val=min.vals,low.val=low.vals,max.val=max.vals,high.val=high.vals, 
       group=c(rep("A",6),rep("B",6)),cycle=rep(c("c0","c0","c1","c1","c2","c2"),2), 
       stringsAsFactors = F) 

在这个例子中,有两个因素我想facet

facet.factors <- c("group","cycle") 
for(f in 1:length(facet.factors)) df[,facet.factors[f]] <- factor(df[,facet.factors[f]],levels=unique(df[,facet.factors[f]])) 
levels.vec <- sapply(facet.factors,function(f) length(levels(df[,f]))) 

但在其他情况下,我可能只有一个或两个以上的因素。

有没有办法传递给facet_wrap要素的面向量和列数?

这里是我试过了,这里除了我创造了我自己的颜色对于每个因子水平:

library(RColorBrewer,quietly=T) 
library(scales,quietly=T) 
level.colors <- brewer.pal(sum(levels.vec),"Set2") 

require(ggplot2) 
ggplot(df,aes_string(x="replicate",ymin="min.val",lower="low.val",middle="val",upper="high.val",ymax="max.val",col=facet.factors,fill=facet.factors))+ 
    geom_boxplot(position=position_dodge(width=0),alpha=0.5,stat="identity")+ 
    facet_wrap(~facet.factors,ncol=max(levels.vec))+ 
    labs(x="Replicate",y="Val")+ 
    scale_x_continuous(breaks=unique(df$replicate))+ 
    scale_color_manual(values=level.colors,labels=unname(unlist(sapply(facet.factors,function(f) levels(df[,f])))),name="factor level")+scale_fill_manual(values=level.colors,labels=unname(unlist(sapply(facet.factors,function(f) levels(df[,f])))),name="factor level")+ 
    theme_bw()+theme(legend.position="none",panel.border=element_blank(),strip.background=element_blank(),axis.title=element_text(size=8)) 

这显然引发此错误:

Error in combine_vars(data, params$plot_env, vars, drop = params$drop) : 
    At least one layer must contain all variables used for facetting 

显然,这种工作原理:

ggplot(df,aes_string(x="replicate",ymin="min.val",lower="low.val",middle="val",upper="high.val",ymax="max.val",col=facet.factors,fill=facet.factors))+ 
    geom_boxplot(position=position_dodge(width=0),alpha=0.5,stat="identity")+ 
    facet_wrap(group~cycle,ncol=max(levels.vec))+ 
    labs(x="Replicate",y="Val")+ 
    scale_x_continuous(breaks=unique(df$replicate))+ 
    scale_color_manual(values=level.colors,labels=unname(unlist(sapply(facet.factors,function(f) levels(df[,f])))),name="factor level")+scale_fill_manual(values=level.colors,labels=unname(unlist(sapply(facet.factors,function(f) levels(df[,f])))),name="factor level")+ 
    theme_bw()+theme(legend.position="none",panel.border=element_blank(),strip.background=element_blank(),axis.title=element_text(size=8)) 

enter image description here

但它忽略了我传递的颜色,并没有添加图例,我想因为我无法通过vectorcolfillaesthetics,显然我必须对这个方面进行硬编码。

这并不无论是对磨制问题工作:

ggplot(df,aes_string(x="replicate",ymin="min.val",lower="low.val",middle="val",upper="high.val",ymax="max.val",col=facet.factors,fill=facet.factors))+ 
     geom_boxplot(position=position_dodge(width=0),alpha=0.5,stat="identity")+ 
     facet_wrap(facet.factors[1]~facet.factors[2],ncol=max(levels.vec))+ 
     labs(x="Replicate",y="Val")+ 
     scale_x_continuous(breaks=unique(df$replicate))+ 
     scale_color_manual(values=level.colors,labels=unname(unlist(sapply(facet.factors,function(f) levels(df[,f])))),name="factor level")+scale_fill_manual(values=level.colors,labels=unname(unlist(sapply(facet.factors,function(f) levels(df[,f])))),name="factor level")+ 
     theme_bw()+theme(legend.position="none",panel.border=element_blank(),strip.background=element_blank(),axis.title=element_text(size=8)) 

所以我的问题是:1。 有没有办法通过一个载体来facet_wrap? 2.是否有一种方法可以通过一个因子矢量而不是单一的因子进行着色和填充?

回答

3

我们不能指定两种颜色的着色/填充到单个框,我建议刻面变量粘贴在一起作为着色/填充比例:

df$col.fill <- Reduce(paste, df[facet.factors]) 

facet_wrapfacets同时接受字符向量或一个双边公式:

facet.formula <- as.formula(paste('~', paste(facet.factors, collapse = '+'))) 

因此,代码最后看起来是这样的:

ggplot(df, 
     aes_string(
      x = "replicate", ymin = "min.val", ymax = "max.val", 
      lower = "low.val", middle = "val", upper = "high.val", 
      col = "col.fill", fill = "col.fill" 
     )) + 
    geom_boxplot(position = position_dodge(width = 0), 
       alpha = 0.5, 
       stat = "identity") + 
    facet_wrap(facet.factors, ncol = max(levels.vec)) + 
    # alternatively: facet_wrap(facet.formula, ncol = max(levels.vec)) + 
    labs(x = "Replicate", y = "Val") + 
    scale_x_continuous(breaks = unique(df$replicate)) + 
    theme_bw() + 
    theme(
     #legend.position = "none", 
     panel.border = element_blank(), 
     strip.background = element_blank(), 
     axis.title = element_text(size = 8) 
    ) 

不显示图例是因为您添加了legend.position = "none",enter image description here

顺便说一句,如果你在你的代码添加一些空间和换行符它肯定会提高可读性。