SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Newbie needs help with DGElist for EdgeR analysis Germwise Bioinformatics 2 01-01-2013 02:46 PM
Help for a newbie in here please! andresroca Bioinformatics 3 12-14-2012 02:43 AM
Newbie in Bioinformatics hirohaha Introductions 0 09-09-2012 08:45 AM
Newbie Help roshanpatel95 Bioinformatics 5 10-31-2011 07:39 AM
hello from a newbie kathryn Introductions 0 08-13-2008 12:36 AM

Reply
 
Thread Tools
Old 01-24-2013, 11:51 AM   #1
Rbc
Junior Member
 
Location: USA

Join Date: Jan 2013
Posts: 2
Default newbie edgeR

Hello,
I am new to next Gen sequencing and R. We got the .bam files from the core and my boss wants me to find the differential expression of the genes in 2 time points. So I was wondering is it possible to do edgeR from the .bam or the rpkm files. If so can you please forward me to that website.
Thank you
Rbc is offline   Reply With Quote
Old 01-25-2013, 01:00 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

edgeR doesn't use the RPKM values as input, but rather the raw counts per gene. The normal work-flow is to run htseq-count on the various BAM files and then use the outputs in edgeR (or DESeq or one of the many other analysis programs). To do that, just invoke something like the following for each BAM file:

Code:
samtools view SOMEFILE.BAM | htseq-count - GTF > counts.txt
where SOMEFILE.BAM is a BAM file and GTF is whatever annotation GTF is appropriate for your experiments. If all of your BAM files are in the same directory, remember to change the name of counts.txt with every file. You may also need to change the options used with htseq-count (often -m or -i).

Also, there are a few lines of summary information printed at the end of the htseq-count input, you can remove those now in a text editor, or do so later in R (whatever you prefer).
dpryan is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:51 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO