单细胞测序FindAllmarkers

执行此代码后Rscript $scripts/seurat_FindAllMarkers.r --rds subset.Epithelial.cells.0.2.rds \

>  -p subset.Epithelial.cells.0.2.FindAllMarkers --test.use  DESeq2   --logfc.threshold 0.25,报错,FindAllMarkers

Calculating cluster 0

converting counts to integer mode

Calculating cluster 1

converting counts to integer mode

Calculating cluster 2

converting counts to integer mode

Calculating cluster 3

converting counts to integer mode

Calculating cluster 4

converting counts to integer mode

Calculating cluster 5

converting counts to integer mode

Calculating cluster 6

converting counts to integer mode

Calculating cluster 7

converting counts to integer mode

Calculating cluster 8

converting counts to integer mode

Warning: No DE genes identified

Warning: The following tests were not performed: 

Warning: When testing 0 versus all:

        every gene contains at least one zero, cannot compute log geometric means

Warning: When testing 1 versus all:

        every gene contains at least one zero, cannot compute log geometric means

Warning: When testing 2 versus all:

        every gene contains at least one zero, cannot compute log geometric means

Warning: When testing 3 versus all:

        every gene contains at least one zero, cannot compute log geometric means

Warning: When testing 4 versus all:

        every gene contains at least one zero, cannot compute log geometric means

Warning: When testing 5 versus all:

        every gene contains at least one zero, cannot compute log geometric means

Warning: When testing 6 versus all:

        every gene contains at least one zero, cannot compute log geometric means

Warning: When testing 7 versus all:

        every gene contains at least one zero, cannot compute log geometric means

Warning: When testing 8 versus all:

        every gene contains at least one zero, cannot compute log geometric means

data frame with 0 columns and 0 rows

Error in `group_by_prepare()`:

! Must group by variables found in `.data`.

• Column `cluster` is not found.

Backtrace:

    ▆

 1. ├─obj.markers %>% group_by(cluster) %>% ...

 2. ├─dplyr::top_n(., n = opt$top_n, wt = avg_log2FC)

 3. │ └─dplyr::filter(...)

 4. ├─dplyr::group_by(., cluster)

 5. └─dplyr:::group_by.data.frame(., cluster)

 6.   └─dplyr::group_by_prepare(.data, ..., .add = .add, caller_env = caller_env())

 7.     └─rlang::abort(c("Must group by variables found in `.data`.", glue("Column `{unknown}` is not found.")))

Execution halted

想请假老师,我应该怎么调整?谢谢老师
请先 登录 后评论

2 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

没有找到差异基因,降低标准试试的: --logfc.threshold 0.15

请先 登录 后评论


老师,我降低了依然是原来的错误。我把之前可以正常运行的文件重新执行代码,也是报这样的错误?

请先 登录 后评论
  • 2 关注
  • 0 收藏,136 浏览
  •  提出于 2025-03-23 14:12