Programs used in RNA-seq workflow
2
0
Entering edit mode
@cirisarri95-10027
Last seen 6.1 years ago

Good morning everyone, 

I was interested in knowing which programs is the bioconductor package RNA-seq workflow: gene exploratory analysis and differential expression for:

-Quality control

-Mapping

-Counting

-Read Analysis

Thanks everyone!

Regards

deseq2 bioconductor workflows differential gene expression rnaseqgene • 2.2k views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 19 hours ago
United States

hi,

I'm not sure I follow your question.

There are a number of R packages used in the workflow.

Every time you see code, such as library("Rsamtools"), that is loading an R package called, "Rsamtools". We describe each package used in each respective section, and when there are multiple options, we describe the one we use as well as provide links to alternatives (counting, differential expression analysis, distance calculation, etc.).

You can just search the page here for "library":

http://www.bioconductor.org/help/workflows/rnaseqGene/

And at the end of the page you can see all the libraries which were loaded throughout the R session, under "other attached packages".

Additionally we describe the use of the STAR read aligner in the workflow, with citation and a link to the STAR aligner software. This software is not an R package, but a standalone tool which is run on the command line.

ADD COMMENT
0
Entering edit mode

Hi,

Then, if I want to run the whole package workflow, do I need to first run STAR read aligner? Or will it be runned with the rest of the package?

Thanks.

ADD REPLY
0
Entering edit mode
STAR is not actually part of the workflow which is run "live". The live workflow contains only commands in R. We just show you example commands that were used for STAR to make the BAM files used in the workflow. You should prepare the alignment yourself using whatever alignment software you prefer to generate the BAM files.
ADD REPLY
0
Entering edit mode
@cirisarri95-10027
Last seen 6.1 years ago

Hi,

Thanks for your help. 

Also, for running this workflow, is it set up so that I input data once and it runs the whole process with no more help? Or should I be inputting for each step the data again and again?

ADD COMMENT
1
Entering edit mode

You only build the DESeqDataSet (dds) object once, and then the different functions do something and either return the object with new information stored inside of it, or they return a different object.

You'll have to just look at the workflow and follow along. You'll see by the names of the R variables when things are created, when they are modified and when you get a new object as an output of a function.

If you are brand new to R, you may benefit from taking a beginner's R course. I have some R resources collected here: 

http://genomicsclass.github.io/book/pages/resources.html

ADD REPLY

Login before adding your answer.

Traffic: 196 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6