I know that we can use differentialGeneTest to find genes that are most differentially expressed across the states by using the formula "~State". However, in the resulting data frame, there is no indication of the state number w.r.t. a given gene. Is it possible to get an extra column with a state number for a gene, essentially indicating which state a given gene represents as a marker? Something like the function "FindMarkers" outputs in the Seurat package?
Thanks in advance for any replies.
P.S: Original post