This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Merge featurecounts lists

I have a very large RNA-Seq sample size (n=100), so while performing quantification through featurecounts, Rstudio crashed. So, I decided to perform the analysis in chunks (n=20 bam files). Hence, I generated 5 lists through featurecounts (fc1,fc2....). Now I tried to merge all these lists in one list fc, but the dimension of fc is list of 20 (ideally should have been- list of 4=> annotation;counts;target;stats).

Functions I used to merge lists- c(fc1,fc2), mapply(fc1,fc2), rbind(fc1,fc2), Map(list,fc1,fc2).

What mistakes I might be doing? Any alternative for merging lists so that I get only one list (ideal list of 4)? Any other way so that I can process all files at once in R?

merge-lists featurecounts rstudio dimension
• 1,405 views • link
written by imaparna27 • imaparna27

I would suggest running featureCounts on the command line and using multiple threads (depending on the local machine you are using). If you have access to a server, even better. That should probably avoid it from crashing and, this way, directly obtaining the full count matrix.

• 0 views • link
written by Marco Pannone • Marco Pannone

I was already using 20 threads in Rstudio. Although, I didn't realise command line option at first. Thanks a lot.

• 0 views • link
written by imaparna27 • imaparna27

0 answers

No answers yet.

Log in to answer this question.

More posts like this