golang mysql 简书_非模式物种GO/KEGG富集分析

前言:

微博参与话题 #给你四年时间你也学不会生信#

先前的富集分析教程

本文主要针对非模式物种，但是有参考基因组可用

1. R包安装及database下载

# non-model, but have the genome

> source("https://bioconductor.org/biocLite.R")

> biocLite("AnnotationHub")

> biocLite("biomaRt")

# load package

> library(AnnotationHub)

> library(biomaRt)

# make a orgDb

> hub

这里以桔小实蝇为例

# fruit fly = bactrocera dorsalis

> query(hub, "bactrocera")

搜索后结果如下：

> query(hub, "bactrocera")

AnnotationHub with 9 records

# snapshotDate(): 2018-04-30

# $dataprovider: ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/

# $species: Bactrocera (Bactrocera)_dorsalis, Bactrocera (Bactrocera)_latifrons, Bactrocera (Dacul...

# $rdataclass: OrgDb

# additional mcols(): taxonomyid,