I have a rather silly question.
I am using DESeq2 to analyze raw counts data from RNAseq data (totally 150 samples; 50 genotypes; 3 replicates for each genotype).
I create the DESeqDataSet object in R through
dds = DESeqDataSetFromMatrix(countData = counts_m, colData = samples_info, design = ~1)
In the samples_info file, I have information of which sample belongs to which genotype.
Now I want to calculate a mean value of all 3 replicates for each genotype, how could I do that? Does DESeq2 have a function to do this easily? I have googled around, but didn't find the answer.
Thanks a lot.