请问各位大佬，“RegProcess('pisces-1-net-final.tsv', mat1, out.dir = 'tutorial/', out.name = 'pisces-1-net-')” 这里的mat1是不是 pisces1 <- readRDS('pisces-clust_aracne-mats/pisces-1-arac.rds')？。如果是，在运行过程中连续出现WARNING! NOT CONVERGENT! number of iterations= 1000 是哪里出了问题？
我顺便把LouvainKRange 也写一下，方便后来人复现。感谢楼主
LouvainKRange <- function(data.obj, kmin = 5, kmax = 50, kstep = 5) {
require(igraph, quietly = TRUE)
# check if seurat object
if (class(data.obj)[1] == "Seurat") {
dist.mat <- data.obj@assays[[data.obj@active.assay]]@misc$dist.mat
} else {
dist.mat <- data.obj
}
# create lists
cluster.list <- list()
sil.list <- list()
# setup iteration
k <- kmin
while (k <= kmax) {
print(paste("Clustering with k = ", k, "...", sep = ''))
# generate clustering and silhouette score
clust.vec <- LouvainClust(dist.mat, k)
sil.score <- cluster::silhouette(clust.vec, dist.mat)
# add to list
k.ind <- paste('k', k, sep = '.')
cluster.list[[k.ind]] <- clust.vec
sil.list[[k.ind]] <- mean(sil.score[,3])
# iterate
k <- k + kstep
}
# identify optimal cluster
opt.clust <- cluster.list[[which.max(sil.list)]]
# add to data.object
if (class(data.obj)[1] == "Seurat") {
data.obj@assays[[data.obj@active.assay]]@misc[['pisces.cluster']] <- opt.clust
data.obj@assays[[data.obj@active.assay]]@misc[['clustering.obj']] <- list('clusterings' = cluster.list, 'sils' = sil.list)
return(data.obj)
} else {
return(list('pisces.cluster' = opt.clust, 'clustering.obj' = list('clusterings' = cluster.list, 'sils' = sil.list)))
}
}
写在前面 感谢califano-lab[https://github.com/califano-lab]团队将代码无私的分享出来=。=！ 本流程需要一点点R语言与linux基...
最近帮师弟做一个GSVA分析，发现已经忘的差不多了，顺便整理一下吧。 基因集变异分析（Gene Set Variation Analysis, GSVA）可以将一个基因表达矩...
对SQL语句不熟悉欢迎查看我整理的笔记：[SQL] MySQL基础 + python交互转载请注明：陈熹 chenx6542@foxmail.com （简书号：半为花间酒）若...
生信技能树B站 GEO挖掘教学视频
使用SCENIC分析singlecell-RNAseq实战演练SCENIC全称：Single-Cell rEgulatory Network Inference and Clustering实战之前需要准备的学习资料：1.阅读关于SCE...
内容不是教程，均为个人学习笔记，以备后日查询学习资源均来自生信技能树论坛，生信技能树B站视频，微信公众号：生信技能树 为啥是我最不想用的方法？ 常用平台可通过biocondu...
SCENIC全称：Single-Cell rEgulatory Network Inference and Clustering实战之前需要准备的学习资料：1.阅读关于SCE...
内容不是教程，均为个人学习笔记，以备后日查询学习资源均来自生信技能树论坛，生信技能树B站视频，微信公众号：生信技能树 画热图我们需要给R两个数据：表达矩阵与样品分组信息 读入...
内容不是教程，均为个人学习笔记，以备后日查询学习资源均来自生信技能树论坛，生信技能树B站视频，微信公众号：生信技能树 需要用到的R包为：org.Hs.eg.db 首先看看包内...
