RCTD算法
用途
使用RCTD做解卷积细胞注释
运行方式
SDAS cellAnnotation rctd -i st.h5ad -o outdir --reference sc.h5ad --bin_size 20 --label_key annotation2 \
--filter_rare_cell 0
输入参数说明
| 参数 | 是否必须 | 默认值 | 描述 |
|---|---|---|---|
| -i / --input | 是 | Stereo-seq h5ad，要求有原始表达矩阵 | |
| -o / --output | 是 | 输出文件夹 | |
| --reference | 是 | 单细胞ref h5ad，要求有原始表达矩阵 | |
| --label_key | 是 | 单细胞ref h5ad.obs中表示细胞类型的列的名称 | |
| --bin_size | 是 | Bin大小，用于控制图中点的大小，不用于计算,比如20,50,100, cellbin (等效于20) | |
| --input_layer | 否 | Stereo-seq h5ad存放raw counts的layer | |
| --ref_layer | 否 | 单细胞ref h5ad存放raw counts的layer | |
| --input_gene_symbol_key | 否 | real_gene_name | Stereo-seq h5ad.var中表示基因名(symbol)的列的名称 (_index 表示使用h5ad.var.index) |
| --ref_gene_symbol_key | 否 | _index | 单细胞ref h5ad.var中表示基因名(symbol)的列的名称 (_index 表示使用h5ad.var.index) |
| --slice_key | 否 | sampleID | 多片h5ad.obs中表示片编号的列的名称，用于画图 |
| --filter_rare_cell | 否 | 100 | 如果某些细胞类型在单细胞ref中细胞数小于此值，则过滤掉这些细胞类型 |
| --mode | 否 | full | RCTD模式。选项：doublet, multi, full |
| --UMI_min | 否 | 100 | 小于此UMI数的空转spot会被过滤 |
| --counts_MIN | 否 | 10 | 小于此基因数的空转spot会被过滤 |
| --seed | 否 | 42 | 随机种子设置 |
| --n_cpus | 否 | 8 | 使用的进程数 |
输出结果展示
| 结果文件 | 描述 |
|---|---|
<input_name>_anno_rctd.csv |
每个spot的注释结果，包括每种细胞类型的分数 |
<input_name>_anno_rctd.h5ad |
输入h5ad+注释结果。每个细胞类型的分数存在obsm['anno_score_rctd']中，分数最高的类型存在obs['anno_rctd']中 |
<input_name>_anno_rctd.png/pdf |
总体注释结果图，多片情况下每片画一张图，同时输出png和pdf |
<input_name>_anno_rctd_split.png/pdf |
每个细胞类型分开展示图，多片情况下每片画一张图，同时输出png和pdf |
<input_name>_anno_score_rctd.png/pdf |
每个细胞类型的分数图，多片情况下每片画一张图，同时输出png和pdf |
详细说明与具体结果展示可参考以下链接。(cell2location算法-->细胞注释-->输出结果展示)。