seurat对象处理 找锚点

2023-10-08 07:40
文章标签 对象 处理 锚点 seurat

本文主要是介绍seurat对象处理 找锚点,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

在找锚点合并之前,需要把每个seurat对象的细胞名改变成唯一
getwd() #改名字

#教程地址
#https://cloud.tencent.com/developer/article/1697249#https://bioconductor.org/packages/release/data/experiment/vignettes/scRNAseq/inst/doc/scRNAseq.html
#https://mp.weixin.qq.com/mp/appmsgalbum?__biz=MzI1Njk4ODE0MQ==&action=getalbum&album_id=1326587538303434752&scene=173&from_msgid=2247484689&from_itemidx=1&count=3&nolastread=1#wechat_redirect
rm(list = ls()) 
Sys.setenv(R_MAX_NUM_DLLS=999)
options(stringsAsFactors = F)

##########三对三的数据
1#准备原始分析数据  先手动下载  去浏览器下载文件到自己的文件夹下,然后解压
#https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi  
# 3个样本,共9个文件。需要分到9个文件夹里,并且重命名
getwd()
path="G:/silicosis/geo/GSE128033_SnRNAseq-idiopathic pulmonary fibrosis/GSE128033" #空间转录组
dir.create(path)
setwd(path)
getwd()
fs=list.files('./','^GSM')  #得到当前目录下所有以GSM开头的文件名称
fs

在这里插入图片描述

# 自行下载GSE164621数据集的GSE164621_RAW压缩包并且解压哦,这样上面的代码就可以运行啦# 然后获取3个样本信息,因为是批量,所以下面的代码可能不好理解,需要熟练掌握R语言哦
library(stringr)
samples=str_split(fs,'_',simplify = T)[,1] #取出可以分组的样本名 有三个重复GSM5015042  GSM5015043 GSM5015044 得到样本分组
samples

在这里插入图片描述

getwd()
###注意workshop的位置,决定下一个语句是否成功!
setwd("G:/silicosis/geo/GSE128033_SnRNAseq-idiopathic pulmonary fibrosis/")lapply(unique(samples),function(x){#x=unique(samples)[1]y=fs[grepl(x,fs)]folder=paste0("GSE128033/", str_split(y[1],'_',simplify = T)[,1])dir.create(folder,recursive = T)#为每个样本创建子文件夹file.rename(paste0("GSE128033/",y[1]),file.path(folder,"barcodes.tsv.gz"))#重命名文件,并移动到相应的子文件夹里file.rename(paste0("GSE128033/",y[2]),file.path(folder,"features.tsv.gz"))file.rename(paste0("GSE128033/",y[3]),file.path(folder,"matrix.mtx.gz"))
})
getwd()
samples=list.files("GSE128033/")
samples# 是两个文件夹的名字哦

在这里插入图片描述

samples=list.files("GSE128033/") %>% grep(pattern = "gz",invert = TRUE,value = TRUE)
samples# 是两个文件夹的名字哦2#创建Seurat对library(Seurat)#测试用
if(1==1){myexam=  CreateSeuratObject(counts = Read10X(folder),project = "pro" ,min.cells = 3,min.features=200)myexam=subset(myexam)
}
dim(myexam)
getwd()
# 循环读取两个文件夹下面的10x的的3个文件
sceList = lapply(samples,function(pro){#pro="GSM3660641"folder=file.path("GSE128033/",pro)CreateSeuratObject(counts = Read10X(folder),  #此处一定要加上筛选条件project = pro ,min.cells=3,min.features=200)##print(paste0("进行到第",pro,"个样本")) 不可以打印这一排
})

添加print回发生错误,所有把print语句注释掉
在这里插入图片描述
在这里插入图片描述

