$Rscript ../scripts/exp_immu_cor_plot.r -h
usage: ../scripts/exp_immu_cor_plot.r [-h] -e expset -i immu [-m method]
[-a alternative] [-n ncol] [-l]
[-x label.x.npc] [-y label.y.npc]
[-o outdir] [-H height] [-W width]
scattor correlateion plot: https://www.omicsclass.com/article/1488
optional arguments:
-h, --help show this help message and exit
-e expset, --expset expset
set event column name [required]
-i immu, --immu immu input data file path[required]
-m method, --method method
a character string indicating which correlation
coefficient (or covariance) is to be computed. One of
"pearson", "kendall", or "spearman". [default pearson]
-a alternative, --alternative alternative
a character string specifying the alternative
hypothesis, must be one of "two.sided" , "greater" or
"less". You can specify just the initial letter.
[default two.sided]
-n ncol, --ncol ncol column num [default [default 4]
-l, --log2 Whether to do log2 transformation of gene expression
[default FALSE]
-x label.x.npc, --label.x.npc label.x.npc
value should be between 0 and 1. Coordinates to be
used for positioning the label, expressed in
"normalized parent coordinates". [default 0]
-y label.y.npc, --label.y.npc label.y.npc
value should be between 0 and 1. Coordinates to be
used for positioning the label, expressed in
"normalized parent coordinates". [default 1]
-o outdir, --outdir outdir
output file directory [default cwd]
-H height, --height height
the height of pic inches [default 5]
-W width, --width width
the width of pic inches [default 10]
Rscript ../scripts/exp_immu_cor_plot.r -e gene_exp1.tsv -i ../03.TIME/immu/timer.res.tsv -o timer_cor1 --log2
参数说明:
-e 指定输入指定基因的表达量,行为 基因,列为样本名称
| ID | TCGA-D7-A74A-01A-11R-A32D-31 | TCGA-BR-7704-01A-11R-2055-13 | TCGA-VQ-A91N-01A-11R-A414-31 | TCGA-CD-A4MH-01A-11R-A251-31 |
| NUP50 | 18.65505 | 31.59232 | 28.23382 | 28.76485 |
| CXCR4 | 64.85805 | 125.123 | 56.35244 | 69.98976 |
| NT5E | 111.4818 | 69.8587 | 79.37382 | 25.05824 |
| EFNA3 | 8.247857 | 42.03308 | 43.46432 | 26.66024 |
| STC1 | 4.781111 | 21.36327 | 40.81077 | 19.51568 |
| ZBTB7A | 95.51678 | 103.4768 | 158.3024 | 126.2677 |
| CLDN9 | 1.187456 | 2.476138 | 0.366081 | 7.347344 |
-i 输入样本免疫侵润结果, 注意脚本会自动对两个文件里的样本取交集
| cell_type | TCGA-B7-A5TK-01A-12R-A36D-31 | TCGA-BR-7959-01A-11R-2343-13 | TCGA-IN-8462-01A-11R-2343-13 | TCGA-BR-A4CR-01A-11R-A24K-31 |
| B cell | 0 | 0.013159742 | 0.000907507 | 0.020234 |
| T cell CD4+ | 0 | 0.114882797 | 0.118061151 | 0.152274 |
| T cell CD8+ | 0.771121767 | 0.130798349 | 0.048964712 | 0.005392 |
| Neutrophil | 0.22692792 | 0.125194736 | 0.106562343 | 0.057013 |
| Macrophage | 0 | 0.213365358 | 0 | 0.05444 |
| Myeloid dendritic cell | 0.935644191 | 0.457844399 | 0.404774813 | 0.296719 |

如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!
