Entering edit mode
6 weeks ago
Chris ▴ 70
Hi admins and all,
I try to run FindConservedMarkers() but I got this message:
Warning: Identity: 8 not present in group B. Skipping VVWarning: Identity: 8 not present in group A. Skipping NCError in marker.test[[i]] : subscript out of bounds.
This error appears in many clusters I chose. Would you have a suggestion? Thank you so much!
What is the actual line of code your are running? Your function parentheses are empty.
Hi @seidel. I run this one:
Have you checked to make sure that you have cells from each condition in each cluster?
Sorry if this is a basic question but how can I check that? I would like to perform pseudo-bulk RNA-Seq but I still struggle to get the count matrix. My seurat objects haven't had sample column and also cell-type annotated yet. Admins on this site are highly experienced so how can my questions get their attention? I could not tag any admins, right?
Admins are just people who've answered lots of questions and gaze over questions that may interest them and that they may know the answer to. You have their attention. Often people don't know the answer directly, so much as they know they can figure it out, but they need the pressure of an impending deadline to come here and seek distraction by solving puzzles for relief. :)
Did you create your Seurat objects, or were they handed to you from an unknown source? If you're following some workflow such as one of the Seurat Tutorials, you can access the metadata.table associated with your cells. For instance:
Thank you for the info! I created my own seurat object from features.tsv.gz, barcode.tsv.gz, and matrix.mtx.gz. And unlike the sample data which is well-annotated cell type, my seurat object is like this
I can get the info as you wrote above 🙂. However my data is not balance for each condition: 17900 vs 5579.
Balance may or may not matter. Sometimes people sample down (randomly select 5579 cells from the 17900 data set) to match. However, this means you have two seurat objects, correct? And you are trying to integrate them and find conserved markers, right? Have you successfully integrated them to a single object? And if so, do you have cell counts per condition?
Thank you for following up on my question. Yes, 2 seurat objects. I tried to find conserved markers. For a reason, I used merge function to integrate into a single object. Cell counts per condition are 17900 vs 5579.
Could other admins have a look?