首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往
您找到你想要的搜索结果了吗?
是的
没有找到

使用pyscenic做转录因子分析

时间消耗对比 一般来说pyscenic是在Linux相关服务器上面操作,所以我也是在自己的ubuntu操作系统演示: 安装自己的conda,每个用户独立操作 安装方法代码如下: # 首先下载文件,20M...然后,有意思的地方就来了,我测试了五台服务器,发现并不是所有的服务器都可以成功安装它!后面我们会讲解一个docker的解决方案,假如你确实只有的conda无法安装,也可以等明天的教程哈!...提取单细胞表达量矩阵csv忘记并且导入Linux服务器 首先我们对文章《Single-cell RNA sequencing highlights the role of inflammatory cancer-associated...<- rownames(fibo_sce) fibo_sce <- ScaleData(fibo_sce, features = all.genes) #去除MT,重新进行标准化 fibo_sce <...(object = fibo_sce), verbose = FALSE) #聚类 fibo_sce <- FindNeighbors(fibo_sce, dims = 1:20) fibo_sce <

11K32

人人都需掌握的单细胞分群聚类分析&Marker基因的可视化

) View(sce$orig.ident) sce$orig.ident<-annotation$Patient #修改sce 矩阵样本的分组信息 View(sce$orig.ident) #识别表达量高变的基因...(sce) sce = ScaleData(sce, vars.to.regress = c("nFeature_RNA", "percent_mito")) sce = RunPCA(sce, npcs...sce = RunTSNE(sce, npcs = 20) sce = Runtsne(sce, dims = 1:10) sce = RunUMAP(sce, dims = 1:10) sce = FindNeighbors...强烈建议你推荐给身边的博士后以及年轻生物学PI,多一点数据认知,让他们的科研上一个台阶: 生物信息学马拉松授课(买一得五) ,你的生物信息学入门课 时隔5年,我们的生信技能树VIP学徒继续招生啦 144线程640Gb内存服务器共享一年仍然是仅需...800 千呼万唤始出来的独享生物信息学云服务器 生信技能树知识整理实习生又又又开放申请啦 生信共享办公室出租

2.4K20

OSCA单细胞数据分析笔记-3 SingleCellExperiment数据结构

四组成 3.1 标准化数据之assays 相关函数命令 assays(sce) 查看当前sce对象的所有assays' name【暂且可理解一种表达矩阵称之为一个assay】 assay(sce,"name...sce=logNormCounts(sce) assays(sce) # List of length 2 # names(2): counts logcounts assay(sce,"logcounts...尤其特殊一点是colData(sce)foo等价于sce #当前sce有一条sizeFactor cell-info,为scater::logNormCounts时自动添加的 colData(sce)...利用scater包进行PCA降维 sce <- runPCA(sce) reducedDims(sce) # List of length 1 # names(1): PCA reducedDim(sce...---- 如果你对单细胞转录组研究感兴趣,但又不知道如何入门,也许你可以关注一下下面的课程 数据挖掘(GEO,TCGA,单细胞)2021第2期 生信爆款入门-2021第2期 96核心384G内存的超级服务

3.7K54

小鼠的5个样品的10x技术单细胞转录组上游定量(文末赠送全套代码)

这在你需要在远程服务器上运行长时间的任务或同时运行多个任务时特别有用。 以下是一些基本的screen命令: screen:启动一个新的screen会话。...使用screen命令可以让你在单个终端窗口中更有效地管理你的工作,特别是在处理远程服务器或长时间运行的任务时。...使用初探 单细胞实战(四) Cell Ranger流程概览 单细胞实战(五) 理解cellranger count的结果 理论上应该是去10x的官网下载cellranger软件和数据库文件,不过我们的共享服务器已经给大家准备了...cellranger软件和数据库文件,以及前面的kingfisher下载了PRJNA853539里面的小鼠的5个样品的10x技术单细胞转录组上游测序数据,接下来就可以跑cellranger啦,因为是共享服务器...虽然说5个样品一起运行会很快,但是大家的CPU和内存可能会不够,如果是共享服务器或者云服务器,另外一个门槛是硬盘。。。。 有可能是硬盘也不够哦。。。。

50920

猪的单细胞分析如何过滤线粒体基因

