Refgenie下载参考基因组
refgenie:参考基因组(阿拉丁)商店 http://refgenie.databio.org/ Here we provide a web interface and a RESTful API to access genome assets for popular reference genome assemblies. 该平台由位于弗吉尼亚大学公共卫生基因组学中心的计算生物学和...
refgenie:参考基因组(阿拉丁)商店 http://refgenie.databio.org/ Here we provide a web interface and a RESTful API to access genome assets for popular reference genome assemblies. 该平台由位于弗吉尼亚大学公共卫生基因组学中心的计算生物学和...
1 2 3 4 5 if (!require("BiocManager", quietly = TRUE)) install.packages("BiocManager") BiocManager::install("STRINGdb") library(STRINGdb) 1、定义要使用的STRING版本、物种,以及PPI阈值分数 1 2 3 4 string_db <- STRINGdb$new(version="11", species=9606, score_threshold=200, input_directory="") 2、示例基因 1 2 3 4 5 6 7 8 9 10 11 12...
1、GO 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66...
一、基因组大小 (1)人类基因组主要由细胞核的23对染色体组成(核基因组),还包括线粒体中的小DNA分子(线粒体基因组)。 (2)单倍体基因组大...
1、不同基因ID转换 1.1 org.Hs.eg.db包 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 library(dplyr) library(org.Hs.eg.db) keytypes(org.Hs.eg.db) # [1] "ACCNUM" "ALIAS" "ENSEMBL" "ENSEMBLPROT" "ENSEMBLTRANS" "ENTREZID" # [7] "ENZYME" "EVIDENCE" "EVIDENCEALL"...
以前通路富集分析直接使用clusterprofiler包,阅读文献发现GSEA分析及可视化较多使用Broad团队研发的工具,现简要学习其(w...
1、背景知识 (1)两种富集分析 基于超几何检验的ORA(over representation analysis)富集分析 ① 假设对转录组分组测序的10000个基因表达数据进...
0、示例数据 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 library(clusterProfiler) library(org.Hs.eg.db) data(geneList, package="DOSE") str(geneList) # Named num [1:12495] 4.57 4.51 4.42 4.14 3.88 ... # - attr(*, "names")= chr [1:12495] "4312" "8318" "10874" "55143" ... gene_ids<-AnnotationDbi::select(org.Hs.eg.db, keys=as.character(names(gene_list)), columns="SYMBOL", #目标格式 keytype="ENTREZID")...
1、关于WGCNA原理 1.1 建立共表达网络 在基因共表达网络中,节点node代表基因,边edge代表两个基因间共表达关系。 若一个基因同时与多个基因...
1、GSVA函数 1 2 3 4 5 6 7 # BiocManager::install("GSVA") library(GSVA) ?gsva gsva(expr = , #metrix格式表达矩阵(行--基因,列--样本) gset.idx.list = , #list格式基因集 method=c("gsva", "ssgsea", "zscore", "plage"), # defaul:gsva kcdf=c("Gaussian", "Poisson", "none"))...