# #设置,回声= FALSE,结果=“隐藏 "------------------------------------------ suppressPackageStartupMessages ({suppressMessages({库(BiocStyle)图书馆(HumanTranscriptomeCompendium)图书馆(beeswarm)图书馆(SummarizedExperiment)图书馆(DT ) }) }) ## ---- lklo --------------------------------------------------------------------- 图书馆(HumanTranscriptomeCompendium) genelev = htx_load () genelev如果(rhdf5client: check_hsds())测定(genelev) # #——lksi --------------------------------------------------------------------- 唱= grep(“单一。细胞”,genelev study_title美元,ignore.case = TRUE)(唱)# #——chkspec长度 ------------------------------------------------------------------ sa = genelev study_accession美元(唱)sing2 =唱(——(复制(sa)))长度(sing2)数据表(as.data.frame (colData (genelev [, sing2])),列表选项= (lengthMenu = c(3、5、10、50100 ))) ## ---- lkchar ------------------------------------------------------------------- 大部分= genelev[唱]kpglio = grep(“胶质母细胞瘤”,大部分study_title美元,ignore.case = TRUE) glioGene =散装,kpglio glioGene # #——lkmat -------------------------------------------------------------------- 如果(rhdf5client:: check_hsds ()) {beeswarm(如。matrix(assay(glioGene["ENSG00000138413.13"，1:100]))， pwcol=as.numeric(factor(glioGene$study_title[1:100]))， ylab="IDH1 expression")15000,传说=独特(glioGene study_accession美元[1:10]),坳= 1:2,pch = c (1 1 )) } ## ---- dotx eval = FALSE --------------------------------------------------------- # txlev = htx_load # txlev (genesOnly = FALSE)