千家信息网

R语言的univariate_cox_batch.r怎么用

发表于:2026-11-03 作者:千家信息网编辑
千家信息网最后更新 2026年11月03日,本文小编为大家详细介绍"R语言的univariate_cox_batch.r怎么用",内容详细,步骤清晰,细节处理妥当,希望这篇"R语言的univariate_cox_batch.r怎么用"文章能帮助
千家信息网最后更新 2026年11月03日R语言的univariate_cox_batch.r怎么用

本文小编为大家详细介绍"R语言的univariate_cox_batch.r怎么用",内容详细,步骤清晰,细节处理妥当,希望这篇"R语言的univariate_cox_batch.r怎么用"文章能帮助大家解决疑惑,下面跟着小编的思路慢慢深入,一起来学习新知识吧。

univariate_cox_batch.r 基因表达量批量单因素cox回归分析

使用说明:

输入生存数据与基因表达量可以做批量单因素cox回归分析

$ Rscript $scriptdir/univariate_cox_batch.r --husage:  univariate_cox_batch.r [-h] -m metadata -g                                                         expset [-t time]                                                         [-e event]                                                         [-l pvalue]                                                         [-b blocksize]                                                         [--log2] [-o outdir]                                                         [-p prefix]batch unvariate cox regression gene expressionoptional arguments:  -h, --help            show this help message and exit  -m metadata, --metadata metadata                        input metadata file path with suvival time [required]  -g expset, --expset expset                        input gene expression set file [required]  -t time, --time time  set suvival time column name in metadata [default                        TIME]  -e event, --event event                        set event column name in metadata [default EVENT]  -l pvalue, --pvalue pvalue                        pvalue cutoff to choose sig gene [default 0.01]  -b blocksize, --blocksize blocksize                        Number of variables Parallel to test in each [default                        2]  --log2                whether do log2 transfrom for expression data                        [optional, default: False]  -o outdir, --outdir outdir                        output file directory [default cwd]  -p prefix, --prefix prefix                        out file name prefix [default cox]

参数说明:

-m 输入生存数据:

event 列: 0表示事件没有发生,1表示事件发生; 0表示alive,1表示死亡;

barcodeTIMEEVENT
TCGA-B7-A5TK-01A-12R-A36D-312880
TCGA-BR-7959-01A-11R-2343-1310100
TCGA-IN-8462-01A-11R-2343-135720
TCGA-CG-4443-01A-01R-1157-139120
TCGA-KB-A93J-01A-11R-A39E-3111240
TCGA-HU-A4H3-01A-21R-A251-318820
TCGA-RD-A8MV-01A-11R-A36D-3137200
TCGA-VQ-A91X-01A-12R-A414-312891
TCGA-D7-8575-01A-11R-2343-135541
TCGA-BR-8485-01A-11R-2402-132800
TCGA-D7-A748-01A-12R-A32D-311321
TCGA-VQ-A91Z-01A-11R-A414-3116900

-g 输入基因表达量文件

IDTCGA-B7-A5TK-01A-12R-A36D-31TCGA-BR-7959-01A-11R-2343-13TCGA-IN-8462-01A-11R-2343-13TCGA-BR-A4CR-01A-11R-A24K-31TCGA-CG-4443-01A-01R-1157-13TCGA-KB-A93J-01A-11R-A39E-31TCGA-BR-4371-01A-01R-1157-13TCGA-IN-A6RO-01A-12R-A33Y-31TCGA-HU-A4H3-01A-21R-A251-31
FGR16.3440811.967395.3508462.2093511.5380215.240164.5011182.6024376.261761
CD3886.8677215.794513.1113421.2407070.86295513.30473.7287081.6739522.675173
ITGAL40.269037.3585663.7691252.3878692.3735138.085918.3052833.6227817.025886
CX3CL1603.013226.9135320.222384.19526219.0409714.1529513.758856.6753744.050271
CEACAM211.8685362.5719170.6108390.6745581.0921273.4835591.1343094.4712740.584159
MATK2.283420.8641160.5197762.4420930.7603483.1929511.1618810.3478821.039336
CD79B3.4531981.8799572.8221920.5235871.9265923.6517420.8312880.8836431.979214
MMP2513.728293.4511481.1065631.1312170.87873510.431861.4758521.9142842.312993
TRAF3IP35.244011.8801860.8752640.7561530.6032513.3250132.3474730.5704621.315916
CD477.7469151.8371922.7707611.0781135.20445122.557831.1010715.0661915.41347
BTK6.8562354.3622611.4826881.3715991.9812366.911543.1878480.9554991.48269
FMO17.1685677.7118173.2231740.9790340.4503071.0934121.0018080.9102041.558515
SYT71.15310581.940682.673384191.611282.493940.5103734.4704821.285060.91944
TYROBP591.7796338.0271184.813369.18483150.6397480.5691121.09672.4588116.9793
CD220.8192952.5216071.5885050.412590.3872881.1236330.4882440.2580940.713988


