forked from Nik-Zainal-Group/signature.tools.lib
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsignature.tools.lib.develop.R
80 lines (70 loc) · 3.15 KB
/
signature.tools.lib.develop.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
#file with commands to set up the signature-tools-lib R package
#setwd("~/sandbox/git")
# install.packages("devtools")
# install.packages("roxygen2")
#
# devtools::create("signature.tools.lib")
#setwd("~/sandbox/git/signature.tools.lib")
devtools::use_package("VariantAnnotation")
devtools::use_package("SummarizedExperiment")
devtools::use_package("BiocGenerics")
devtools::use_package("GenomeInfoDb")
devtools::use_package("BSgenome")
devtools::use_package("BSgenome.Hsapiens.UCSC.hg38")
devtools::use_package("BSgenome.Hsapiens.1000genomes.hs37d5")
devtools::use_package("BSgenome.Mmusculus.UCSC.mm10")
devtools::use_package("BSgenome.Cfamiliaris.UCSC.canFam3")
devtools::use_package("NMF")
devtools::use_package("foreach")
devtools::use_package("doParallel")
devtools::use_package("doMC")
devtools::use_package("lpSolve")
devtools::use_package("ggplot2")
devtools::use_package("methods")
devtools::use_package("cluster")
devtools::use_package("stats")
devtools::use_package("NNLM")
devtools::use_package("nnls")
devtools::use_package("GenSA")
devtools::use_package("gmp")
devtools::use_package("plyr")
devtools::use_package("RCircos")
devtools::use_package("scales")
devtools::use_package("GenomicRanges")
devtools::use_package("IRanges")
devtools::load_all()
#add internal data
RS.Breast560 <- read.table("data/Breast560_rearrangement.signatures.txt", sep="\t", header=T, as.is=T, check.names = FALSE)
cosmic30 <- read.table("data/COSMIC30_subs_signatures.txt", sep="\t", header=T, as.is=T, check.names = FALSE)
all_organ_sigs_subs <- read.table("data/2019_01_10_all_PCAWG_sigs_subs.tsv", sep="\t", header=T, as.is=T, check.names = FALSE)
all_organ_sigs_rearr <- read.table("data/2019_01_10_all_PCAWG_sigs_rearr.tsv", sep="\t", header=T, as.is=T, check.names = FALSE)
conversion_matrix_subs <- read.table("data/2019_01_10_ConversionMatrix_subs.tsv", sep="\t", header=T, as.is=T, check.names = FALSE)
conversion_matrix_rearr <- read.table("data/2019_01_10_ConversionMatrix_rearr.tsv", sep="\t", header=T, as.is=T, check.names = FALSE)
load("data/chrominfo.RData")
load("data/chrominfo.snp6.hg19.RData")
usethis::use_data(RS.Breast560,
cosmic30,
chrominfo,
chrominfo.snp6,
all_organ_sigs_subs,
all_organ_sigs_rearr,
conversion_matrix_subs,
conversion_matrix_rearr,
internal = TRUE,overwrite = TRUE)
devtools::document()
devtools::install()
#test all
devtools::test()
#some individual tests
devtools::test(pkg = ".",filter = "ascatToHRDLOH")
devtools::test(pkg = ".",filter = "tabToIndelsClassification")
devtools::test(pkg = ".",filter = "vcfToIndelsClassification")
devtools::test(pkg = ".",filter = "HRDetect")
devtools::test(pkg = ".",filter = "tabToSNVcatalogue")
devtools::test(pkg = ".",filter = "vcfToSNVcatalogue")
devtools::test(pkg = ".",filter = "bedpeToRearrCatalogue")
devtools::test(pkg = ".",filter = "genomePlot")
devtools::test(pkg = ".",filter = "SignatureExtraction")
devtools::test(pkg = ".",filter = "SignaturePartialExtraction")
devtools::test(pkg = ".",filter = "SignatureFit")
devtools::test(pkg = ".",filter = "snvTabToDNVcatalogue")