库(BiocStyle)库(HPAanalyze)库(dplyr)
如果您想从HPA系统地下载图像,该怎么做?
CCNB1xml <- hpaXmlGet("ENSG00000134057") CCNB1_ab <- hpaXmlAntibody(CCNB1xml) CCNB1_ab #> # A tibble: 4 x 4 #> id releaseDate releaseVersion RRID #> #> 1 CAB000115 2006-03-13 1.2 #> 2 CAB003804 2006-10-30 2 AB_562272 #> 3 HPA030741 2013-12-05 12 AB_2673586 # b> 4 HPA061448 2016-12-04 16 AB_2684522 CCNB1_expr <- hpaXmlTissueExpr(CCNB1xml) str(CCNB1_expr[[1]]) #>类'tbl_df', 'tbl'和'data.frame': 331 obs。#> $ patientId: chr“1653”“1721”“1725”“598”…#> $年龄:CHR“53”“60”“57”“7”…#> $ sex: chr“男性”“女性”“男性”“男性”…#> $染色:chr NA NA NA…#> $ intensity: chr NA NA NA NA…#> $ quantity: chr NA NA NA NA…#> $ location: chr NA NA NA NA…#> $ imageUrl: chr“http://v18.proteinatlas.org/images/115/2043_B_4_5.jpg”“http://v18.proteinatlas.org/images/115/2043_B_6_5.jpg”“http://v18.proteinatlas.org/images/115/2043_B_5_5.jpg”“http://v18.proteinatlas.org/images/115/2043_A_2_2.jpg”…#> $ snomedCode1: chr“M-00100”“M-00100”“M-00100”“M-00100”… #> $ snomedCode2 : chr "T-93000" "T-93000" "T-93000" "T-66000" ... #> $ snomedCode3 : chr NA NA NA NA ... #> $ snomedCode4 : chr NA NA NA NA ... #> $ snomedCode5 : chr NA NA NA NA ... #> $ tissueDescription1: chr "Normal tissue, NOS" "Normal tissue, NOS" "Normal tissue, NOS" "Normal tissue, NOS" ... #> $ tissueDescription2: chr "Adrenal gland" "Adrenal gland" "Adrenal gland" "Appendix" ... #> $ tissueDescription3: chr NA NA NA NA ... #> $ tissueDescription4: chr NA NA NA NA ... #> $ tissueDescription5: chr NA NA NA NA ...
dir.create("img") for (i in 1:nrow(CCNB1_expr[[1]])){下载。文件(CCNB1_expr[[1]]美元imageUrl[我],destfile = paste0(“img / CCNB1_ab $ id[1],“_”,CCNB1_expr[[1]]美元patientId[我],“_”,CCNB1_expr[[1]]美元tissueDescription2[我],“jpg”),模式=“白平衡”)}
陈安,2018-2022
请列举:Tran, a.n., Dussaq, a.m., Kennell, T.等人。HPAanalyze:一个R包,便于检索和分析人类蛋白质图谱数据。BMC生物信息学20,463 (2019)https://doi.org/10.1186/s12859-019-3059-z