本篇内容主要讲解“R语言的multi_cox.r怎么用”，感兴趣的朋友不妨来看看。本文介绍的方法操作简单快捷，实用性强。下面就让小编来带大家学习“R语言的multi_cox.r怎么用”吧!
multi_cox.r 多因素cox分析
使用方法：
指定多个基因表达量做多因素cox回归分析，并构建预后模型，以及评估预后模型：
$Rscript $scriptdir/multi_cox.r -h usage: /share/nas1/huangls/test/TCGA_immu/scripts/multi_cox.r [-h] -i data -t time -e event -v variate [variate ...] [-P predict.time [predict.time ...]] [-c cut.score] [-s seed] [-o outdir] [-p prefix] multi variate cox regression analysis using gene expression optional arguments: -h, --help show this help message and exit -i data, --data data input data file path[required] -t time, --time time set suvival time column name [required] -e event, --event event set event column name must 0 or 1 code format [required] -v variate [variate ...], --variate variate [variate ...] variate for cox analysis [required] -P predict.time [predict.time ...], --predict.time predict.time [predict.time ...] Time point to draw the ROC curve [default 365 1095 1825] -c cut.score, --cut.score cut.score set cut score value to divide high and low risk groups [default median] -s seed, --seed seed set random seed [default 2021] -o outdir, --outdir outdir output file directory [default cwd] -p prefix, --prefix prefix out file name prefix [default cox]
使用举例：
Rscript $scriptdir/multi_cox.r -i imm.unicox.metadata-exp.tsv -e EVENT -t TIME \ -v PDGFRL CXCR4 PAK3 CSF1R PDCD1 -P 365 1095 1825 \ -o multicox -p multicox
参数说明：
-i 输入生存数据与基因表达文件
| barcode | TIME | EVENT | FGR | CD38 | ITGAL | CX3CL1 | CEACAM21 | MATK | CD79B | MMP25 |
| TCGA-B7-A5TK-01A-12R-A36D-31 | 288 | 0 | 16.34408 | 86.86772 | 40.26903 | 603.0132 | 1.868536 | 2.28342 | 3.453198 | 13.72829 |
| TCGA-BR-7959-01A-11R-2343-13 | 1010 | 0 | 11.96739 | 15.79451 | 7.358566 | 26.91353 | 2.571917 | 0.864116 | 1.879957 | 3.451148 |
| TCGA-IN-8462-01A-11R-2343-13 | 572 | 0 | 5.350846 | 3.111342 | 3.769125 | 20.22238 | 0.610839 | 0.519776 | 2.822192 | 1.106563 |
| TCGA-CG-4443-01A-01R-1157-13 | 912 | 0 | 1.53802 | 0.862955 | 2.37351 | 19.04097 | 1.092127 | 0.760348 | 1.926592 | 0.878735 |
| TCGA-KB-A93J-01A-11R-A39E-31 | 1124 | 0 | 15.24016 | 13.3047 | 38.08591 | 14.15295 | 3.483559 | 3.192951 | 3.651742 | 10.43186 |
| TCGA-HU-A4H3-01A-21R-A251-31 | 882 | 0 | 6.261761 | 2.675173 | 7.025886 | 4.050271 | 0.584159 | 1.039336 | 1.979214 | 2.312993 |
| TCGA-RD-A8MV-01A-11R-A36D-31 | 3720 | 0 | 27.07415 | 20.15885 | 34.91309 | 34.71821 | 4.113112 | 2.615557 | 16.51946 | 17.72674 |
| TCGA-VQ-A91X-01A-12R-A414-31 | 289 | 1 | 1.062341 | 0.752018 | 2.380513 | 4.415815 | 0.518142 | 0.212197 | 1.239203 | 0.582114 |
到此，相信大家对“R语言的multi_cox.r怎么用”有了更深的了解，不妨来实际操作一番吧！这里是云搜网网站，更多相关内容可以进入相关频道进行查询，关注我们，继续学习！
