Some Thing Wrong with DESeq2 and Limma?

[hjh-air-cond]

HI! I've been using Libra for months, but today I found something wrong with the DESeq2 and Limma.

I'm running Libra on my dataset, each condition containing 3 replicates.

But When I using Deseq2, It returns: And When I using Limma, It returns:

Also, for now run_de() can work with edgeR. And the fun fact is the same code can work on the the other Seurat Object but not the one I'm analysing now. The difference between this dataset and others is this is bigger and containg other assays( including the 'intergrated' assay and the splicing assays for velocyto. )

Could someone offer a soulation? Thanks!

[jordansquair]

If you have multiple assays and you pass a Seurat object you need to correctly set the default assay, or else pass in the matrix directly to avoid errors. DefaultAssay(DEG.data) = 'RNA'

[hjh-air-cond]

Thanks for you reply so fast!

So do you mean the multiple assays cause that error?

I do try the DefaultAssay(), bu it seem don't work. and It still doesn't work when I remove all the other assays away from the Seurat Object in my code, which make me really confusing LoL.

[jordansquair]

Given there is no error with other datasets I do not think this is a bug with Libra. Perhaps try passing the matrix you want to do your analysis with directly in, instead of through the Seurat object. Alternatively, you can check the Seurat github for tutorials on assay manipulation.

[hjh-air-cond]

Sure, I will try your advice! Thanks for all the help!