[AppleScript] 纯文本查看 复制代码
library(org.Hs.eg.db)
allSymbols=unique(mappedRkeys(org.Hs.egSYMBOL))
sirt_genes = grep('SIRT',allSymbols,value = T)
GeneList <- sirt_genes
eg2symbol=toTable(org.Hs.egSYMBOL)
eg2name=toTable(org.Hs.egGENENAME)
eg2map=toTable(org.Hs.egMAP)
eg2start=toTable(org.Hs.egCHRLOC)
eg2end=toTable(org.Hs.egCHRLOCEND)
eg2alias=toTable(org.Hs.egALIAS2EG)
eg2alis_list=lapply(split(eg2alias,eg2alias$gene_id),function(x){paste0(x[,2],collapse = ";")})
#eg2accum=toTable(org.Hs.egACCNUM)
#eg2accum_list=lapply(split(eg2accum,eg2accum$gene_id),function(x){paste0(x[,2],collapse = ";")})
if( GeneList[1] %in% eg2symbol$symbol ){
symbols=GeneList
geneIds=eg2symbol[match(symbols,eg2symbol$symbol),'gene_id']
}else{
geneIds=GeneList
symbols=eg2symbol[match(geneIds,eg2symbol$gene_id),'symbol']
}
geneNames=eg2name[match(geneIds,eg2name$gene_id),'gene_name']
geneMaps=eg2map[match(geneIds,eg2map$gene_id),'cytogenetic_location']
geneStart=eg2start[match(geneIds,eg2name$gene_id),'start_location']
geneEnd=eg2end[match(geneIds,eg2name$gene_id),'end_location']
geneChr=eg2start[match(geneIds,eg2name$gene_id),'Chromosome']
geneAlias=sapply(geneIds,function(x){ifelse(is.null(eg2alis_list[[x]]),"no_alias",eg2alis_list[[x]])})
#geneAssnum=sapply(geneIds,function(x){ifelse(is.null(eg2accum_list[[x]]),"no_accnum",eg2accum_list[[x]])})
createLink <- function(base,val) {
sprintf('<a href="%s" class="btn btn-link" target="_blank" >%s</a>',base,val) ##target="_blank"
}
gene_info=data.frame( symbols=symbols,
geneIds=createLink(paste0("http://www.ncbi.nlm.nih.gov/gene/",geneIds),geneIds),
geneNames=geneNames,
geneAlias=geneAlias,
#accnum=geneAssnum,
map=geneMaps,
chr=geneChr,
start=geneStart,
end=geneEnd,
stringsAsFactors = F
)
#library("xtable")
#print(xtable(gene_info), type="html",include.rownames=F, file='all_gene.anno',sanitize.text.function = force)
file='gene_info_bioconductor.html'
y <- DT::datatable(gene_info,escape = F,rownames=F)
DT::saveWidget(y,file)