Hi there,

I'm at a beginner level using R/Bioconductor. I'm currently trying to read in RNA-seq bed files from local and web (UCSC) sources and then hand over the score of every chromosome separately to a wavelet transform (packages "Wavelets", "WaveletComp"). Could you give me a brief instruction how to proceed? I would be very grateful for your help.

Cheers! SNP

The Bioconductor package rtracklayer provides for reading in bed files.

source("https://bioconductor.org/biocLite.R")

biocLite("rtracklayer")

library("rtracklayer")

?BEDFile

I'm not exactly sure what you mean by score of every chromosome, is that something that is provided in the file or is it something that requires calculation?