I am still deeply irritated by the fact, that limma also calculates statistics for a one-sample t-test (so, no contrasts & designs involved) with a n=1 datarow - should this not be a classical case for a "NA"? #limma #data

@haojiawu @biorxivpreprint this is very interesting. I have been telling people that one of the main reasons for using #RStats instead of #python for gene expression analysis, is the lack of methods such as #DESeq2 or #limma as python libraries. This might tip the balance for a lot of people.

Of course there are many reasons to prefer R, the excellent #Bioconductor ecosystem one of them, and in fairness, for #scRNA analysis python has very strong ecosystem and community.

Wow, this paper on #scRNAseq and differential expression methods is an eye opener. https://www.nature.com/articles/s41467-021-25960-2

Nice overview of methods and relatively easy to understand explanation of what is wrong with certain, especially single-cell-specific, methods.

Of great help in my own scRNA-seq efforts. Should probably have read this earlier. Now, back to R I go. 😅

#statistics #rnaseq #limma #edgeR #DESeq2