admin管理员组

文章数量:1650776

在运行此代码时报错，请问如何解决？感谢！

expr_matrix <- as(as.matrix(seurat_object_cluster@assays$RNA@counts), 'sparseMatrix')
f_data <- data.frame(gene_short_name = row.names(seurat_object_cluster),row.names = row.names(seurat_object_cluster))
pd <- new('AnnotatedDataFrame', data =seurat_object_cluster@meta.data) 
fd <- new('AnnotatedDataFrame', data = f_data)
#将p_data和f_data从data.frame转换AnnotatedDataFrame对象。
cds <- newCellDataSet(expr_matrix,
                      phenoData = pd,
                      featureData = fd,
                      lowerDetectionLimit = 0.5,
                      expressionFamily = negbinomial.size())

我发现在seurat_object_cluster中的counts是33231 x 45183，而f_data提取出来的是2000，是否是这里出现的问题，求教如何解决？

 

 

本文标签： 报错featurenewCellDataSetassayDatafeatureData

版权声明：本文标题：求助！使用monocle2 创建newCellDataSet时，报错feature numbers differ between assayData and featureData 内容由热心网友自发贡献，该文观点仅代表作者本人， 转载请联系作者并注明出处：https://www.elefans.com/xitong/1729531010a1204804.html， 本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容，一经查实，本站将立刻删除。