$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 | 

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