This is a test version of Biostars. For the public version, visit https://www.biostars.org.
FIltering Seurat clusters

Hello I am trying to pull out CD8 T cell clusters from my seurat object. I am using the code

CD8 <- CH005@meta.data[“seurat_clusters” %in% c(1,2,3,7,8),]

But it keeps resturning an empty table. Does anyone have any suggestions? Thanks!

rna-seq seurat r 10x
• 4,715 views • link
updated by rpolicastro rpolicastro • written by roberts • roberts

1 answer

Instead of accessing the object slots directly, you can use the subset function.

CD8 <- subset(CH005, subset=seurat_clusters %in% c(1, 2, 3, 7, 8))
• 0 views • link
written by rpolicastro rpolicastro

Log in to answer this question.

More posts like this
  • retaining only the clusters of interest
    written by shamza • shamza
    Hey I am trying to subset only the clusters of my interest by the following command ```r diff2 <- subset(merged_seurat_filtered, idents = c(0:16),) ``` There …
  • Splitting a Seurat Object HTO classification
    written by ecl75 • ecl75
    Hello, I am trying to split a Seurat object based on HTO classification. I have 3 different HTO based on different stimulation conditions. I want …
  • Gini Index Box plot
    written by roberts • roberts
    Hello, I want to make a box plot of Gini indices for each CD8 T cell clusters/CD4 T cell clusters showing significant clonal expansion within …
  • Issue renaming clusters in Python
    written by V V
    Hello, I've analysed my single cell RNAseq dataset in Seurat (v4). I am now wanting to run velocity analysis on it using CellRank (in python). …
  • Adding UMAP coordinates to metadata
    written by roberts • roberts
    Hello, I want to store my umap coordinates in my seurat objects meta.data. The umap coordinates are saved in my Object under reductions in cell.embeddings. …
  • Top 10 clonotype frequency
    written by roberts • roberts
    Hello I am using top10_clonotype <- table(AD007S$clonotype_id) top10_clonotype <- as.data.frame(top10_clonotype) top10_clonotype <- slice_max(top10_clonotype, order_by = top10$freq, n=11) I am trying to pick out the top …
  • Mapping clonotypes on UMAP
    written by roberts • roberts
    Hello! I have a suerat object and I would like to map the top 10 clonotypes in it. I am using the code top10 <- …
  • Two different library preps
    written by roberts • roberts
    Hello I have to combine two data sets from 10x in seurat. one data set was prepared following 5' library prep and the other was …
  • pulling out TCR data
    written by roberts • roberts
    Hello! I have successfully loaded my VDJ data into my seurate object. Once I have identified the cluster of CD4 and CD8 I want to …
  • Analyzing TCR clonotype in Seurat
    written by roberts • roberts
    Hello, I have merged my VDJ data into my Seurat object by adding it as metadata. Now is there a way to for example pull …