Hi all,
I am fairly new in using R for loading custom tracks into UCSC genome browser.
Here is my code and error:
The error said:
Does anyone know why? Also, I cannot find the GRangesForGenome function... Thanks!
I am fairly new in using R for loading custom tracks into UCSC genome browser.
Here is my code and error:
Code:
library("rtracklayer")
y0=scan("coverage.txt")
wins=import.bed("ref.bed")
peakTrack<-GenomicData(wins$ranges,chrom="chr16",score=y0,genome="hg19")
session=browserSession("UCSC")
track(session,"y0")=peakTrack
region <- range(peakTrack[,])
view <- browserView(session, region, full=c("y0"))
Code:
Error in GRangesForGenome(genome, chrom, range, seqinfo = seqinfo(session)) : 'genome' must be a single string identifying a genome