溫馨提示×

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

密碼登錄×
登錄注冊(cè)×
其他方式登錄
點(diǎn)擊 登錄注冊(cè) 即表示同意《億速云用戶服務(wù)條款》

R語(yǔ)言如何在行上取數(shù)據(jù)子集

發(fā)布時(shí)間:2022-03-17 16:37:47 來(lái)源:億速云 閱讀:263 作者:iii 欄目:開發(fā)技術(shù)

本篇內(nèi)容介紹了“R語(yǔ)言如何在行上取數(shù)據(jù)子集”的有關(guān)知識(shí),在實(shí)際案例的操作過(guò)程中,不少人都會(huì)遇到這樣的困境,接下來(lái)就讓小編帶領(lǐng)大家學(xué)習(xí)一下如何處理這些情況吧!希望大家仔細(xì)閱讀,能夠?qū)W有所成!

subset_row.r  提供列表,在行上取數(shù)據(jù)子集

使用方法:

usage: /share/nas1/huangls/test/TCGA_GEO_immu/scripts/subset_row.r
       [-h] -i input -m idlist -c colname [-o outdir] [-p prefix]
get subset by ID
optional arguments:
  -h, --help            show this help message and exit
  -i input, --input input
                        input data matrix data [required]
  -m idlist, --idlist idlist
                        input idlist [required]
  -c colname, --colname colname
                        match column name in idlist and input [required]
  -o outdir, --outdir outdir
                        output file directory [default cwd]
  -p prefix, --prefix prefix
                        output file name prefix [default demo]

參數(shù)說(shuō)明:

-i 輸入要取子集的數(shù)據(jù),

IDTCGA-D7-A74A-01A-11R-A32D-31TCGA-BR-7704-01A-11R-2055-13TCGA-VQ-A91N-01A-11R-A414-31TCGA-CD-A4MH-01A-11R-A251-31TCGA-BR-7715-11A-01R-2055-13TCGA-EQ-8122-01A-11R-2343-13TCGA-BR-7704-11A-01R-2055-13TCGA-VQ-A91Y-01A-11R-A414-31
TSPAN65.6074344.57416314.3965427.346677.3813427.0949494.6385393.431114
TNMD0.036605000.101440.0696770.0139610.0883280
DPM147.2760269.1515739.0786754.3847634.7598753.7432917.2287645.18002
SCYL32.4808831.1205444.5047552.0202542.4544892.8067471.9564082.751801
C1orf1121.9728380.6577914.0215061.7966650.868562.4898930.5683091.631315
FGR0.6403955.7056550.7682221.8972575.5282512.1286352.7079952.874438
CFH1.0057374.2310423.060052.49926317.276623.24941412.79788.515792
FUCA229.3294935.9760131.4095633.7312921.7049245.0462216.312617.08673
GCLC8.6697327.5759797.3870316.2070494.92036611.92634.3986246.320377
NFYA7.2603468.4350749.56661710.992196.499417.38717.0037426.259676


-m 指定ID列表,-c 指定 ID列表中的列名:

IDlog2FClogCPMPValueFDRregulate
FGR1.151873.5463733.56E-412.75E-39Up
CFH1.1484866.2921132.77E-255.74E-24Up
CD381.6294474.3209281.23E-283.14E-27Up
HSPB63.0133337.4689488.24E-384.57E-36Up
PDK41.058446.5365841.98E-099.88E-09Up
MEOX11.5952762.0779712.50E-361.21E-34Up
COPZ21.0160283.0431141.38E-191.85E-18Up
PRKAR2B1.5292854.0160246.52E-352.73E-33Up
ITGAL1.5956755.1187292.03E-473.33E-45Up

使用舉例:

Rscript $scriptdir/subset_row.r -m ../04.deg/S1_vs_S2.edger.DEG.final.tsv -c ID \
  -i ../02.sample_select/TCGA-STAD_gene_expression_TPM.tsv -p deg_gene_exp_tpm

“R語(yǔ)言如何在行上取數(shù)據(jù)子集”的內(nèi)容就介紹到這里了,感謝大家的閱讀。如果想了解更多行業(yè)相關(guān)的知識(shí)可以關(guān)注億速云網(wǎng)站,小編將為大家輸出更多高質(zhì)量的實(shí)用文章!

向AI問(wèn)一下細(xì)節(jié)

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

AI