举个例子: mito_genes=rownames(sce.all)[grep("^MT-", rownames(sce.all))] mito_genes #13个线粒体基因 # character...文件,希望在里面能找到些信息,解压并打开文件,果不其然,文件中的gene命名是以ENSSSCG开头的; 找一个基因测试下ENSSSCG00000018060 > mito_genes=rownames(sce.all...这个代码来自于我们《生信技能树》的单细胞数据分析工程师 林枫,欢迎大家前来咨询我们团队的单细胞数据分析服务哦 接下来我们进行实战 前面的PercentageFeatureSet()函数有两种指定模式:...)) kp=mt.genes %in% rownames(sce.all) table(kp) mt.genes=mt.genes[kp] C<-GetAssayData(object = sce.all...(sce.all, percent.mito, col.name = "percent.mito") sce.all[["percent.mito"]] fivenum(sce.all@meta.data

1.9K20

跟学单细胞周更(一)

acc=GSE175687 获取数据: 这里给大家介绍一下我个人下载时发现的一个小技巧 我想要下载这个数据集下四个矩阵文件 复制打包好的tar文件的http链接 发现是一个跳转链接,无法直接服务器中下载...(如果点击下载下到自己电脑上,再传给服务器,有点麻烦) 检索资料,发现有人说使用curl下载 -L参数可以跟随重定向进行下载 但经curl -L -O测试使用,发现下载到的文件和前面直接wget是一样的...<- RunPCA(sce, features = VariableFeatures(object = sce))可以发现,并不是拿所有的基因进行PCA,而是前面识别高变基因找到的VariableFeatures...1) } colnames(sce.all@meta.data) apply(sce.all@meta.data[,grep("RNA_snn",colnames(sce.all@meta.data)...=sce[,sce$celltype=='epi'] #sce.epi=sce[,sce$RNA_snn_res.0.8 %in% c(1,3,4,7,18,22,24,25,26,32,34,35)]

20620

什么,需要5Tb内存啊!

最近我就接到一个粉丝咨询,说他想处理一个公共数据集,只有8个原位肿瘤+3个转移肿瘤的10X单细胞转录组样品,但是数据处理的过程发现系统提示说需要5Tb内存,虽然说他自己有一个512G内存的服务器,但是也承受不起...额,给他配置一个5Tb内存服务器倒是简单,我自己就有2.5T内存的服务器,不就是加倍嘛!不过,我注意到他就是11个10X转录组样品,理论上不可能是需要5Tb内存的,所以让他把代码发过来我检查看看....<- FindIntegrationAnchors(object.list = sceList, dims = 1:30) # load('sce.anchors.RDATA') Sys.time()...sce.integrated <- IntegrateData(anchorset = sce.anchors, dims = 1:30) Sys.time() 看起来中规中矩,就是我一直教学使用的代码...见:共享云服务器又又又又又来了(只有50个名额)

1K30

跟学单细胞周更(一)

acc=GSE175687 获取数据: 这里给大家介绍一下我个人下载时发现的一个小技巧 我想要下载这个数据集下四个矩阵文件 复制打包好的tar文件的http链接 发现是一个跳转链接,无法直接服务器中下载...(如果点击下载下到自己电脑上,再传给服务器,有点麻烦) 检索资料,发现有人说使用curl下载 -L参数可以跟随重定向进行下载 但经curl -L -O测试使用,发现下载到的文件和前面直接wget是一样的...<- RunPCA(sce, features = VariableFeatures(object = sce))可以发现,并不是拿所有的基因进行PCA,而是前面识别高变基因找到的VariableFeatures...1) } colnames(sce.all@meta.data) apply(sce.all@meta.data[,grep("RNA_snn",colnames(sce.all@meta.data)...=sce[,sce$celltype=='epi'] #sce.epi=sce[,sce$RNA_snn_res.0.8 %in% c(1,3,4,7,18,22,24,25,26,32,34,35)]

17330

给你一个病毒感染与否的时间序列采样的转录组你能挖掘出什么

$celltype) sce.all=sce.all.int[,sce.all.int$celltype=='mac'] av <-AverageExpression(sce.all ,...$orig.ident), function(x){ #x = unique(sce.all$orig.ident)[1] kp = sce.all$orig.ident %in%...x print(table(kp)) rowSums( as.matrix( sce.all@assays$RNA@counts[, kp] )) }) ) colnames(ct...强烈建议你推荐给身边的博士后以及年轻生物学PI,多一点数据认知,让他们的科研上一个台阶: 生物信息学马拉松授课(买一得五) ,你的生物信息学入门课 时隔5年,我们的生信技能树VIP学徒继续招生啦 144线程640Gb内存服务器共享一年仍然是仅需...800 千呼万唤始出来的独享生物信息学云服务器 生信技能树知识整理实习生又又又开放申请啦(不招了,谢谢) 生信共享办公室出租

23530
领券