-->

GoogleSearch



Scientist. Husband. Daddy. --- TOLLE. LEGE
외부자료의 인용에 있어 대한민국 저작권법(28조)과 U.S. Copyright Act (17 USC. §107)에 정의된 "저작권물의 공정한 이용원칙 | the U.S. fair use doctrine" 을 따릅니다. 저작권(© 최광민)이 명시된 모든 글과 번역문들에 대해 (1) 복제-배포, (2) 임의수정 및 자의적 본문 발췌, (3) 무단배포를 위한 화면캡처를 금하며, (4) 인용 시 URL 주소 만을 사용할 수 있습니다. [후원 | 운영] [대문으로] [방명록] [옛 방명록] [티스토리 (백업)]

이 블로그 검색

[BioC] GenomeGraph

라벨:




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))






라벨:





Scientist. Husband. Daddy. --- TOLLE. LEGE
외부자료의 인용에 있어 대한민국 저작권법(28조)과 U.S. Copyright Act (17 USC. §107)에 정의된 "저작권물의 공정한 이용원칙 | the U.S. fair use doctrine" 을 따릅니다. 저작권(© 최광민)이 명시된 모든 글과 번역문들에 대해 (1) 복제-배포, (2) 임의수정 및 자의적 본문 발췌, (3) 무단배포를 위한 화면캡처를 금하며, (4) 인용 시 URL 주소 만을 사용할 수 있습니다. [후원 | 운영] [대문으로] [방명록] [옛 방명록] [티스토리 (백업)] [신시내티]

-