【无标题】

rm(list = ls())
library(dplyr)
library(Seurat)
library(patchwork)
library(ggplot2)
setwd("D:/脑谱图/newdata")
data <- readRDS("D:/脑谱图/newdata/data.rds")
data <- readRDS("data.rds")
data <- UpdateSeuratObject(object = data)####seurat版本不够

 

install.packages("SeuratObject")

Selection: 16
Warning message:
package 'SeuratObject' is not available (for R version 3.6.3)

#使用read.table读取meta文件
meta <- read.table("D:/脑谱图/newdata/meta.tsv", header=T, sep="\t")
# Convert between row names and column --------------------------------
colnames(meta)##使用colnames获取列名
meta <- meta %>% column_to_rownames(var = "cellId") 
colnames(meta)

 思路：使用UCSC下载标准10X文件再添加meta信息再提取表达的假基因

######从UCSC下载标准10X文件再添加meta信息
bra <- Read10X(data.dir = "D:/脑谱图/CS12_Hindbrain.mex/CS12_Hindbrain")
bra <- CreateSeuratObject(counts = bra,Project="", min.cells = 0, min.features = 0)
#使用read.table读取meta文件
meta <- read.table("D:/脑谱图/newdata/meta.tsv", header=T, sep="\t")
# Convert between row names and column --------------------------------
colnames(meta)##使用colnames获取列名
meta <- meta %>% column_to_rownames(var = "cellId") 
colnames(meta)
####在bra的基础上，使用AddMetaData函数进行添加
bra <- AddMetaData(
  object = bra,
  metadata = meta
)

 2.提取表达的假基因