R语言的univariate_cox_batch.r怎么用
发表于:2025-01-21 作者:千家信息网编辑
千家信息网最后更新 2025年01月21日,本文小编为大家详细介绍"R语言的univariate_cox_batch.r怎么用",内容详细,步骤清晰,细节处理妥当,希望这篇"R语言的univariate_cox_batch.r怎么用"文章能帮助
千家信息网最后更新 2025年01月21日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表示死亡;
barcode | TIME | EVENT |
TCGA-B7-A5TK-01A-12R-A36D-31 | 288 | 0 |
TCGA-BR-7959-01A-11R-2343-13 | 1010 | 0 |
TCGA-IN-8462-01A-11R-2343-13 | 572 | 0 |
TCGA-CG-4443-01A-01R-1157-13 | 912 | 0 |
TCGA-KB-A93J-01A-11R-A39E-31 | 1124 | 0 |
TCGA-HU-A4H3-01A-21R-A251-31 | 882 | 0 |
TCGA-RD-A8MV-01A-11R-A36D-31 | 3720 | 0 |
TCGA-VQ-A91X-01A-12R-A414-31 | 289 | 1 |
TCGA-D7-8575-01A-11R-2343-13 | 554 | 1 |
TCGA-BR-8485-01A-11R-2402-13 | 280 | 0 |
TCGA-D7-A748-01A-12R-A32D-31 | 132 | 1 |
TCGA-VQ-A91Z-01A-11R-A414-31 | 1690 | 0 |
-g 输入基因表达量文件
ID | TCGA-B7-A5TK-01A-12R-A36D-31 | TCGA-BR-7959-01A-11R-2343-13 | TCGA-IN-8462-01A-11R-2343-13 | TCGA-BR-A4CR-01A-11R-A24K-31 | TCGA-CG-4443-01A-01R-1157-13 | TCGA-KB-A93J-01A-11R-A39E-31 | TCGA-BR-4371-01A-01R-1157-13 | TCGA-IN-A6RO-01A-12R-A33Y-31 | TCGA-HU-A4H3-01A-21R-A251-31 |
FGR | 16.34408 | 11.96739 | 5.350846 | 2.209351 | 1.53802 | 15.24016 | 4.501118 | 2.602437 | 6.261761 |
CD38 | 86.86772 | 15.79451 | 3.111342 | 1.240707 | 0.862955 | 13.3047 | 3.728708 | 1.673952 | 2.675173 |
ITGAL | 40.26903 | 7.358566 | 3.769125 | 2.387869 | 2.37351 | 38.08591 | 8.305283 | 3.622781 | 7.025886 |
CX3CL1 | 603.0132 | 26.91353 | 20.22238 | 4.195262 | 19.04097 | 14.15295 | 13.75885 | 6.675374 | 4.050271 |
CEACAM21 | 1.868536 | 2.571917 | 0.610839 | 0.674558 | 1.092127 | 3.483559 | 1.134309 | 4.471274 | 0.584159 |
MATK | 2.28342 | 0.864116 | 0.519776 | 2.442093 | 0.760348 | 3.192951 | 1.161881 | 0.347882 | 1.039336 |
CD79B | 3.453198 | 1.879957 | 2.822192 | 0.523587 | 1.926592 | 3.651742 | 0.831288 | 0.883643 | 1.979214 |
MMP25 | 13.72829 | 3.451148 | 1.106563 | 1.131217 | 0.878735 | 10.43186 | 1.475852 | 1.914284 | 2.312993 |
TRAF3IP3 | 5.24401 | 1.880186 | 0.875264 | 0.756153 | 0.603251 | 3.325013 | 2.347473 | 0.570462 | 1.315916 |
CD4 | 77.74691 | 51.83719 | 22.77076 | 11.07811 | 35.20445 | 122.5578 | 31.10107 | 15.06619 | 15.41347 |
BTK | 6.856235 | 4.362261 | 1.482688 | 1.371599 | 1.981236 | 6.91154 | 3.187848 | 0.955499 | 1.48269 |
FMO1 | 7.168567 | 7.711817 | 3.223174 | 0.979034 | 0.450307 | 1.093412 | 1.001808 | 0.910204 | 1.558515 |
SYT7 | 1.153105 | 81.94068 | 2.673384 | 191.6112 | 82.49394 | 0.510373 | 4.470482 | 1.28506 | 0.91944 |
TYROBP | 591.7796 | 338.0271 | 184.8133 | 69.18483 | 150.6397 | 480.5691 | 121.096 | 72.4588 | 116.9793 |
CD22 | 0.819295 | 2.521607 | 1.588505 | 0.41259 | 0.387288 | 1.123633 | 0.488244 | 0.258094 | 0.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分析 结果:
Variable | Term | Beta | StandardError | Z | P | LRT | Wald | LogRank | HR | HRlower | HRupper |
SYT12 | SYT12 | 0.091121 | 0.019495 | 4.674035 | 2.95E-06 | 0.000128 | 2.95E-06 | 1.80E-06 | 1.095402 | 1.054336 | 1.138067 |
CDH2 | CDH2 | 0.013266 | 0.003014 | 4.401803 | 1.07E-05 | 0.001993 | 1.07E-05 | 2.11E-07 | 1.013354 | 1.007386 | 1.019357 |
GPNMB | GPNMB | 0.002759 | 0.000768 | 3.590118 | 0.000331 | 0.001313 | 0.000331 | 0.000409 | 1.002762 | 1.001253 | 1.004274 |
TMIGD3 | TMIGD3 | 0.06788 | 0.019314 | 3.514468 | 0.000441 | 0.001248 | 0.000441 | 0.000419 | 1.070237 | 1.03048 | 1.111528 |
LINC01094 | LINC01094 | 0.132441 | 0.040375 | 3.280293 | 0.001037 | 0.002242 | 0.001037 | 0.001026 | 1.141611 | 1.054754 | 1.235621 |
SLC22A20P | SLC22A20P | 0.049415 | 0.01583 | 3.12165 | 0.001798 | 0.012065 | 0.001798 | 0.000736 | 1.050656 | 1.018559 | 1.083765 |
IGHV4-61 | IGHV4-61 | 0.001791 | 0.000582 | 3.077573 | 0.002087 | 0.00859 | 0.002087 | 0.001742 | 1.001793 | 1.000651 | 1.002937 |
IGHV2-5 | IGHV2-5 | 0.002236 | 0.000737 | 3.034961 | 0.002406 | 0.007782 | 0.002406 | 0.002205 | 1.002238 | 1.000792 | 1.003686 |
SERPINA5 | SERPINA5 | 0.007681 | 0.002558 | 3.002428 | 0.002678 | 0.009799 | 0.002678 | 0.002064 | 1.007711 | 1.00267 | 1.012776 |
MS4A4A | MS4A4A | 0.01446 | 0.00505 | 2.863218 | 0.004194 | 0.008427 | 0.004194 | 0.004722 | 1.014565 | 1.004572 | 1.024657 |
FAM83A | FAM83A | 0.006237 | 0.002368 | 2.633445 | 0.008452 | 0.028523 | 0.008452 | 0.007295 | 1.006256 | 1.001596 | 1.010938 |
IGLV3-9 | IGLV3-9 | 0.000547 | 0.00021 | 2.608281 | 0.0091 | 0.039427 | 0.0091 | 0.010244 | 1.000547 | 1.000136 | 1.000958 |
STARD3 | STARD3 | 0.000928 | 0.000358 | 2.588913 | 0.009628 | 0.030258 | 0.009628 | 0.006743 | 1.000928 | 1.000225 | 1.001631 |
读到这里,这篇"R语言的univariate_cox_batch.r怎么用"文章已经介绍完毕,想要掌握这篇文章的知识点还需要大家自己动手实践使用过才能领会,如果想了解更多相关内容的文章,欢迎关注行业资讯频道。
语言
基因
文章
分析
输入
事件
内容
因素
数据
妥当
使用说明
参数
思路
文件
新知
更多
步骤
知识
知识点
篇文章
数据库的安全要保护哪些东西
数据库安全各自的含义是什么
生产安全数据库录入
数据库的安全性及管理
数据库安全策略包含哪些
海淀数据库安全审计系统
建立农村房屋安全信息数据库
易用的数据库客户端支持安全管理
连接数据库失败ssl安全错误
数据库的锁怎样保障安全
安全服务器分类
广州多迪网络技术有限公司
中小学网络安全答题
稳定的二手服务器回收
矿机服务器
网络安全知识儿歌大全
数字化应用备案数据库收录信息
商标软件开发和运行是第几类
谈一下对网络安全网络文明的认识
软件开发各阶段文档标识
网络安全领导小组专题会议
数据库模型图导出
魔力宝贝 服务器架设
宝卡服务器怎么绑定域名
vb 远程连接数据库代码
网络安全工作方面的汇报材料
读取文件放入数据库
网络技术与电子信息工程的简介
倩女幽魂新服务器
一个软件开发公司多少人
安全服务器分类
双线服务器
海量数据库四大技术
服务器管理ppt素材
课程表下载软件开发
歌莉娅软件开发
服务器开机连续四声短鸣报警
数据库mc升级中啥意思
服务器管理员怎么输代码
一个软件开发公司多少人