Skip to content

Instantly share code, notes, and snippets.

View grimbough's full-sized avatar

Mike Smith grimbough

View GitHub Profile
fullXmlQuery <- "<?xml version='1.0' encoding='UTF-8'?> <!DOCTYPE Query>
<Query virtualSchemaName = 'default' uniqueRows = '1' count = '0' datasetConfigVersion = '0.6' header='1' requestid= 'biomaRt'>
<Dataset name = 'hsapiens_gene_ensembl'>
<Attribute name = 'phenotype_description'/>
<Attribute name = 'ensembl_gene_id'/>
<Filter name = 'ensembl_gene_id' value = 'ENSG00000094804'/>
</Dataset>
</Query>"
rcurl_return <- RCurl::postForm(uri = "http://www.ensembl.org:80/biomart/martservice?",
@grimbough
grimbough / BioMart_query.R
Last active October 28, 2017 22:16
Query BioMart with httr and RCurl
fullXmlQuery <- "<?xml version='1.0' encoding='UTF-8'?> <!DOCTYPE Query>
<Query virtualSchemaName = 'default' uniqueRows = '1' count = '0' datasetConfigVersion = '0.6' header='1' requestid= 'biomaRt'>
<Dataset name = 'hsapiens_gene_ensembl'>
<Attribute name = 'ensembl_transcript_id'/>
<Attribute name = 'ensembl_exon_id'/>
<Attribute name = 'rank'/>
<Attribute name = 'genomic_coding_start'/>
<Attribute name = 'cds_start'/>
<Attribute name = '5_utr_start'/>
</Dataset>
@grimbough
grimbough / fastq-memory.R
Last active August 29, 2015 14:18
Examining memory usage in ShortReadQ class
library(ShortRead)
## create a single string with the number of characters defined by 'length'
## and the letter sampled from 'alphabet'
generateString <- function(length = 50, alphabet = c("A","C","G","T") ) {
paste0(sample(alphabet, size = length, replace = TRUE), collapse = "")
}
## create a ShortReadQ object containing 'nreads' reads, each of which has length 'length'
createFastq <- function(nreads, length) {