Hi , I'm trying to perform rrbs dme analysis using edgeR on 50 samples (5 replicates per timepoint /sex). It seems to be running forever (>8 hours) and not produce any output. Design matrix has 50 samples , 10 conditions. Is there anyway I can speed up the dispersion estimation?
y <- estimateDisp(y, design, trend="none",subset=5000)