单细胞聚类方法SC3疯狂报错——Error in [.default(clusts, , paste0(“sc3_“, hash.table[i, 1], “_clusters“))

SC3 - consensus clustering of single-cell RNA-Seq data这是scRNA-seq聚类方法里面比较经典的聚类方法,在跑数据集的时候,几个一千多个细胞的小数据集没有报错,但是运行一个5000+细胞以上的数据集时,在这步:

sce <- sc3(sce, ks = 9 ,n_cores = 30)

疯狂报错!!!
在这里插入图片描述
就很无语,难道这种经典方法我都跑不通嘛!到了推理原因的时刻了
First,小数据集上我能跑通,所以问题不是在这个方法上,应该是因为数据集大小的改变,哪些参数不合适了,second, 我去看了SC3的手册:https://bioconductor.org/packages/release/bioc/manuals/SC3/man/SC3.pdf
在这里插入图片描述
发现kmeans_nstart这个参数以2000个数据为分界,小于2000,参数设置为1000,大于2000设置为50,finally, 将代码改成:

sce <- sc3(sce, ks = 9 ,n_cores = 30,kmeans_nstart = 50)

即可解决问题!
遇到问题,推理出问题原因哇!

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值