温馨提示×

温馨提示×

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

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

R语言的univariate_cox.r怎么使用

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

这篇文章主要介绍了R语言的univariate_cox.r怎么使用的相关知识,内容详细易懂,操作简单快捷,具有一定借鉴价值,相信大家阅读完这篇R语言的univariate_cox.r怎么使用文章都会有所收获,下面我们一起来看看吧。

单因素cox回归分析

使用方法:

指定单个变量做cox回归分析,可以是临床数据

$Rscript $scriptdir/univariate_cox.r -h
usage: /share/nas1/huangls/test/TCGA_immu/scripts/univariate_cox.r
       [-h] -i data -t time -e event -v variate [variate ...] [-b blocksize]
       [-o outdir] [-p prefix]

unvariate cox regression

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]
  -b blocksize, --blocksize blocksize
                        Number of variables Parallel to test in each [default
                        2]
  -o outdir, --outdir outdir
                        output file directory [default cwd]
  -p prefix, --prefix prefix
                        out file name prefix [default cox]

使用举例:

Rscript $scriptdir/univariate_cox.r -i nomogram_metadata.tsv -e EVENT -t TIME \
  -v gender  riskScore age_at_index    ajcc_pathologic_stage \
  -o univariate_cox \
  -p  univariate_cox

参数说明:

-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











结果展示:

VariableTermBetaHRHRlowerHRupperP
riskScoreriskScore1.3124753.7153562.3299775.9244683.53E-08
ajcc_pathologic_stageajcc_pathologic_stageStage IV1.4190894.1333531.9438288.7891540.000227
ajcc_pathologic_stageajcc_pathologic_stageStage III0.8551992.3518431.2037484.5949530.012327
age_at_indexage_at_index0.0181241.0182891.0007321.0361540.04111
ajcc_pathologic_stageajcc_pathologic_stageStage II0.5336181.705090.8381573.4687170.140829
gendergendermale0.2178161.2433580.8579621.8018740.249868







关于“R语言的univariate_cox.r怎么使用”这篇文章的内容就介绍到这里,感谢各位的阅读!相信大家对“R语言的univariate_cox.r怎么使用”知识都有一定的了解,大家如果还想学习更多知识,欢迎关注亿速云行业资讯频道。

向AI问一下细节

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

AI