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

相关文章

Python FastAPI+Celery+RabbitMQ实现分布式图片水印处理系统

《PythonFastAPI+Celery+RabbitMQ实现分布式图片水印处理系统》这篇文章主要为大家详细介绍了PythonFastAPI如何结合Celery以及RabbitMQ实现简单的分布式... 实现思路FastAPI 服务器Celery 任务队列RabbitMQ 作为消息代理定时任务处理完整

C#使用SQLite进行大数据量高效处理的代码示例

《C#使用SQLite进行大数据量高效处理的代码示例》在软件开发中,高效处理大数据量是一个常见且具有挑战性的任务,SQLite因其零配置、嵌入式、跨平台的特性,成为许多开发者的首选数据库,本文将深入探... 目录前言准备工作数据实体核心技术批量插入:从乌龟到猎豹的蜕变分页查询:加载百万数据异步处理:拒绝界面

Springboot处理跨域的实现方式(附Demo)

《Springboot处理跨域的实现方式(附Demo)》:本文主要介绍Springboot处理跨域的实现方式(附Demo),具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不... 目录Springboot处理跨域的方式1. 基本知识2. @CrossOrigin3. 全局跨域设置4.

在java中如何将inputStream对象转换为File对象(不生成本地文件)

《在java中如何将inputStream对象转换为File对象(不生成本地文件)》:本文主要介绍在java中如何将inputStream对象转换为File对象(不生成本地文件),具有很好的参考价... 目录需求说明问题解决总结需求说明在后端中通过POI生成Excel文件流,将输出流(outputStre

python+opencv处理颜色之将目标颜色转换实例代码

《python+opencv处理颜色之将目标颜色转换实例代码》OpenCV是一个的跨平台计算机视觉库,可以运行在Linux、Windows和MacOS操作系统上,:本文主要介绍python+ope... 目录下面是代码+ 效果 + 解释转HSV: 关于颜色总是要转HSV的掩膜再标注总结 目标:将红色的部分滤

Python实现自动化接收与处理手机验证码

《Python实现自动化接收与处理手机验证码》在移动互联网时代,短信验证码已成为身份验证、账号注册等环节的重要安全手段,本文将介绍如何利用Python实现验证码的自动接收,识别与转发,需要的可以参考下... 目录引言一、准备工作1.1 硬件与软件需求1.2 环境配置二、核心功能实现2.1 短信监听与获取2.

Python使用date模块进行日期处理的终极指南

《Python使用date模块进行日期处理的终极指南》在处理与时间相关的数据时,Python的date模块是开发者最趁手的工具之一,本文将用通俗的语言,结合真实案例,带您掌握date模块的六大核心功能... 目录引言一、date模块的核心功能1.1 日期表示1.2 日期计算1.3 日期比较二、六大常用方法详

利用Go语言开发文件操作工具轻松处理所有文件

《利用Go语言开发文件操作工具轻松处理所有文件》在后端开发中,文件操作是一个非常常见但又容易出错的场景,本文小编要向大家介绍一个强大的Go语言文件操作工具库,它能帮你轻松处理各种文件操作场景... 目录为什么需要这个工具?核心功能详解1. 文件/目录存javascript在性检查2. 批量创建目录3. 文件

Java使用多线程处理未知任务数的方案介绍

《Java使用多线程处理未知任务数的方案介绍》这篇文章主要为大家详细介绍了Java如何使用多线程实现处理未知任务数,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 知道任务个数,你可以定义好线程数规则,生成线程数去跑代码说明:1.虚拟线程池:使用 Executors.newVir

一文带你深入了解Python中的GeneratorExit异常处理

《一文带你深入了解Python中的GeneratorExit异常处理》GeneratorExit是Python内置的异常,当生成器或协程被强制关闭时,Python解释器会向其发送这个异常,下面我们来看... 目录GeneratorExit:协程世界的死亡通知书什么是GeneratorExit实际中的问题案例