DESeq2 input file from featureCount
Entering edit mode
12 months ago

Hello. I am analysing RNA-Seq data for 8 samples (2 treatments, 4 repetition each). I just finished the gene count using featureCount (Subread). How do I generate a matrix to use on DeSeq2 from the files (.RDS) genereted from featureCount? Do I just create a file with the gene_id and the counts? Also, what do I have to put (what are the columns) in the Metadata file? Thanks in advance.

Deseq2 matrix featureCount • 794 views
Entering edit mode


Assuming that you're running featureCounts inside R, a similar question was already answered:

Output from featureCounts() as input to DESeq2

I used previously featureCounts command-line tool, but not inside R. The output of the command-line tool is a text file per sample or one for all the samples, depending on how you run the tool, with several columns (with gene id, read-length, genomic coordinates, counts, etc).

From your description above, I'm assuming that the output of featureCounts R function is a rds file (read more about rds files in R). If so, you can import this rds file into R by doing the following:

counts <- readRDS("path_to/rds_file/file.rds")

Where "path_to/rds_file/file.rds" is a string to the output .RDS file that you have. This will import this object back to R. Although I'm not familiar with the data structure of this and you may need to parse this data in order to provide it to DESeq2.

Perhaps looking into the previous link on Biostars that answers to a similar question is the best way to go.

I hope this helps,


Entering edit mode

Yes, I'm using R. I imported the count table to excel already. Thank you for your answer. I'll check the other post.


Login before adding your answer.

Traffic: 2292 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6