0
$\begingroup$

I have a lists of DEGs identified by DESeq2, EdgeR and Limma. I would like to correlate the the gene rankings in the lists to decide on a package to use in downstream analysis. I am havig a few problems with this, would anyone be able to help?

Thanks

$\endgroup$
  • $\begingroup$ Hello Lyd, could you be a bit more specific? I think all these software have their advantages and disadvantages, it's impossible to give you a good advice without more information. Also, did you tried to figure it by your own? What are the results of your research and why you are not convinced yet? $\endgroup$ – Kamil S Jaron Aug 8 at 20:16
  • $\begingroup$ Im doing differential expresion analysis in R. Im comparing KD samples rto control samples using all 3 packages above. It was suggeested to me to do a correlation test between the results but Im no t sure what way to go about it. Not sure what the inputs should be! $\endgroup$ – Lyd Aug 9 at 8:42
1
$\begingroup$

I would use the p-value/FDR which each method returns to rank the gene list for that method in order from 'most likely to be DE' to least likely. You would then have three ranked lists of genes - you would put the actual rankings in for each gene - 1 for the top ranked, 2 for the next, etc. You could then carry out, say Spearman's rank correlation on the lists of rankings to assess agreement between the methods.

You could also decide to keep all 3 methods, and accept genes as DE if, for example, at least 2 out of 3 methods give an acceptable p-value/FDR.

From my experience, correlation should be reasonably high at the top of the list, where the methods tend to agree more on the rankings, and when the correlation gets lower, further down the list, would be where I would start to question whether the genes are really differentially expressed.

$\endgroup$
  • $\begingroup$ This worked! Thank you :) $\endgroup$ – Lyd Aug 15 at 12:04

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.