第一步：從GEO數(shù)據(jù)庫下載.h5文件
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM3892578
第二步:Rstudio加載包窗悯，讀入數(shù)據(jù)
? ??????library(Seurat)
????????library(ggplot2)
????????library(cowplot)
????????library(Matrix)
????????library(dplyr)
????????library(ggsci)
????????library(scater)
????????library(SingleR)
????????library(hdf5r)
###用read10x_h5函數(shù)讀取數(shù)據(jù)
? ??????Read10X_h5(GSM3892578_SKIN_HV1_F2_filtered_gene_bc_matrices_h5.h5)
Error in Read10X_h5(GSM3892578_SKIN_HV1_F2_filtered_gene_bc_matrices_h5.h5) :? Please install hdf5r to read HDF5 files ###顯示要安裝包
????????install.packages('hdf5r')
????????library(hdf5r)
????????h1= Read10X_h5(GSM3892578_SKIN_HV1_F2_filtered_gene_bc_matrices_h5.h5)
? ? ? ? #再用creatseuratobject函數(shù)創(chuàng)建seurat對象就可以了
? ??????HC_1 <- CreateSeuratObject(counts = hc1, project = "HC_1", min.cells = 3, min.features = 200)
HC1就是一個seurat對象啦，再開啟下游分析。
