https://github.com/oncogenetics/oncofunco
DNA selection from plates...
Collapse transcripts into gene symbols. Output is a Granges object for defined region in a chromosome.
Required packages:
require(dplyr)
require(GenomicFeatures)
require(TxDb.Hsapiens.UCSC.hg19.knownGene)
require(org.Hs.eg.db)
Get all genes for chr17
allGenesChr17 <- GeneSymbol(chrom = "chr17")
# Define region of interest - ROI
roi <- GRanges(seqnames="chr17",
IRanges(start=41150000,
end=41300000,
names="myRegion"))
# Plot using ggbio - geom_alignment
ggplot(subsetByOverlaps(allGenesChr17,roi)) +
geom_alignment(aes(group=gene_id,fill=strand,col=strand))
Get all genes for chr17:41150000,41300000
subsetGenesChr17 <- GeneSymbol("chr17",41150000,41300000)
ggplot(subsetGenesChr17) +
geom_alignment(aes(group=gene_id,fill=strand,col=strand))
Calculate Odds ratio weighted SNP risk... In progress...