小编aur*_*ien的帖子

当使用broom我用结合dplyr::group_by,并dplyr::do进行分组数据由于上行动@drob.例如,根据齿轮系统将线性模型拟合到汽车:

library("dplyr")
library("tidyr")
library("broom")

# using do()
mtcars %>%
  group_by(am) %>%
  do(tidy(lm(mpg ~ wt, data = .)))

# Source: local data frame [4 x 6]
# Groups: am [2]

#     am        term  estimate std.error statistic      p.value
#   (dbl)       (chr)     (dbl)     (dbl)     (dbl)        (dbl)
# 1     0 (Intercept) 31.416055 2.9467213 10.661360 6.007748e-09
# 2     0          wt -3.785908 0.7665567 -4.938848 1.245595e-04
# 3     1 (Intercept) 46.294478 3.1198212 14.838824 1.276849e-08
# 4     1          wt -9.084268 1.2565727 -7.229401 …

使用这个很棒的ggplot包,我想要一个条形图,其中fill美学被映射到一个连续变量,实际上是qvalues,在它上面是一个纹理,实际上是条纹和交叉阴影.

颜色渐变很重要,因为它代表重要性,而纹理则显示类别"A","B"及其重叠.所以对角线分别以一种方式,相反的方式和交叉的阴影线.我知道维恩图可以完成这项工作,但我认为我们更容易看到35个样本和比较.

这fill是微不足道的,然而,纹理是棘手的.在SO上,感谢@baptise(见这里和这里)我设法用虚拟数据得到这个输出:

问题是十字架背景的透明度.如果有可能alpha在grid.patternFill()功能上的后台工作那将是伟大的.不幸的是,在我手中它不起作用.

任何帮助都感激不尽.

这里可以加载虚拟数据:

dfso <- structure(list(Sample = c("S1", "S2", "S1", "S2", "S1", "S2"), 
  qvalue = c(14.704287341, 8.1682824035, 13.5471896224, 6.71158432425, 
  12.3900919038, 5.254886245), type = structure(c(1L, 1L, 2L, 
  2L, 3L, 3L), .Label = c("A", "overlap", "B"), class = "factor"), 
  value = c(897L, 1082L, 503L, 219L, 388L, 165L)), class = c("tbl_df", 
  "tbl", "data.frame"), row.names = c(NA, -6L), .Names = c("Sample", 
  "qvalue", "type", …