## ----setup, include = FALSE---------------------------------------------------
library(knitr)
# knitr::opts_knit$set(root.dir = normalizePath("/Users/mbuljan/Documents/Paket_op/DominoEffect_radna/inst/doc"))
#library(devtools)
#knitr::opts_knit$set(root.dir = normalizePath(inst("DominoEffect"))) 

## ----eval=FALSE---------------------------------------------------------------
#  if (!require("BiocManager"))
#      install.packages("BiocManager")
#  BiocManager::install("DominoEffect")

## -----------------------------------------------------------------------------
library(DominoEffect)

## -----------------------------------------------------------------------------
data("TestData", package = "DominoEffect")
data("SnpData", package = "DominoEffect")
data("DominoData", package = "DominoEffect")

## ----eval=FALSE---------------------------------------------------------------
#  DominoEffect(TestData, DominoData, SnpData)

## ----eval=FALSE---------------------------------------------------------------
#  mutation_dataset = read.table ("user_file_with_mutations.txt", header = T)
#  gene_data = read.table ("user_ensembl_gene_list.txt", header = T)
#  snp_data = read.table ("user_population_SNPs_with_frequency.txt", header = T)

## ----message = FALSE, results = "hide"----------------------------------------
hotspot_mutations <- DominoEffect(mutation_dataset = TestData, 
                                  gene_data = DominoData, snp_data = SnpData)

## ----eval= FALSE--------------------------------------------------------------
#  data("TestData", package = "DominoEffect")

## ----echo = FALSE-------------------------------------------------------------
kable(head(TestData), row.names = FALSE)

## ----echo = FALSE-------------------------------------------------------------
kable(head(hotspot_mutations), row.names = FALSE)

## ----eval = FALSE-------------------------------------------------------------
#  hotspot_mutations <- DominoEffect(mutation_dataset, gene_data, snp_data, min_n_muts, MAF_thresh, flanking_region, poisson.thr, percentage.thr, ratio.thr, approach, write_to_file)

## -----------------------------------------------------------------------------
min_n_muts <- 5

## -----------------------------------------------------------------------------
MAF.thr <- 0.01

## -----------------------------------------------------------------------------
flanking_region <- c(200, 300)
flanking_region <- c(300)

## -----------------------------------------------------------------------------
poisson.thr <- 0.01

## -----------------------------------------------------------------------------
percentage.thr <- 0.15

## -----------------------------------------------------------------------------
ratio.thr <- 40

## -----------------------------------------------------------------------------
approach = "percentage"

## ----eval = FALSE-------------------------------------------------------------
#  write_to_file = "YES"

## ----message = FALSE, results = "hide", eval = FALSE--------------------------
#  hotspot_mutations <- DominoEffect(mutation_dataset = TestData,
#                                    gene_data = DominoData, snp_data = SnpData)

## ----message = FALSE----------------------------------------------------------
hotspot_mutations <- identify_hotspots(mutation_dataset = TestData, 
                                       gene_data = DominoData, 
                                       snp_data = SnpData, min_n_muts = 5, 
                                       MAF_thresh = 0.01, 
                                       flanking_region = c(200, 300), 
                                       poisson.thr = 0.01, 
                                       percentage.thr = 0.15, ratio.thr = 45, 
                                       approach = "percentage")

## ----message = FALSE, eval = FALSE--------------------------------------------
#  results_w_annotations <- map_to_func_elem(hotspot_mutations,
#                                            write_to_file = "NO",
#                                            ens_release = "73")

## ----echo = FALSE-------------------------------------------------------------
kable(head(DominoData), row.names = FALSE)

## ----echo = FALSE-------------------------------------------------------------
kable(head(SnpData), row.names = FALSE)

## ----echo = FALSE-------------------------------------------------------------
hotspot_mutations.GPo <- GPo_of_hotspots(hotspot_mutations)
head(hotspot_mutations.GPo)

## -----------------------------------------------------------------------------
sessionInfo()