使用clusterProfiler给非模式物种做富集分析的一些参
2021-02-12 本文已影响0人
小明的数据分析笔记本
玉米虽然是非模式物种，但是研究的人也挺多的，注释信息也非常完整了
- 2 https://bioconductor.statistik.tu-dortmund.de/packages/2.12/bioc/vignettes/clusterProfiler/inst/doc/clusterProfiler_for_unsupported_organisms.pdf
这个是一个细菌的例子 - 3-1 构建自己物种的OrgDB
https://rdrr.io/bioc/AnnotationForge/man/makeOrgPackage.html
代码复制过来
if(interactive()){
## Makes an organism package for Zebra Finch data.frames:
finchFile <- system.file("extdata","finch_info.txt",package="AnnotationForge")
finch <- read.table(finchFile,sep="\t")
## not that this is how it should always be, but that it *could* be this way.
fSym <- finch[,c(2,3,9)]
fSym <- fSym[fSym[,2]!="-",]
fSym <- fSym[fSym[,3]!="-",]
colnames(fSym) <- c("GID","SYMBOL","GENENAME")
fChr <- finch[,c(2,7)]
fChr <- fChr[fChr[,2]!="-",]
colnames(fChr) <- c("GID","CHROMOSOME")
finchGOFile <- system.file("extdata","GO_finch.txt",package="AnnotationForge")
fGO <- read.table(finchGOFile,sep="\t")
fGO <- fGO[fGO[,2]!="",]
fGO <- fGO[fGO[,3]!="",]
colnames(fGO) <- c("GID","GO","EVIDENCE")
makeOrgPackage(gene_info=fSym, chromosome=fChr, go=fGO,
version="0.1",
maintainer="Some One <so@someplace.org>",
author="Some One <so@someplace.org>",
outputDir = ".",
tax_id="59729",
genus="Taeniopygia",
species="guttata",
goTable="go")
## then you can call install.packages based on the return value
install.packages("./org.Tguttata.eg.db", repos=NULL)
}
- 3-2 构建自己物种的OrgDB
https://rdrr.io/bioc/AnnotationForge/src/inst/script/makeNoSchemaDb.R
代码复制过来
## test data lets load up a couple tab files from /extdata...
## file finch_info.txt now has data I can use (for examples)
library(AnnotationForge)
finchFile <- system.file("extdata","finch_info.txt",package="AnnotationForge")
finch <- read.table(finchFile,sep="\t")
## not that this is how it should always be, but that it *could* be this way.
fSym <- finch[,c(2,3,9)]
fSym <- fSym[fSym[,2]!="-",]
fSym <- fSym[fSym[,3]!="-",]
colnames(fSym) <- c("GID","SYMBOL","GENENAME")
fChr <- finch[,c(2,7)]
fChr <- fChr[fChr[,2]!="-",]
colnames(fChr) <- c("GID","CHROMOSOME")
finchGOFile <- system.file("extdata","GO_finch.txt",package="AnnotationForge")
fGO <- read.table(finchGOFile,sep="\t")
fGO <- fGO[fGO[,2]!="",]
fGO <- fGO[fGO[,3]!="",]
colnames(fGO) <- c("GID","GO","EVIDENCE")
## Now make a list
## data <- list(gene_info=fSym, chromosome=fChr, go=fGO)
genus <- "Taeniopygia"
species <- "guttataTestingSubset"
dbName <- AnnotationForge:::.generateOrgDbName(genus,species)
## this becomes the file name for the DB
dbfile <- paste(dbName, ".sqlite", sep="")
## Test DB building:
## AnnotationForge:::makeOrgDbFromDataFrames(data, genus, species, dbfile)
# debug(AnnotationForge:::.makeAnnDbPkg)
# debug(AnnotationForge:::.createAnnotPackage)
## or test pkg building
AnnotationForge:::makeOrgPackage(gene_info=fSym, chromosome=fChr, go=fGO,
version="0.1",
maintainer="Some One <so@someplace.org>",
author="Some One <so@someplace.org>",
outputDir = ".",
tax_id="59729",
genus=genus,
species=species)
AnnotationForge:::makeOrgPackage(gene_info=fSym, chromosome=fChr, go=fGO,
version="0.1",
maintainer="Some One <so@someplace.org>",
author="Some One <so@someplace.org>",
outputDir = ".",
tax_id="59729",
genus=genus,
species=species,
goTable="go")
## then you can install on the return value
install.packages("./org.Tguttata.eg.db", repos=NULL)
## NEXT UP: lets make an actual template in AnnotationDbi so that I
## can start making these things as packages.
library(org.Tguttata.eg.db)
#debug(AnnotationDbi:::.noSchemaCols)
columns(org.Tguttata.eg.db)
## so that change will also work for keytypes
keytypes(org.Tguttata.eg.db)
## Now for keys I need to read carefully.
## I want to do the same thing for .keys that I did for
library(org.Tguttata.eg.db)
# debug(AnnotationDbi:::.noSchemaKeys)
keys(org.Tguttata.eg.db, "CHROMOSOME")
head(keys(org.Tguttata.eg.db, "GID"))
head(keys(org.Tguttata.eg.db, "SYMBOL", pattern="BDNF"))
head(keys(org.Tguttata.eg.db, "GID", pattern="BDNF", column="SYMBOL"))
head(keys(org.Tguttata.eg.db, "SYMBOL", column="GID"))
## Now I just need select() to work...
library(org.Tguttata.eg.db)
## debug(AnnotationDbi:::.noSchemaSelect)
## TODO: add check to the keys argument to make sure its a character()
select(org.Tguttata.eg.db, keys="100008579", columns="SYMBOL", keytype="GID")
##
select(org.Tguttata.eg.db, keys="100008579", columns=c("SYMBOL","CHROMOSOME"), keytype="GID")
## This one should warn about the number or rows - it does
select(org.Tguttata.eg.db, keys="100008579", columns="GO", keytype="GID")
## now fixed
select(org.Tguttata.eg.db, keys="100008579", columns=c("GO","EVIDENCE"), keytype="GID")
## What if there is only one table to visit?
select(org.Tguttata.eg.db, keys="BDNF", columns="GENENAME", keytype="SYMBOL")
## TODO:
## 1) Add checks for all arguments (especially data) - DONE
## 2) Add code to look for and extract ALL extra GIDs - DONE
## 3) Add code to filter GO IDs :
## A - add arg to function to say which table should get GO treatment. - DONE
## B - add special checks for GO frame (special checker)- DONE
## - only three cols required (GID, GO, EVIDENCE)- DONE
## - those cols must have names: GID, GO, EVIDENCE- DONE
## C - add function to drop old go table, and make 2 new ones: -DONE
## D - use .filterGOFrame to remove stuff that's too new from go->populate - DONE
## E - make newer .expandGOFRAME to make a go_all table -DONE
## - (GID, GOALL, EVIDENCEALL, ONTOLOGYALL) - yep
## - require(GO.db); select(GO.db, go$GO, "ONTOLOGY") -yep
## 4) Add indexing to relevant fields (_id, go_id and all others) - DONE
## 5.5) change the argument on the external function from a list to a ... Then internally convert that to a list use: extraArgs <- list(...) - DONE
## 6) Add manual page - DONE
## 5) Add unit tests for new helpers - Need some help
## 6.5) Add more checks:
##1) GO terms should look like: "GO:DDDDDDD";
## 2) no row names on the data.frames. - DONE
## 7) Add vignette for makeOrgPackage() and for makeOrgPackageFromNCBI()
欢迎大家关注我的公众号
小明的数据分析笔记本
公众号 主要分享：1、R语言和python做数据分析和数据可视化的简单小例子；2、园艺植物相关转录组学、基因组学、群体遗传学文献阅读笔记！
