Monocle3 直接用 Seurat 的 Umap+Cluster (2)?

文章预览

生信益站，一点就有益 ！ 祝友友们天天开心，月月发 CNS~ 设我为🌟星标，每天就能第一时间看到推送啦~ Seurat + monocle3 第 2 弹 ~ ❝ 上一篇文章 《 Monocle3直接用Seurat的Umap+Cluster(1) 》， 其实， 没有把 seurat 的 cluster 、 gene_short_name 映射回 monocle3 的 cds 对象 ， 这篇文章提供了一种尝试 。 load packages library (Seurat) library (tidyverse) library (dplyr) library (monocle3) library (SeuratWrappers) library (ggplot2) library (ggridges) Read 10X dataset(feature barcode matrix) and create a Seurat Object data  < - Read10X(data.dir =  "aaa/tes.filtered_feature_bc_matrix" ) data  < - CreateSeuratObject(counts = data, project= 'aaa' , min.cells =  3 , min.features =  200 ) Check and Visualize QC metrics like mitochonrial RNA reads. data[[ "percent.mt" ]]  < - PercentageFeatureSet(data, pattern =  "^MT-" ) VlnPlot(data, features = c( "nFeature_RNA" ,  "nCount_RNA" ,  "percent.mt" ), ………………………………