这篇文章主要介绍“R语言的immune_infiltrates_rnaseq.r RNA-seq怎么用”的相关知识,小编通过实际案例向大家展示操作过程,操作方法简单快捷,实用性强,希望这篇“R语言的immune_infiltrates_rnaseq.r RNA-seq怎么用”文章能帮助大家解决问题。
immune_infiltrates_rnaseq.r RNA-seq 转录组数据免疫侵润分析
$ Rscript immune_infiltrates_rnaseq.r -h usage: immune_infiltrates_rnaseq.r [-h] -i expset [-g gene.info] [-t type] [--tpm] [-o outdir] RNA seq 免疫侵润分析. In general values should be TPM-normalized ,not log- transformed. optional arguments: -h, --help show this help message and exit -i expset, --expset expset input gene expression set matrix from RNA-seq data tsv format [required] -g gene.info, --gene.info gene.info input gene info data [required] -t type, --type type TIMER uses indication-specific reference profiles. [optional] --tpm whether convert fpkm to tpm [optional, default: False] -o outdir, --outdir outdir output file directory [default cwd]
-i 输入基因表达量文件,建议用TPM标准化之后的数据,如果是FPKM可以设置--tpm进行转换;第一列ID为基因NAME
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 |
-t 指定癌症类型, TIMER计算需要,可以指定类型,如果不指定 不输出TIMER的分析结果:
# "kich" "blca" "brca" "cesc" "gbm" "hnsc" "kirp" "lgg" "lihc" "luad"
# "lusc" "prad" "sarc" "pcpg" "paad" "tgct" "ucec" "ov" "skcm" "dlbc"
# "kirc" "acc" "meso" "thca" "uvm" "ucs" "thym" "esca" "stad" "read"
# "coad" "chol"
使用R包:immunedeconv
可以输出以下方法的免疫侵润分析结果,mcp_counter 需要联网才可以运行,有时报错不出结果;
quantiseq timer cibersort cibersort_abs mcp_counter xcell epic
关于“R语言的immune_infiltrates_rnaseq.r RNA-seq怎么用”的内容就介绍到这里了,感谢大家的阅读。如果想了解更多行业相关的知识,可以关注亿速云行业资讯频道,小编每天都会为大家更新不同的知识点。
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。