I am using DESeq2 for differential gene expression. I have done alignment using STAR aligner and GRCh38 reference genome. I have generated count files using HTSeq. I have 3 controls and 4 cases and the size of each count file varies from 8 GB to 11 GB which makes to total 73 GB of data (7 Count files).
ddsHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable, directory = directory, design = ~ condition)
directory where I kept the count files.
This command is taking more than 12 hours and consumed memory of 46 GB. I have a maximum memory limit of 64 GB in my System.
Could anyone help how to solve this problem? ![Screenshot is given here...] : https://ibb.co/3s48HRw