for (i in 1:length(sceList)) {#计算线粒体比例sceList[[i]][["percent.mt"]]=PercentageFeatureSet(sceList[[i]],pattern = "^MT-")
}
sceList
DefaultAssay(sceList)getwd()
save(sceList,file = 

这篇关于seurat对象处理 找锚点的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/163761

相关文章

Java对象转换的实现方式汇总

《Java对象转换的实现方式汇总》:本文主要介绍Java对象转换的多种实现方式,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录Java对象转换的多种实现方式1. 手动映射(Manual Mapping)2. Builder模式3. 工具类辅助映

Python Transformers库(NLP处理库)案例代码讲解

《PythonTransformers库(NLP处理库)案例代码讲解》本文介绍transformers库的全面讲解,包含基础知识、高级用法、案例代码及学习路径,内容经过组织,适合不同阶段的学习者,对... 目录一、基础知识1. Transformers 库简介2. 安装与环境配置3. 快速上手示例二、核心模

一文详解Java异常处理你都了解哪些知识

《一文详解Java异常处理你都了解哪些知识》:本文主要介绍Java异常处理的相关资料,包括异常的分类、捕获和处理异常的语法、常见的异常类型以及自定义异常的实现,文中通过代码介绍的非常详细,需要的朋... 目录前言一、什么是异常二、异常的分类2.1 受检异常2.2 非受检异常三、异常处理的语法3.1 try-

Python使用getopt处理命令行参数示例解析(最佳实践)

《Python使用getopt处理命令行参数示例解析(最佳实践)》getopt模块是Python标准库中一个简单但强大的命令行参数处理工具,它特别适合那些需要快速实现基本命令行参数解析的场景,或者需要... 目录为什么需要处理命令行参数?getopt模块基础实际应用示例与其他参数处理方式的比较常见问http

Java Response返回值的最佳处理方案

《JavaResponse返回值的最佳处理方案》在开发Web应用程序时,我们经常需要通过HTTP请求从服务器获取响应数据,这些数据可以是JSON、XML、甚至是文件,本篇文章将详细解析Java中处理... 目录摘要概述核心问题:关键技术点:源码解析示例 1:使用HttpURLConnection获取Resp

Java中Switch Case多个条件处理方法举例

《Java中SwitchCase多个条件处理方法举例》Java中switch语句用于根据变量值执行不同代码块,适用于多个条件的处理,:本文主要介绍Java中SwitchCase多个条件处理的相... 目录前言基本语法处理多个条件示例1:合并相同代码的多个case示例2:通过字符串合并多个case进阶用法使用

Java实现优雅日期处理的方案详解

《Java实现优雅日期处理的方案详解》在我们的日常工作中,需要经常处理各种格式,各种类似的的日期或者时间,下面我们就来看看如何使用java处理这样的日期问题吧,感兴趣的小伙伴可以跟随小编一起学习一下... 目录前言一、日期的坑1.1 日期格式化陷阱1.2 时区转换二、优雅方案的进阶之路2.1 线程安全重构2

Python处理函数调用超时的四种方法

《Python处理函数调用超时的四种方法》在实际开发过程中,我们可能会遇到一些场景,需要对函数的执行时间进行限制,例如,当一个函数执行时间过长时,可能会导致程序卡顿、资源占用过高,因此,在某些情况下,... 目录前言func-timeout1. 安装 func-timeout2. 基本用法自定义进程subp

Java字符串处理全解析(String、StringBuilder与StringBuffer)

《Java字符串处理全解析(String、StringBuilder与StringBuffer)》:本文主要介绍Java字符串处理全解析(String、StringBuilder与StringBu... 目录Java字符串处理全解析:String、StringBuilder与StringBuffer一、St

Python中判断对象是否为空的方法

《Python中判断对象是否为空的方法》在Python开发中,判断对象是否为“空”是高频操作,但看似简单的需求却暗藏玄机,从None到空容器,从零值到自定义对象的“假值”状态,不同场景下的“空”需要精... 目录一、python中的“空”值体系二、精准判定方法对比三、常见误区解析四、进阶处理技巧五、性能优化