温馨提示×

温馨提示×

您好,登录后才能下订单哦!

密码登录×
登录注册×
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》

R语言的multi_cox.r怎么用

发布时间:2022-03-21 10:35:25 来源:亿速云 阅读:129 作者:iii 栏目:开发技术

本篇内容主要讲解“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 输入生存数据与基因表达文件 

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











到此,相信大家对“R语言的multi_cox.r怎么用”有了更深的了解,不妨来实际操作一番吧!这里是亿速云网站,更多相关内容可以进入相关频道进行查询,关注我们,继续学习!

向AI问一下细节

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

AI