R绘图3
2021-07-21 本文已影响0人
rong酱
做喜欢做的事情 心情还是好多啦
library(Seurat)
library(dplyr)
library(ggplot2)
library(Matrix)
library(reshape2)
dpi <- 300
data.integrated <- readRDS('Merge.Rdata')
cbmc <- FindNeighbors(data.integrated)
head(cbmc)
#levels(Idents(data.integrated))
# draw
#head(data.integrated)
#png(file=paste('n03.diff/',"Cell_UMAP_DimPlot3.png",sep="/"), width = dpi*10, height = dpi*8, units = "px",res = dpi,type='cairo')
#pp_temp=DimPlot(data.integrated,reduction='umap',label=TRUE,pt.size=0.5)+NoLegend() # umap 可视化
#pp_temp=DimPlot(data.integrated,reduction='umap',label=TRUE,repel = TRUE)
#pp_temp=DimPlot(data.integrated,reduction='umap',group.by = "ident")pp_temp=DimPlot(data.integrated,reduction='umap',group.by = "patient")
#pp_temp=DimPlot(data.integrated,reduction='umap',group.by = "seurat_clusters",label=TRUE,repel = TRUE)
head(data.integrated1)
png(file=paste('./',"Cell_split_Dimplot2.png",sep="/"), width = dpi*10, height = dpi*8, units = "px",res = dpi,type='cairo')
pp_temp = DimPlot(object = data.integrated1, reduction = 'umap',label = FALSE, label.size = 6,split.by = 'type', ncol = 3,repel =TRUE,combine = TRUE)
print(pp_temp)
dev.off()