2015-05-18 39 views
0

我想根据示例代码here使用ggplot2中的geom_tile制作标准化读取丰度值的热图。我当前的代码为所需的范围生成了热图,但由于某些原因,7个范围中只有4个显示在热图中,我无法弄清楚问题所在。当我在原始链接中遵循示例时,它工作正常,所以我必须在代码中错误地更改了某些内容。任何人都可以请帮我找出导致这个错误的代码吗? 我想有以下颜色方案:传说不能正确显示在与ggplot2的热图中

-Inf < value <= 0 -> white 
    0 < value <=1 -> yellow 
    1< value <=10 -> orange 
    10< value <= 100 -> darkorange2 
    100< value <= 1000 -> red 
    1000 <value <= 10000 -> red3 
    10000 < value <= 32000 -> red4 

这里是我的代码:

#re-order the labels in the order of appearance in the data frame 
    df$label <- factor(df$X1, as.character(df$X1)) 
    # make the cuts 
    df$value1 <-cut(df$value,breaks=c(Inf,0,1,10,100,1000,10000,32000),right = T) 
    ggplot(data = df, aes(x = label, y = X2)) + geom_tile(aes(fill=value1), colour= "black") + scale_fill_manual(breaks=c("(-Inf,0]", "(0,1]", "(1,10]", "(10,100]", "(100,1000]", "(1000,10000]", "(10000,32000]"),values =c("white","yellow","orange","darkorange2","red","red3","red4")) 

这里是我的数据的预览(实际数据有228行每百万值配备读取38点的ID在6个不同的实验):

head(df) 
       X1  X2  value        label  value1 
     1 merged_read_17785-997_aka_156_aka_21 RPM.MT1 91.783028 merged_read_17785-997_aka_156_aka_21 (10,100] 
     2 merged_read_133362-79_aka_156_aka_21 RPM.MT1 6.403467 merged_read_133362-79_aka_156_aka_21  (1,10] 
     3 merged_read_147828-69_aka_156_aka_20 RPM.MT1 4.268978 merged_read_147828-69_aka_156_aka_20  (1,10] 
     4 merged_read_162443-60_aka_156_aka_21 RPM.MT1 0.000000 merged_read_162443-60_aka_156_aka_21 (-Inf,0] 
     5 merged_read_262156-32_aka_156_aka_21 RPM.MT1 5.691971 merged_read_262156-32_aka_156_aka_21  (1,10] 
     6 merged_read_22905-759_aka_159_aka_21 RPM.MT1 140.164780 merged_read_22905-759_aka_159_aka_21 (100,1e+03] 

这里是我从上面的数据得到的情节: heatmap

回答

0

我想我想通了这一点,如果我拿出从scale_fill_manual休息时间参数,那么所有的传说都显示:

ggplot(data = df, aes(x = label, y = X2)) + geom_tile(aes(fill=value1), colour= "black") + scale_fill_manual(values =c("white","yellow","orange","darkorange2","red","red3","red4"))