[BioC] GenomeGraphs question

Iain Gallagher iaingallagher at btopenworld.com
Wed Feb 8 16:53:36 CET 2012



Dear List

I am plotting some snps using GenomeGraphs to show genomic context. I have a simple question: how do I change the size of the plotting character?

# SCRIPT BEGINS

require(biomaRt) # get and load the library
require(GenomeGraphs) 

# get the snps we need info on
snps <- c('rs136617760', 'rs42494357', 'rs110465273', 'rs42494342',  'rs109809949', 'rs134449285', 'rs137562332', 'rs132841890')

# set up the database we need to query
snpmart <- useMart('snp', dataset = 'btaurus_snp')

# we want to get the snp id, allele, chromosome position and chromosome
# in biomaRt attributes are what you want, filters are what you're using to query the db
# there's load of info available; use listAttributes(snpmart) to see what you can get back
# and listFilters(snpmart) to see what you can put in as a query

snpInfo <- getBM(attributes = c('refsnp_id', 'allele', 'chrom_start', 'chr_name'), filters = 'snp_filter', values = snps, mart = snpmart)

# plot these using GenomeGraphs
# subset out the chromosome 13 snps
thirteenSnps <- subset(snpInfo, snpInfo$chr_name == 13)

# make pretend 'intensity values' for later
ints <- seq(1,6) # could replace with 1/log10p perhaps!
names(ints) <- thirteenSnps$refsnp_id
ints <- as.matrix(ints)

# set up mart
cowmart <- useMart('ensembl', dataset = 'btaurus_gene_ensembl')

# make a plot of the genes on bovine chromosome 13 around the SNPs Mairead found on chromosome 13
ideogram <- makeIdeogram(chromosome = 13)

plusStrand <- makeGeneRegion(chromosome = 13, start = min(thirteenSnps$chrom_start)-5e5, end = min(thirteenSnps$chrom_start)+5e5, strand = "+", biomart = cowmart)

minStrand <- makeGeneRegion( chromosome = 13, start = min(thirteenSnps$chrom_start)-5e5, end = min(thirteenSnps$chrom_start)+5e5, strand = "-", biomart = cowmart)

genomeAxis <- makeGenomeAxis(add53 = TRUE, add35 = TRUE)

snpPlot <- makeGenericArray(intensity = ints, probeStart = thirteenSnps$chrom_start,  dp = DisplayPars(color="darkred", type="point", cex = 2))

gdPlot(list(makeTitle('Chromosome 13'), a = ideogram, b = plusStrand, c = genomeAxis, d = minStrand, e = snpPlot), labelCex = 2, labelRot = 0)

# SCRIPT ENDS

I would like to make the dots representing the snps larger. I tried e.g.

dp = DisplayPars(cex = 2) and dp = DisplayPars(size = 2) in the snpPlot definition without success and I can't work out the answer from the documentation.

Incidentally I'm aware that the above script is perhaps abusing makeGenericArray - it works for me. If anyone has a more 'proper' way of doing this then that would be educational as well.

Best

Iain

> sessionInfo()
R version 2.14.1 (2011-12-22)
Platform: x86_64-pc-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=en_GB.utf8       LC_NUMERIC=C             
 [3] LC_TIME=en_GB.utf8        LC_COLLATE=en_GB.utf8    
 [5] LC_MONETARY=en_GB.utf8    LC_MESSAGES=en_GB.utf8   
 [7] LC_PAPER=C                LC_NAME=C                
 [9] LC_ADDRESS=C              LC_TELEPHONE=C           
[11] LC_MEASUREMENT=en_GB.utf8 LC_IDENTIFICATION=C      

attached base packages:
[1] grid      stats     graphics  grDevices utils     datasets  methods  
[8] base     

other attached packages:
[1] GenomeGraphs_1.14.0 biomaRt_2.10.0      ggplot2_0.9.0      

loaded via a namespace (and not attached):
 [1] colorspace_1.1-1   dichromat_1.2-4    digest_0.5.1       MASS_7.3-16       
 [5] memoise_0.1        munsell_0.3        plyr_1.7.1         proto_0.3-9.2     
 [9] RColorBrewer_1.0-5 RCurl_1.9-5        reshape2_1.2.1     scales_0.1.0      
[13] stringr_0.6        tools_2.14.1       XML_3.9-2         
> 




More information about the Bioconductor mailing list