sct:表示是否使用了SCTransform方法进行标准化。 scRNA <- doubletFinder_v3(scRNA, PCs = 1:15, pN = 0.25, pK = pK_bcmvn, nExp = nExp_poi.adj, reuse.pANN = F, sct = T) ## 运行结果储存在scRNA@meta.data中 #可视化如图 #将双细胞剔除后生成新的对象scRNA.singlet,便于我们后续分析 scRNA....
1.数据的预处理 代码语言:javascript 复制 library(Seurat)library(DoubletFinder)library(hdf5r)pbmc<-Read10X_h5("5k_pbmc_v3_nextgem_filtered_feature_bc_matrix.h5")pbmc<-CreateSeuratObject(pbmc)pbmc<-SCTransform(pbmc)pbmc<-RunPCA(pbmc,verbose=F)pc.num=1:30pbmc<-RunUMAP(pbmc,dims=pc.num)pbmc<...
nExp_poi<-round(DoubletRate*nrow(scRNA@meta.data))# 使用同源双细胞比例对推算出期望的异源双细胞数量nExp_poi.adj <- round(nExp_poi*(1-homotypic.prop))##使用DoubletFinder算法来检测和去除单细胞数据中的异源双细胞。doubletFinder_v3函数会对单细胞对象scRNA进行双细胞检测和鉴定。这个函数有以下参数:...
#开始鉴定双细胞 pbmc <- doubletFinder_v3(pbmc, PCs= 1:10, pN = 0.25, pK = pK_bcmvn, nExp = nExp_poi.adj, reuse.pANN = F, sct = F) 主要参数如下所述: http://events.jianshu.io/p/6770c6a05287 从下面的结果,可以看出预测完之后,亏有pANN和DF.classifications两个量。 ###下面,我...
(假设双细胞形成率为 7.5%)nExp_poi<-round(0.075*nrow(scRNA_harmony@meta.data))nExp_poi.adj<-round(nExp_poi*(1-homotypic.prop))# 计算异源双细胞数量# 使用确定好的参数鉴定doubletsscRNA_harmony<-doubletFinder_v3(scRNA_harmony,PCs=1:20,pN=0.25,pK=pK_bcmvn,nExp=nExp_poi.adj,reuse.pANN...
软件包安装很简单,github模式安装即可: ## 载入需要的程辑包:DoubletFinder if(!require(DoubletFinder)) remotes::install_github('chris-mcginnis-ucsf/DoubletFinder') 数据读取 这里我们使用 pbmc3k 的数据集来做例子,完成整个分析过程,这里我们需要调用 Seurat 软件包来完成数据的处理工作。数据读入如下: ...
软件安装 代码语言:javascript 复制 remotes::install_github('chris-mcginnis-ucsf/DoubletFinder') 10X PBMC数据实测 测试数据来源于10x genomics官网的示例数据集。 代码语言:javascript 复制 数据下载链接:https://cf.10xgenomics.com/samples/cell-exp/3.0.2/5k_pbmc_v3_nextgem/5k_pbmc_v3_nextgem_filtered_...
软件包安装很简单,github模式安装即可: ## 载入需要的程辑包:DoubletFinder if(!require(DoubletFinder)) remotes::install_github('chris-mcginnis-ucsf/DoubletFinder') 数据读取 这里我们使用 pbmc3k 的数据集来做例子,完成整个分析过程,这里我们需要调用 Seurat 软件包来完成数据的处理工作。数据读入如下: ...
软件安装 代码语言:javascript 复制 remotes::install_github('chris-mcginnis-ucsf/DoubletFinder') 10X PBMC数据实测 测试数据来源于10x genomics官网的示例数据集。 代码语言:javascript 复制 数据下载链接:https://cf.10xgenomics.com/samples/cell-exp/3.0.2/5k_pbmc_v3_nextgem/5k_pbmc_v3_nextgem_filtered_...
安装: devtools::install_github('chris-mcginnis-ucsf/DoubletFinder') 2.1 寻找最优pK值 library(DoubletFinder)#这是一个测试最佳参数的过程,运行速度慢sweep.res.list<-paramSweep_v3(scRNA,PCs=1:pcSelect,sct=F)#使用log标准化,sct参数设置为 sct = F(默认 ),如使用SCT标准化方法,设置为Tsweep.stats...