Hi, I have the samples for 4 plant types (WT, KO1, KO2 and OE1) in three replicates, Mock treated and after stress treatment. So in total 24 samples (4 Mock and 4 Treated, 3 replicates of each). I would like to analyse the differential gene expression in each sample type before treatment (Mock) and after treatment. I used STAR for alignment of reads on reference genome. Then I obtained count.txt files for each sample by using the tool featureCounts. Now I would like to analyse the differential gene expression using edgeR. But I am having problem in creating the matrix as I have a large number of samples. Could you please help me in that?