r - how to add the third legend in ggplot graph -


the problem want add multiple legend different layer of layoutkaryogram function in ggbio plot. 

library(genomicranges) library(ggbio) df<-data.frame(chromosome=c(9,9,9,9,9,9),start=c(12229,12244,12245,12248,19576961,19792417)            ,end=c(12229,12244,12245,12248,19576961,19792417),sample1=c("kdnl","ct","ir","ir","dh","dh"),            marker=c('snp',"snp",'snp',"snp",'ssr',"ssr"))  myideo<-data.frame(chromosome=9,start=1,end=21757032)

then change granges object

> dfgr  granges object 6 ranges , 2 metadata columns:       seqnames               ranges strand |   sample   marker          <rle>            <iranges>  <rle> | <factor> <factor>   [1]        9 [   12229,    12229]      * |     kdnl      snp   [2]        9 [   12244,    12244]      * |       ct      snp   [3]        9 [   12245,    12245]      * |       dh      snp   [4]        9 [   12248,    12248]      * |       dh      snp   [5]        9 [19576961, 19576961]      * |     kdnl      ssr   [6]        9 [19792417, 19792417]      * |     kdnl      ssr   -------   seqinfo: 1 sequence unspecified genome; no seqlengths   > myideogr granges object 1 range , 0 metadata columns:       seqnames        ranges strand          <rle>     <iranges>  <rle>   [1]        9 [1, 21757032]      *   -------   seqinfo: 1 sequence unspecified genome; no seqlengths

and try plot , give me nice graph

q1 <- ggplot() + layout_karyogram(myideogr, cytoband = false)  q2=q1 + layout_karyogram(dfgr[dfgr$marker == "snp",], geom ="rect", ylim =c(0, 10), aes(fill=sample,colour=sample)) +   layout_karyogram(dfgr[dfgr$marker == "ssr",], ylim =c(11, 19), aes(fill=sample,colour=sample)) +   scale_color_manual(values=c("#22b700","#00c0ba", "#f564e3","#ff6c91" ),name="snp",breaks=c("ct","ir","kdnl")) +   scale_fill_manual(values=c("#22b700","#00c0ba", "#f564e3","#ff6c91"),name="ssr",breaks=c("dh","kdnl")) +   guides(fill = guide_legend(order=2), color = guide_legend(order=1))

however, have information want add plot

gene_overlapgr    granges object 7 ranges , 2 metadata columns:       seqnames               ranges strand |      sample       trait          <rle>            <iranges>  <rle> | <character> <character>   [1]        9 [10128837, 10131086]      * |    sample_5         mxa   [2]        9 [10128837, 10131086]      * |    sample_5         mxa   [3]        9 [10128837, 10131086]      * |    sample_5         mxa   [4]        9 [10128837, 10131086]      * |    sample_5         mxa   [5]        9 [10338558, 10340355]      * |    sample_5         mxa   [6]        9 [10338558, 10340355]      * |    sample_5         mxa   [7]        9 [10328808, 10333199]      * |    sample_5         mxa

i try add new plot code below legend 

q2+ layout_karyogram(gene_overlapgr, geom ="rect",ylim=c(-10,-1),aes(fiil=trait))

any suggestion add third legend mxa lebeling?


-

Comments

Post a Comment

Popular posts from this blog

c# - Binding a comma separated list to a List<int> in asp.net web api -

the default binder in web api expecst http://url.com/webapi/report/?pageids=3243&pageids=2365 to bind public ihttpactionresult report(list<int> pageids){ // exciting webapi code} i wish bind http://url.com/webapi/report/?pageids=3243,2365 as running out of space in url get . i have created class public class commaseparatedmodelbinder : system.web.http.modelbinding.imodelbinder { public bool bindmodel(httpactioncontext actioncontext, modelbindingcontext bindingcontext) { //binding in here } } and registered in webapiconfig.cs var provider = new simplemodelbinderprovider( typeof(list<int>), new commaseparatedmodelbinder()); config.services.insert(typeof(modelbinderprovider), 0, provider); i have altered method signature use model binder so public ihttpactionresult report( [modelbinder] list<int> pageids){ // exciting webapi code} however break point in binder not being ...
Read more

how to prompt save As Box in Excel Interlop c# MVC 4 -

i want export excel sheet using excel interloop... i have datatable convert excel sheet... now there way save file using save dialogue box in mvc 4 here code public string datatabletoexcel(datatable dt, string htmlheading, string key) { try { microsoft.office.interop.excel.application excel; microsoft.office.interop.excel.workbook excelworkbook; microsoft.office.interop.excel.worksheet excelsheet; microsoft.office.interop.excel.range excelcellrange; excel = new microsoft.office.interop.excel.application(); excel.visible = false; excel.displayalerts = false; excelworkbook = excel.workbooks.add(type.missing); excelsheet = (microsoft.office.interop.excel.worksheet)excelworkbook.activesheet; excelsheet.name = "sheet1"; dataset ...
Read more

xslt 1.0 - How to access or retrieve mets content of an item from another item? -

my use case have item has link in item. example, item 123456789/152 has metadata field dc.relation.hasversion=123456789/717 . within item view of 123456789/152 , how can retrieve values of metadata of 123456789/717 ? example, want retrieve dc.language.iso value of 123456789/717 123456789/152 . i looked @ related items feature in dspace don't know how metadata displayed in related items list pulled items in list. i using dspace version 5.3 mirage 2 theme . edit based on schweerelos answer, below actual code. using custom metadata field dc.relation.languageversion . basically, want link other version instead of displaying value, display dc.language.iso of other version in item-view.xsl . i've incorporated answer of schweerelos question in code just displaying value of dc.relation.languageversion . sample values of dc.relation.languageversion 10665.1/9843 ; 10665.1/9844 ie not complete uri handle. thanks in advance! actual code <xsl:template name=...
Read more