溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊(cè)×
其他方式登錄
點(diǎn)擊 登錄注冊(cè) 即表示同意《億速云用戶服務(wù)條款》

R語(yǔ)言的multi_cox.r怎么用

發(fā)布時(shí)間:2022-03-21 10:35:25 來(lái)源:億速云 閱讀:121 作者:iii 欄目:開(kāi)發(fā)技術(shù)

本篇內(nèi)容主要講解“R語(yǔ)言的multi_cox.r怎么用”,感興趣的朋友不妨來(lái)看看。本文介紹的方法操作簡(jiǎn)單快捷,實(shí)用性強(qiáng)。下面就讓小編來(lái)帶大家學(xué)習(xí)“R語(yǔ)言的multi_cox.r怎么用”吧!

multi_cox.r 多因素cox分析

使用方法:

指定多個(gè)基因表達(dá)量做多因素cox回歸分析,并構(gòu)建預(yù)后模型,以及評(píng)估預(yù)后模型:

$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

參數(shù)說(shuō)明:

-i 輸入生存數(shù)據(jù)與基因表達(dá)文件 

barcodeTIMEEVENTFGRCD38ITGALCX3CL1CEACAM21MATKCD79BMMP25
TCGA-B7-A5TK-01A-12R-A36D-31288016.3440886.8677240.26903603.01321.8685362.283423.45319813.72829
TCGA-BR-7959-01A-11R-2343-131010011.9673915.794517.35856626.913532.5719170.8641161.8799573.451148
TCGA-IN-8462-01A-11R-2343-1357205.3508463.1113423.76912520.222380.6108390.5197762.8221921.106563
TCGA-CG-4443-01A-01R-1157-1391201.538020.8629552.3735119.040971.0921270.7603481.9265920.878735
TCGA-KB-A93J-01A-11R-A39E-311124015.2401613.304738.0859114.152953.4835593.1929513.65174210.43186
TCGA-HU-A4H3-01A-21R-A251-3188206.2617612.6751737.0258864.0502710.5841591.0393361.9792142.312993
TCGA-RD-A8MV-01A-11R-A36D-313720027.0741520.1588534.9130934.718214.1131122.61555716.5194617.72674
TCGA-VQ-A91X-01A-12R-A414-3128911.0623410.7520182.3805134.4158150.5181420.2121971.2392030.582114











到此,相信大家對(duì)“R語(yǔ)言的multi_cox.r怎么用”有了更深的了解,不妨來(lái)實(shí)際操作一番吧!這里是億速云網(wǎng)站,更多相關(guān)內(nèi)容可以進(jìn)入相關(guān)頻道進(jìn)行查詢,關(guān)注我們,繼續(xù)學(xué)習(xí)!

向AI問(wèn)一下細(xì)節(jié)

免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如果涉及侵權(quán)請(qǐng)聯(lián)系站長(zhǎng)郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。

AI