GoogleSearch
이 블로그 검색
[BioC] GenomeGraph
라벨:
Informatics
이메일로 전송BlogThis!X에 공유Facebook에서 공유
http://www.alexaplatform.org/examples/GenomeGraphsExample.R
#GenomeGraphs example
#This example was tested with R2.14.1
#Install package if you have not already done so
source("http://bioconductor.org/biocLite.R")
biocLite("GenomeGraphs")
#Load libraries
library(GenomeGraphs)
#Draw the gene and then just one specific transcript of that gene
mart = useMart(biomart = "ensembl", dataset = "hsapiens_gene_ensembl")
gene = makeGene(id = "ENSG00000171988",type = "ensembl_gene_id", biomart = mart)
transcript = makeTranscript(id = "ENST00000399262", type="ensembl_transcript_id", biomart = mart)
gdPlot(list(gene, transcript))
#Draw the gene and then all transcripts of that gene
gene = makeGene(id = "ENSG00000171988",type = "ensembl_gene_id", biomart = mart)
transcript = makeTranscript(id = "ENSG00000171988", type="ensembl_gene_id", biomart = mart)
gdPlot(list(gene, transcript))
#Draw a more annotated figure (using the gene UMPS as an example)
gene = makeGene(id = "ENSG00000114491",type = "ensembl_gene_id", biomart = mart)
transcript = makeTranscript(id = "ENSG00000114491", type="ensembl_gene_id", biomart = mart, dp=DisplayPars(plotId=TRUE, cex=0.5))
#Create an ideogram for the entire chromosome 3
ideog = new("Ideogram", chromosome = "3")
#Create a highlight of the gene position on the ideogram using "absolute coordinates"
ol1 = makeRectangleOverlay(0.60, 0.65, region=c(0.75, 0.82), coords ="absolute", dp = DisplayPars(alpha = .2, fill = "red"))
#Create a highlight of exon 2 of UMPS spanning across the gene and transcript tracks
ol2 = makeRectangleOverlay(124453840, 124454200, region=c(3,4), coords = c("genomic", "absolute"), dp = DisplayPars(alpha = .2, fill = "olivedrab1"))
#Create the plot
gdPlot(list(makeTitle("UMPS"), "Chr3"=ideog, "Gene"=gene, "Transcripts"=transcript), overlays=list(ol1,ol2))
이메일로 전송BlogThis!X에 공유Facebook에서 공유
라벨:
Informatics
Scientist. Husband. Daddy. --- TOLLE. LEGE
외부자료의 인용에 있어 대한민국 저작권법(28조)과 U.S. Copyright Act (17 USC. §107)에 정의된 "저작권물의 공정한 이용원칙 | the U.S. fair use doctrine" 을 따릅니다. 저작권(© 최광민)이 명시된 모든 글과 번역문들에 대해 (1) 복제-배포, (2) 임의수정 및 자의적 본문 발췌, (3) 무단배포를 위한 화면캡처를 금하며, (4) 인용 시 URL 주소 만을 사용할 수 있습니다. [후원 | 운영] [대문으로] [방명록] [옛 방명록] [티스토리 (백업)] [신시내티]
-