1. Seurat 自带方法 Idents(Seurat_object) <- "Column_name_in_metadata" 2. 自己写的方法
1. Seurat 自带方法 Idents(Seurat_object) <- "Column_name_in_metadata" 2. 自己写的方法
公众号后台回复20210418,获取本文的测试数据 1. 回顾 前面我已经讲解了inferCNV的基本用法[https://www.jianshu.com/p/c01f56b...
挖掘公共单细胞数据集时,会遇到常见各种单细胞测序数据格式。现总结如下,方便自己日后调用,以创建Seurat对象(1)barcodes.tsv.gz、features.tsv....
SeuratData SeuratData是一种使用R的内部软件包和数据管理系统以Seurat对象形式分发数据集的机制。它为用户提供了一种访问Seurat小插曲中使用的数据集...
1. ##从rfam下载所有的序列## wget ftp://ftp.ebi.ac.uk/pub/databases/Rfam/12.0/fasta_files/* cat ...
Rfam是用来鉴定非编码RNA的数据库,常用于注释新的核酸序列或基因组序列。 最新版本的Rfam14.4,已经包含了3941个家族。 Rfam官方网址:https://rfa...
R语言接受命令行参数的三种方式 最近要用一下R的命令行,但是发现命令行参数不是很友好,就总结了一下目前常用的R的命令行参数的接受方法。 平常我们使用R语言的时候常常都是在交互...
在单细胞数据分析中,在确定细胞类型后,除了可以进行差异表达基因分析外,还可以针对单个细胞类型进行分析特定分析,这时就需要我们提取细胞子集分开处理了。 一、Seurat数据格式...
本文涉及:亚群合并,亚群提取,亚群细分,亚群抽样 一般来讲,不同亚群的细胞具有不同的功能,因此亚群分析对于研究十分重要,会影响后续分析结果。 进行细胞亚群的分群的依据: 细胞...
GNU Parallel是一个shell工具,为了在一台或多台计算机上并行的执行计算任务。本文简要介绍GNU Parallel的使用。 1、parallel 用法简介 1.1...
今天将为大家介绍小麦的特有的生物学网站,该数据库是由国际小麦基因组测序联盟(International Wheat Genome Sequencing Consortium,...
注意,这里使用的seurat对象要求已经run过runUMAMP()findCluster等函数,否则也没有必要把seurat的结果弄到monocle3的cds对象里 1. ...
scRNA矩阵存储的文件格式有10X单细胞测序数据、h5、h5ad、loom:10X单细胞测序数据经过cellranger处理后会得到矩阵的三个文件:matrix.mtx...
常用的细胞通讯软件:CellphoneDB[https://www.jianshu.com/p/38a9376f5286]:是公开的人工校正的,储存受体、配体以及两种相互作用...
KEGG pathway 注释整理 获得KEGG注释 通过eggnog-mapper和interproscan两个软件(或数据库),可以获得KEGG ORTHOLOGY(KO...
查看perl 版本 查看perl模块搜索路径 查看某个perl模块是否存在 模块安装方法1 —— 命令行形式安装Set::IntervalTree 模块安装方法2 —— co...
以Pacific Biosciences (PacBio)公司为代表的第三代单分子实时荧光测序平台为以测序为基础的科学研究带来了革命性的突破。目前该技术广泛应用于基因组Den...
目前做GSEA主要有两种办法,一种是使用GSEA的java软件,另一种是使用R语言包,当时还有在线的工具做GSEA的数据也有两种,一种是有基因表达矩阵和分组,另一种是使用两列...
相关系数主要分为Pearson、Spearman相关系数 0.8-1.0 极强相关0.6-0.8 强相关0.4-0.6 中等程度相关0.2-0.4 ...
R中的线性回归函数比较简单,就是lm(),比较复杂的是对线性模型的诊断和调整。这里结合Statistical Learning和杜克大学的Data Analysis and ...