lengshizai 2022-10-26 23:54 采纳率: 39.1%
浏览 856
已结题

使用monocle进行拟时序分析时,在细胞排序orderCells时报错怎么办

在用R语言中的monocle进行拟时序分析时,在画拟时许图的前一步细胞排序时报错,怎么办?

用代码块功能插入代码,请勿粘贴截图
##报错的代码是本篇代码最后一句,这是前面相关联的。前面的东西都没有报错。
##准备细胞轨迹分析所需要的文件
monocle.matrix=as.matrix(pbmc@assays$RNA@data)
monocle.sample=pbmc@meta.data
monocle.geneAnn=data.frame(gene_short_name = row.names(monocle.matrix), row.names = row.names(monocle.matrix))
monocle.clusterAnn=clusterAnn
monocle.markers=sig.markers

#?将Seurat结果转换为monocle需要的细胞矩阵,细胞注释表格和基因注释表格
data <- (as(as.matrix(monocle.matrix), 'sparseMatrix'))
pd<-new("AnnotatedDataFrame", data = monocle.sample)
fd<-new("AnnotatedDataFrame", data = monocle.geneAnn)
cds <- newCellDataSet(data, phenoData = pd, featureData = fd)
names(pData(cds))[names(pData(cds))=="seurat_clusters"]="Cluster"
pData(cds)[,"Cluster"]=paste0("cluster",pData(cds)[,"Cluster"])
#添加细胞聚类的数据
clusterAnn=as.character(monocle.clusterAnn[,2])
names(clusterAnn)=paste0("cluster",monocle.clusterAnn[,1])
pData(cds)$cell_type2 <- plyr::revalue(as.character(pData(cds)$Cluster),clusterAnn)
#细胞轨迹分析流程
cds <- estimateSizeFactors(cds)
cds <- estimateDispersions(cds)
cds <- setOrderingFilter(cds, as.vector(sig.markers$gene))
#plot_ordering_genes(cds)
cds <- reduceDimension(cds, max_components = 2, reduction_method = 'DDRTree')
cds <- orderCells(cds)##就是这句报错,报错提示在下面
运行结果及报错内容
Error in if (class(projection) != "matrix") projection <- as.matrix(projection) : 
  the condition has length > 1
In addition: Warning message:
In graph.dfs(dp_mst, root = root_cell, neimode = "all", unreachable = FALSE,  :
  Argument `neimode' is deprecated; use `mode' instead

我的解答思路和尝试过的方法

这个完全没头绪,非常高兴各位指点一下能学习的地方

  • 写回答

2条回答 默认 最新

  • Bioinfo Guy R语言领域新星创作者 2022-10-27 09:52
    关注

    报错矩阵构建有问题,试试改一下第一句

    monocle.matrix=as.sparse(pbmc@assays$RNA@data)
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(1条)

报告相同问题?

问题事件

  • 系统已结题 4月16日
  • 已采纳回答 4月8日
  • 修改了问题 10月26日
  • 修改了问题 10月26日
  • 展开全部

悬赏问题

  • ¥15 微信会员卡等级和折扣规则
  • ¥15 微信公众平台自制会员卡可以通过收款码收款码收款进行自动积分吗
  • ¥15 随身WiFi网络灯亮但是没有网络,如何解决?
  • ¥15 gdf格式的脑电数据如何处理matlab
  • ¥20 重新写的代码替换了之后运行hbuliderx就这样了
  • ¥100 监控抖音用户作品更新可以微信公众号提醒
  • ¥15 UE5 如何可以不渲染HDRIBackdrop背景
  • ¥70 2048小游戏毕设项目
  • ¥20 mysql架构,按照姓名分表
  • ¥15 MATLAB实现区间[a,b]上的Gauss-Legendre积分