使用举例:

Rscript $scriptdir/univariate_cox_batch.r -m metadata_survival_time.tsv \   -g deg_gene_exp_tpm.tsv  -e EVENT -t TIME -p  imm.unicox  --pvalue 0.01

批量cox分析 结果:

VariableTermBetaStandardErrorZPLRTWaldLogRankHRHRlowerHRupper
SYT12SYT120.0911210.0194954.6740352.95E-060.0001282.95E-061.80E-061.0954021.0543361.138067
CDH2CDH20.0132660.0030144.4018031.07E-050.0019931.07E-052.11E-071.0133541.0073861.019357
GPNMBGPNMB0.0027590.0007683.5901180.0003310.0013130.0003310.0004091.0027621.0012531.004274
TMIGD3TMIGD30.067880.0193143.5144680.0004410.0012480.0004410.0004191.0702371.030481.111528
LINC01094LINC010940.1324410.0403753.2802930.0010370.0022420.0010370.0010261.1416111.0547541.235621
SLC22A20PSLC22A20P0.0494150.015833.121650.0017980.0120650.0017980.0007361.0506561.0185591.083765
IGHV4-61IGHV4-610.0017910.0005823.0775730.0020870.008590.0020870.0017421.0017931.0006511.002937
IGHV2-5IGHV2-50.0022360.0007373.0349610.0024060.0077820.0024060.0022051.0022381.0007921.003686
SERPINA5SERPINA50.0076810.0025583.0024280.0026780.0097990.0026780.0020641.0077111.002671.012776
MS4A4AMS4A4A0.014460.005052.8632180.0041940.0084270.0041940.0047221.0145651.0045721.024657
FAM83AFAM83A0.0062370.0023682.6334450.0084520.0285230.0084520.0072951.0062561.0015961.010938
IGLV3-9IGLV3-90.0005470.000212.6082810.00910.0394270.00910.0102441.0005471.0001361.000958
STARD3STARD30.0009280.0003582.5889130.0096280.0302580.0096280.0067431.0009281.0002251.001631

读到这里,这篇"R语言的univariate_cox_batch.r怎么用"文章已经介绍完毕,想要掌握这篇文章的知识点还需要大家自己动手实践使用过才能领会,如果想了解更多相关内容的文章,欢迎关注行业资讯频道。

语言 基因 文章 分析 输入 事件 内容 因素 数据 妥当 使用说明 参数 思路 文件 新知 更多 步骤 知识 知识点 篇文章 数据库的安全要保护哪些东西 数据库安全各自的含义是什么 生产安全数据库录入 数据库的安全性及管理 数据库安全策略包含哪些 海淀数据库安全审计系统 建立农村房屋安全信息数据库 易用的数据库客户端支持安全管理 连接数据库失败ssl安全错误 数据库的锁怎样保障安全 互联网科技企业 互联网科技名片 百单网互联网科技有限公司怎么样 工业互联网科技新基建 小马过河互联网科技 互联网科技上市公司 中国最先进互联网科技公司 中国的互联网科技企业的崛起 宝灵互联网科技有限公司 互联网与医疗科技 呼和浩特新华互联网科技学院 互联网金融科技展 智帮互联网科技公司最近动向 广西科技产业互联网是什么 互联网科技新闻 云南互联网科技学费 北京军民互联网科技公司 互联网企业和网络科技有限公司 高科技和互联网的现代国家 小牛互联网科技有限公司地址 工业互联网科技新基建 有马互联网科技 互联网 新科技 互联网科技创新的路径 上海科技产业互联网平台 富国互联网金融科技混合基金 互联网科技股为什么大涨 卡尔迅互联网科技有限公司 吉林省码上互联网科技电话 银川新华互联网科技学校原画
0