why Error in medianExpr[runGenes, kpClust] : subscript out of bounds

[peachgong]

Sorry, the previous question is closed too soon. So I open a new issue regarding it.

I actually choose some very reasonable marker genes to start. Those are the genes I got from Seurat analysis. But it seems they are not included in the median table. How can it be?

[jeremymiller]

You must have either (1) had a different starting set of genes in your Seurat analysis as compared to what is in the medianExpr or (2) subset one of the two data matrices at some point. You can always right something like: runGenes <- intersect(runGenes, rownames(medianExpr)) to avoid this issue.