This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to calculate gene expression level from BW and BED file

I am very new to the field.

I have a new dataset , GRO-seq from GEO site, and it only consist of raw files (both positive and negative strand measured at different time) in bw format and bed format. I am not sure how to process/annotate/quantify the dataset. I want to transform the data such that it would be a r object, one column for probe id, and a column for expression level. I've heard that I'll need to use rtracklayer, but I'm not exactly sure what I'm suppose to do. And what is the process called, transforming from the raw dataset to gene expression.

rna-seq r

1 answer

All bioconductor packages have extensive documentation and tutorials. Please have a read: https://bioconductor.org/packages/release/bioc/vignettes/rtracklayer/inst/doc/rtracklayer.pdf https://bioconductor.org/packages/release/bioc/html/rtracklayer.html

rtracklayer is part of a larger ecosystem, a very good online course by the main author of a lot of these packages (Hansen): https://kasperdanielhansen.github.io/genbioconductor/

Good luck.

Log in to answer this question.