## Not run:
# library(TxDb.Hsapiens.UCSC.hg19.knownGene)
#
# # Create a missense mutations dataset in TP53
# set.seed(1)
# locations = sample(7575000:7580000, size = 25)
# mut.df <- data.frame(chrom = "chr17",
# start = locations,
# end = locations,
# strand = "-",
# fill = letters[1:25])
#
# # Generating sample data using TP53
# tx.db <- TxDb.Hsapiens.UCSC.hg19.knownGene
# seqs.hg19 <- seqlengths(tx.db)
# data(genesymbol, package = "biovizBase")
#
# # Sample GRanges data with domains for P53
# sample <- GRanges(seqnames = "chr17",
# IRanges(start = c(7576884, 7577499, 7579707),
# end = c(7577102, 7579403, 7579899)),
# strand = "-",
# domain = c("P53_tetramer", "P53", "P53_TAD"))
#
# # Create PFAM objects
# objs <- makePFAMObjs()
# desc <- objs$desc
# ids <- objs$ids
#
# int.opt <- c("Whole", "Custom")
# disp.opt <- c("All", "Reduce")
# dom.opt <- c("All", "Longest")
#
# # GUI
# mepGUI(m.data = mut.df, tx.db = tx.db, gene.loc = genesymbol,
# seq.lens = seqs.hg19, pfam.desc = desc, pfam.ids = ids,
# pfam.gr = sample, int.opt = int.opt, disp.opt = disp.opt,
# dom.opt = dom.opt)
#
# # Text Interface
# mepTxtInt(m.data = mut.df, tx.db = tx.db, gene.loc = genesymbol,
# seq.lens = seqs.hg19, pfam.desc = desc, pfam.ids = ids,
# pfam.gr = sample, int.opt = int.opt, disp.opt = disp.opt,
# dom.opt = dom.opt)
# ## End(Not run)
Run the code above in your browser using DataLab