您好,登錄后才能下訂單哦!
本篇內容主要講解“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怎么用”有了更深的了解,不妨來實際操作一番吧!這里是億速云網站,更多相關內容可以進入相關頻道進行查詢,關注我們,繼續學習!
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。