您好,登錄后才能下訂單哦!
這篇文章主要介紹“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進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。