我对R语言(和ggplot2)的工作有点新意。我正在尝试复制这张地图的设计:
我已经在y轴上列出了某些基因的热图,我还使用以下方法对基因进行了分组:
<div class="row r1">
<div class="square o1">
<div class="isq i1">adsfas</div>
<div class="isq i2">adfasd</div>
<div class="isq i3"></div>
<div class="isq i4"></div>
<div class="isq i5"></div>
<div class="isq i6"></div>
<div class="isq i7"></div>
<div class="isq i8"></div>
<div class="isq i9"></div>
</div>
<div class="square o2">
<div class="isq i1">adfasdf</div>
<div class="isq i2"></div>
<div class="isq i3"></div>
<div class="isq i4"></div>
<div class="isq i5"></div>
<div class="isq i6"></div>
<div class="isq i7"></div>
<div class="isq i8"></div>
<div class="isq i9"></div>
</div>
<div class="square o3">
<div class="isq i1"></div>
<div class="isq i2"></div>
<div class="isq i3"></div>
<div class="isq i4"></div>
<div class="isq i5"></div>
<div class="isq i6"></div>
<div class="isq i7"></div>
<div class="isq i8"></div>
<div class="isq i9"></div>
</div>
</div>
在我的长数据框架上:'Stuff' variable which contains the data to be plot
可重复的例子是:
interaction(Gene_Type, Gene)
完整的绘图代码是:
Sample<- c(3B.A9HV.01,3B.A9HX.01,3B.A9HY.01,3B.A9HZ.01,3B.A9I0.01,3B.A9I1.01,3B.A9I3.01)
Gene<- c(RB1,RB1,RB1,RB1,RB1,RB1,RB1)
v<- c(0,0,1,1,0,1,0)
Gene_Type<- c(SMG,SMG,Cell Cycle,SMG,Cell Cycle,SMG,SMG)
stuff<-data.frame(Sample, Gene, v, Gene_Type)
产生这个:
正如你所知道的那样,我的图表已经不在我想要的地方了。根据我的需要,基因按其gene_type分组。但是,由于某种原因,标签现在是Gene.Gene_Type而不仅仅是基因名称。此外,我在模型图的另一侧错过了一个显示分组的栏。是否有一些美学功能允许我像模型图一样对基因进行亚组?或者是覆盖和创建我自己的标签的唯一解决方案?