Advertisement
Guest User

Untitled

a guest
Feb 17th, 2016
100
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 1.02 KB | None | 0 0
  1. # replicate highlight box reverseStrand error
  2.  
  3. transcriptID <- "ENST00000395443"
  4.  
  5. mart <- biomaRt::useMart(
  6.   biomart = "ENSEMBL_MART_ENSEMBL",
  7.   dataset = "hsapiens_gene_ensembl",
  8.   host = "dec2013.archive.ensembl.org")
  9.  
  10. # Create transcript track
  11. biomartTrack <- Gviz::BiomartGeneRegionTrack(
  12.   filters = list(
  13.     ensembl_transcript_id = transcriptID),
  14.   biomart = mart)
  15.  
  16. # Create axis track
  17. axisTrack <- Gviz::GenomeAxisTrack()
  18.  
  19. # works as expected
  20. Gviz::plotTracks(
  21.   c(biomartTrack, axisTrack),
  22.   chromosome = "chr3")
  23.  
  24. # works as expected
  25. Gviz::plotTracks(
  26.   c(biomartTrack, axisTrack),
  27.   chromosome = "chr3",
  28.   reverseStrand = TRUE)
  29.  
  30. # Create highlight track
  31. highlightTrack <- Gviz::HighlightTrack(
  32.   trackList = biomartTrack,
  33.   start = biomartTrack@start,
  34.   width = 2000,
  35.   chromosome = "chr3")
  36.  
  37. # works as expected
  38. Gviz::plotTracks(c(highlightTrack, axisTrack))
  39.  
  40. # this only reverses the transcript, not the highlighttrack
  41. Gviz::plotTracks(
  42.   c(highlightTrack, axisTrack),
  43.   reverseStrand = TRUE)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement