In byandell/CCSanger: Utilities to get and use Sanger data for Collaborative Cross

knitr::opts_chunk$set(echo = TRUE)

v2path <- "~/Documents/Research/attie_alan/DO/AttieDOv2"

query_genes <- 
  qtl2db::create_gene_query_func(
    file.path(v2path, "qtl2db", "mouse_genes.sqlite"))
query_variants <- 
  qtl2db::create_variant_query_func(
    file.path(v2path, "qtl2db", "cc_variants.sqlite"))

chr_id <- "1"
start <- 34
stop <- 35

genes <- CCSanger::get_gene(chr_id, start, stop)

summary(genes)

summary(subset(genes, 34.2, 34.6))

byandell/CCSanger documentation built on May 13, 2019, 9:26 a.m.