利用R读取h5ad文件

安装包

install.packages(“remotes”)
remotes::install_github(“mojaveazure/seurat-disk”,force = TRUE)

加载包
library(SeuratDisk)
library(patchwork)
library(Seurat)
library(dplyr)
library(ggplot2)

转换

Convert(‘TS_Eye.h5ad’,”h5seurat”,
overwrite=TRUE,assay=”RNA”) #TS_Eye.h5ad为需要转换的文件
seurat_obj<-LoadH5Seurat(“TS_Eye.h5seurat”,verbose=T)

列重命名

colnames(seurat_obj@meta.data)[7]<-“nCount_RNA”
colnames(seurat_obj@meta.data)[8]<-“nFeature_RNA”
colnames(seurat_obj@meta.data)[9]<-“percent.mt”

本站原创,如若转载,请注明出处:https://www.ouq.net/1800.html

(0)
打赏 微信打赏,为服务器增加50M流量 微信打赏,为服务器增加50M流量 支付宝打赏,为服务器增加50M流量 支付宝打赏,为服务器增加50M流量
上一篇 05/01/2022 20:26
下一篇 05/16/2022 21:32

相关推荐

  • Rstudio/Rstudio Server enable Copilot-Rstudio Server打开Copilot

    Rstudio Server 默认是关闭Copilot Copilot is turned off by default. Copilot is turned off with copilot-enabled=0 in /etc/rstud…

    R 06/05/2025
    110
  • R_Code: KEGG analysis

    library(clusterProfiler) library(org.Hs.eg.db) # 读取输入数据文件 file_path <- “C:/Users/Lamarck/Desktop/UP_genes_ENSEMBL_ENT…

    R 06/02/2025
    133
  • R_Code:GO and Functional GO

    GO analysis: library(AnnotationDbi) library(org.Hs.eg.db) #基因注释包 library(clusterProfiler) #富集包 # 读取CSV文件 file_path <-…

    R 06/02/2025
    135
  • R_Code:WGCNA and WGCNA_Get_Gene_Length

    library(WGCNA) library(DESeq2) # enableWGCNAThreads(nThreads = 10) # 在处理数据框(data.frame)时,不会自动给将String类型转换成factor类型 optio…

    R 06/02/2025
    138
  • R_Code:RNAseq_GSEA_analysis

    library(clusterProfiler) # GSEA 和富集分析主力包 library(org.Hs.eg.db) # 人类注释数据库(ENTREZID 与 SYMBOL 等 ID 转换) library(enrichplot) …

    R 06/02/2025
    137