Search
Question: How to get UCSC url when use rtracklayer, contain the track in the url ,it may be very useful
0
gravatar for zhangqshxxzz
8 days ago by
zhangqshxxzz0 wrote:

Dr Michael Lawrence:

Hellp, I use the rtracklayer to product the UCSC show with track, it work very well. But if the aim region are too much, is there have a option can direct get the UCSC url with the track, not the opened HTML........

It will very convenience to annotation the region file in bath, when use for cnv discover and annotation.

 

Thanks!!

 

JOSH ZHANG

ADD COMMENTlink modified 7 days ago by James W. MacDonald45k • written 8 days ago by zhangqshxxzz0

region<-c("chr2:240029740-240037011")
session <- browserSession()
genome(session) <- "hg19"
export(region, "Track12.bed")
restoredTrack2 <- import("Track12.bed")
track(session, "Track") =  restoredTrack2
browserView(session,restoredTrack2*-10,hide=c("UCSC Genes","GTEx Gene","ENCODE Regulation...",
                                              "Conservation","Common SNPs(150)","RepeatMasker"),
            pack=c("Track","DECIPHER","RefSeq Genes","ClinGen CNVs","ClinVar Variants",
                   "GeneReviews","OMIM Genes","DGV Struct Var")  )

and also when I RUN the command ,it is work well before , now it report the error like this....

Error in names(trackIds) <- sub("^ ", "", sapply(nodes, xmlValue)) : 
  'names' attribute [211] must be the same length as the vector [209]
In addition: Warning message:
In readChar(f, 1) : can only read in bytes in a non-UTF-8 MBCS locale

ADD REPLYlink written 8 days ago by zhangqshxxzz0

Should be fixed as of rtracklayer 1.38.2.

ADD REPLYlink written 5 days ago by Michael Lawrence9.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 307 users visited in the last hour