美文网首页
R Seurat subset 报错 Errors "No ce

R Seurat subset 报错 Errors "No ce

作者: Aji | 来源:发表于2022-10-20 15:13 被阅读0次

参考: https://github.com/satijalab/seurat/issues/1511

好久没更新了,进来发现呜呜有一个关注和一个点赞,实在太随意的更新,还能有人喜欢和点赞,感动!!

1.问题

今天重新找到一篇文章的code, 准备copy它的代码进行分析,发现跑到subset那个代码时,发现报错了
哎想不通啊
然后记得以前遇到过这个问题,然后这次又碰到了,实在给忘了当时怎么解决的,然后就百度呀,发现真得也有同样的人遇到这个问题,我就记录下,给自己方便,也给他人方便。

2. 解决办法

github上提供的解决办法

其实就是active.ident不是factor变量,且这个变量的行名不是细胞名字,所以就没法subset啦
成功解决啦

  temp@active.ident <- factor(temp@active.ident)
  names(temp@active.ident) <- colnames(temp)
  print(dim(temp))
  temp <- subset(temp,  subset = nCount_RNA < 50000 &  nFeature_RNA < 6500 & nFeature_RNA > 500  & percent.rp < 5 & percent.mt < 5 )

写在后面啦
今天买的小林漫画的2023年日历到了,好开心呀,以下封面文字也送给看到这篇博文的你们啦
希望时光与爱 幸运常在!!


小林漫画日历

相关文章

网友评论

      本文标题:R Seurat subset 报错 Errors "No ce

      本文链接:https://www.haomeiwen.com/subject/bclmzrtx.html