91超碰碰碰碰久久久久久综合_超碰av人澡人澡人澡人澡人掠_国产黄大片在线观看画质优化_txt小说免费全本

溫馨提示×

溫馨提示×

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

密碼登錄×
登錄注冊×
其他方式登錄
點擊 登錄注冊 即表示同意《億速云用戶服務條款》

R語言的univariate_cox_batch.r怎么用

發布時間:2022-03-21 10:51:14 來源:億速云 閱讀:188 作者:iii 欄目:開發技術

本文小編為大家詳細介紹“R語言的univariate_cox_batch.r怎么用”,內容詳細,步驟清晰,細節處理妥當,希望這篇“R語言的univariate_cox_batch.r怎么用”文章能幫助大家解決疑惑,下面跟著小編的思路慢慢深入,一起來學習新知識吧。

univariate_cox_batch.r  基因表達量批量單因素cox回歸分析

使用說明:

輸入生存數據與基因表達量可以做批量單因素cox回歸分析

$ Rscript $scriptdir/univariate_cox_batch.r --h
usage:  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 expression
optional 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怎么用”文章已經介紹完畢,想要掌握這篇文章的知識點還需要大家自己動手實踐使用過才能領會,如果想了解更多相關內容的文章,歡迎關注億速云行業資訊頻道。

向AI問一下細節

免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。

AI

兴安县| 卓资县| 尉氏县| 杭锦旗| 武平县| 揭西县| 赤壁市| 霍山县| 丰县| 塔城市| 那曲县| 贵溪市| 南皮县| 色达县| 乐陵市| 沙湾县| 邹平县| 奈曼旗| 嘉兴市| 临澧县| 呼图壁县| 错那县| 长兴县| 沅江市| 弥勒县| 吉林省| 全南县| 营口市| 依兰县| 承德县| 太原市| 舟曲县| 类乌齐县| 吉首市| 临海市| 衡水市| 织金县| 东乡县| 保亭| 全州县| 平安县|