没有找到差异基因,降低标准试试的: --logfc.threshold 0.15
执行此代码后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