Controlling for covariates in DE analysis of single cell rna seq data

Question

[saeedabdurahman]

Dear Community,

I have human data from different ages and gender. After using integration with seurat, how would I best control for these confounding factors during differential gene expression. I see the option of latent.vars in FindMarkers function. Can I give latent.vars = c("Age", "gender") to account for both together? or can I only use one at a time?

Is there alternative packages to do the test better? Kindly advice, Kind regards, Saeed

[ATpoint]

Can you give some details what you are comparing and what the setup is?