溫馨提示×

溫馨提示×

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

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

R語言的univariate_cox.r怎么使用

發(fā)布時間:2022-03-21 10:33:50 來源:億速云 閱讀:265 作者:iii 欄目:開發(fā)技術

這篇文章主要介紹了R語言的univariate_cox.r怎么使用的相關知識,內(nèi)容詳細易懂,操作簡單快捷,具有一定借鑒價值,相信大家閱讀完這篇R語言的univariate_cox.r怎么使用文章都會有所收獲,下面我們一起來看看吧。

單因素cox回歸分析

使用方法:

指定單個變量做cox回歸分析,可以是臨床數(shù)據(jù)

$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

參數(shù)說明:

-i 輸入生存數(shù)據(jù)與臨床數(shù)據(jù) 

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怎么使用”這篇文章的內(nèi)容就介紹到這里,感謝各位的閱讀!相信大家對“R語言的univariate_cox.r怎么使用”知識都有一定的了解,大家如果還想學習更多知識,歡迎關注億速云行業(yè)資訊頻道。

向AI問一下細節(jié)

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

AI