This page was generated on 2018-04-12 13:42:52 -0400 (Thu, 12 Apr 2018).
R version 3.4.4 (2018-03-15) -- "Someone to Lean On"
Copyright (C) 2018 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin15.6.0 (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> BiocGenerics:::testPackage("PureCN")
Attaching package: 'BiocGenerics'
The following objects are masked from 'package:parallel':
clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
clusterExport, clusterMap, parApply, parCapply, parLapply,
parLapplyLB, parRapply, parSapply, parSapplyLB
The following objects are masked from 'package:stats':
IQR, mad, sd, var, xtabs
The following objects are masked from 'package:base':
Filter, Find, Map, Position, Reduce, anyDuplicated, append,
as.data.frame, cbind, colMeans, colSums, colnames, do.call,
duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted,
lapply, lengths, mapply, match, mget, order, paste, pmax, pmax.int,
pmin, pmin.int, rank, rbind, rowMeans, rowSums, rownames, sapply,
setdiff, sort, table, tapply, union, unique, unsplit, which,
which.max, which.min
Attaching package: 'S4Vectors'
The following object is masked from 'package:base':
expand.grid
Welcome to Bioconductor
Vignettes contain introductory material; view with
'browseVignettes()'. To cite Bioconductor, see
'citation("Biobase")', and for packages 'citation("pkgname")'.
Attaching package: 'matrixStats'
The following objects are masked from 'package:Biobase':
anyMissing, rowMedians
Attaching package: 'DelayedArray'
The following objects are masked from 'package:matrixStats':
colMaxs, colMins, colRanges, rowMaxs, rowMins, rowRanges
The following object is masked from 'package:base':
apply
Attaching package: 'Biostrings'
The following object is masked from 'package:DelayedArray':
type
The following object is masked from 'package:base':
strsplit
Attaching package: 'VariantAnnotation'
The following object is masked from 'package:base':
tabulate
Loading required package: GenomicFeatures
Loading required package: AnnotationDbi
'select()' returned many:1 mapping between keys and columns
'select()' returned many:1 mapping between keys and columns
WARN [2018-04-12 08:11:27] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
WARN [2018-04-12 08:11:29] Chromosome naming style of txdb file (UCSC) was different from interval file (NCBI).
'select()' returned many:1 mapping between keys and columns
'select()' returned many:1 mapping between keys and columns
WARN [2018-04-12 08:11:33] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
FATAL [2018-04-12 08:11:36] tumor.coverage.file and gc.gene.file do not align.
FATAL [2018-04-12 08:11:36]
FATAL [2018-04-12 08:11:36] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:11:36] parameters (PureCN 1.8.1).
WARN [2018-04-12 08:11:37] No mappability scores provided.
INFO [2018-04-12 08:11:37] Calculating GC-content...
WARN [2018-04-12 08:11:37] No mappability scores provided.
INFO [2018-04-12 08:11:37] Calculating GC-content...
FATAL [2018-04-12 08:11:37] Interval coordinates should start at 1, not at 0
FATAL [2018-04-12 08:11:37]
FATAL [2018-04-12 08:11:37] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:11:37] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:11:38] Tiling off-target regions to an average width of 200000.
WARN [2018-04-12 08:11:38] No mappability scores provided.
INFO [2018-04-12 08:11:38] Calculating GC-content...
WARN [2018-04-12 08:11:38] Intervals contain off-target regions. Will not change intervals.
WARN [2018-04-12 08:11:38] No mappability scores provided.
INFO [2018-04-12 08:11:38] Calculating GC-content...
INFO [2018-04-12 08:11:39] Calculating GC-content...
INFO [2018-04-12 08:11:39] Calculating GC-content...
WARN [2018-04-12 08:11:39] No mappability scores provided.
INFO [2018-04-12 08:11:39] Calculating GC-content...
WARN [2018-04-12 08:11:39] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2018-04-12 08:11:39] No mappability scores provided.
INFO [2018-04-12 08:11:39] Calculating GC-content...
FATAL [2018-04-12 08:11:39] Chromosome naming style of interval file unknown, should be UCSC.
FATAL [2018-04-12 08:11:39]
FATAL [2018-04-12 08:11:39] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:11:39] parameters (PureCN 1.8.1).
WARN [2018-04-12 08:11:39] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2018-04-12 08:11:39] Chromosome naming style of mappability file (NCBI) was different from reference (UCSC).
INFO [2018-04-12 08:11:40] Calculating GC-content...
Minimum 3 supporting reads.
Expected allelic fraction 0.5.
Minimum 3 supporting reads.
Expected allelic fraction 0.5.
Minimum 3 supporting reads.
Expected allelic fraction 0.125.
Minimum 5 supporting reads.
Expected allelic fraction 0.02.
Minimum 3 supporting reads.
FATAL [2018-04-12 08:11:40] Need either f or purity and ploidy.
FATAL [2018-04-12 08:11:40]
FATAL [2018-04-12 08:11:40] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:11:40] parameters (PureCN 1.8.1).
Minimum 3 supporting reads.
FATAL [2018-04-12 08:11:40] f not in expected range.
FATAL [2018-04-12 08:11:40]
FATAL [2018-04-12 08:11:40] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:11:40] parameters (PureCN 1.8.1).
FATAL [2018-04-12 08:11:40] coverage not in expected range (>=2)
FATAL [2018-04-12 08:11:40]
FATAL [2018-04-12 08:11:40] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:11:40] parameters (PureCN 1.8.1).
Minimum 2 supporting reads.
FATAL [2018-04-12 08:11:40] purity not in expected range.
FATAL [2018-04-12 08:11:40]
FATAL [2018-04-12 08:11:40] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:11:40] parameters (PureCN 1.8.1).
Minimum 2 supporting reads.
FATAL [2018-04-12 08:11:40] ploidy not in expected range.
FATAL [2018-04-12 08:11:40]
FATAL [2018-04-12 08:11:40] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:11:40] parameters (PureCN 1.8.1).
FATAL [2018-04-12 08:11:40] cell.fraction not in expected range.
FATAL [2018-04-12 08:11:40]
FATAL [2018-04-12 08:11:40] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:11:40] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:11:42] No on_target column in gc.gene.file. Recreate this file with IntervalFile.R.
INFO [2018-04-12 08:11:42] No mappability column in gc.gene.file.
INFO [2018-04-12 08:11:45] ------------------------------------------------------------
INFO [2018-04-12 08:11:45] PureCN 1.8.1
INFO [2018-04-12 08:11:45] ------------------------------------------------------------
INFO [2018-04-12 08:11:45] Loading coverage files...
INFO [2018-04-12 08:11:46] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-04-12 08:11:46] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:11:46] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:11:46] Removing 22 small (< 5bp) targets.
INFO [2018-04-12 08:11:46] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-04-12 08:11:46] No normalDB provided. Provide one for better results.
INFO [2018-04-12 08:11:46] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-04-12 08:11:46] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-04-12 08:11:46] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-04-12 08:11:46] No gc.gene.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-04-12 08:11:46] Loading VCF...
INFO [2018-04-12 08:11:46] LIB-02240e4 is tumor in VCF file.
INFO [2018-04-12 08:11:46] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-04-12 08:11:46] Found 2331 variants in VCF file.
INFO [2018-04-12 08:11:46] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:11:46] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:11:46] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:11:46] Removing 16 low quality variants with BQ < 25.
INFO [2018-04-12 08:11:46] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-04-12 08:11:46] 1.0% of targets contain variants.
INFO [2018-04-12 08:11:46] Removing 2122 variants outside intervals.
INFO [2018-04-12 08:11:47] Found SOMATIC annotation in VCF.
INFO [2018-04-12 08:11:47] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-04-12 08:11:47] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-04-12 08:11:47] Sample sex: ?
INFO [2018-04-12 08:11:47] Segmenting data...
INFO [2018-04-12 08:11:47] Setting undo.SD parameter to 1.000000.
INFO [2018-04-12 08:11:50] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-04-12 08:11:51] Using 123 variants.
INFO [2018-04-12 08:11:51] Mean standard deviation of log-ratios: 0.37
INFO [2018-04-12 08:11:51] 2D-grid search of purity and ploidy...
INFO [2018-04-12 08:11:53] Local optima: 0.63/1.9, 0.4/2.2, 0.5/2
INFO [2018-04-12 08:11:53] Testing local optimum 1/3 at purity 0.63 and total ploidy 1.90...
INFO [2018-04-12 08:11:54] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:11:55] Optimized purity: 0.65
INFO [2018-04-12 08:11:55] Testing local optimum 2/3 at purity 0.40 and total ploidy 2.20...
INFO [2018-04-12 08:11:56] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2018-04-12 08:11:57] Optimized purity: 0.40
INFO [2018-04-12 08:11:57] Testing local optimum 3/3 at purity 0.50 and total ploidy 2.00...
INFO [2018-04-12 08:11:59] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:12:00] Optimized purity: 0.65
INFO [2018-04-12 08:12:00] Done.
INFO [2018-04-12 08:12:00] ------------------------------------------------------------
FATAL [2018-04-12 08:12:02] fun.countMutation not a function.
FATAL [2018-04-12 08:12:02]
FATAL [2018-04-12 08:12:02] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:12:02] parameters (PureCN 1.8.1).
FATAL [2018-04-12 08:12:02] exclude not a GRanges object.
FATAL [2018-04-12 08:12:02]
FATAL [2018-04-12 08:12:02] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:12:02] parameters (PureCN 1.8.1).
FATAL [2018-04-12 08:12:02] callable not a GRanges object.
FATAL [2018-04-12 08:12:02]
FATAL [2018-04-12 08:12:02] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:12:02] parameters (PureCN 1.8.1).
FATAL [2018-04-12 08:12:11] No gc_bias column in gc.gene.file.
FATAL [2018-04-12 08:12:11]
FATAL [2018-04-12 08:12:11] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:12:11] parameters (PureCN 1.8.1).
WARN [2018-04-12 08:12:11] tumor.coverage.file and gc.gene.file do not align.
INFO [2018-04-12 08:12:35] Reading Sample1_PureCN.rds...
INFO [2018-04-12 08:12:35] Reading Sample1_PureCN.rds...
INFO [2018-04-12 08:12:35] Reading Sample1_PureCN.rds...
INFO [2018-04-12 08:12:36] Reading Sample1_PureCN.rds...
INFO [2018-04-12 08:12:36] Reading Sample1_PureCN.rds...
INFO [2018-04-12 08:12:36] Reading Sample1_PureCN.rds...
FATAL [2018-04-12 08:12:36] Purity or Ploidy not numeric or in expected range.
FATAL [2018-04-12 08:12:36]
FATAL [2018-04-12 08:12:36] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:12:36] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:12:36] Reading Sample1_PureCN.rds...
INFO [2018-04-12 08:12:36] Reading Sample1_PureCN.rds...
FATAL [2018-04-12 08:12:37] 'Failed' column in Sample1_PureCN.csv not logical(1).
FATAL [2018-04-12 08:12:37]
FATAL [2018-04-12 08:12:37] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:12:37] parameters (PureCN 1.8.1).
WARN [2018-04-12 08:12:38] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:12:38] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:12:41] Pooling example_normal.txt, example_normal2.txt.
INFO [2018-04-12 08:12:42] Coverage file does not contain read count information, using total coverage for calculating log-ratios.
WARN [2018-04-12 08:12:53] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:12:53] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:12:56] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:12:56] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:12:58] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:12:58] Allosome coverage missing, cannot determine sex.
FATAL [2018-04-12 08:12:58] Length of normal.coverage.files and sex different
FATAL [2018-04-12 08:12:58]
FATAL [2018-04-12 08:12:58] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:12:58] parameters (PureCN 1.8.1).
WARN [2018-04-12 08:13:00] tumor.coverage.file and gc.gene.file do not align.
WARN [2018-04-12 08:13:02] Target intervals were not sorted.
WARN [2018-04-12 08:13:03] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:13:03] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:13:03] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:13:05] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:13:05] Sample sex: NA
FATAL [2018-04-12 08:13:06] tumor.coverage.file and normalDB do not align.
FATAL [2018-04-12 08:13:06]
FATAL [2018-04-12 08:13:06] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:06] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:08] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:13:08] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:13:08] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:13:08] Removing 16 low quality variants with BQ < 25.
INFO [2018-04-12 08:13:08] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:13:08] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:13:08] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:13:08] Removing 6 blacklisted variants.
INFO [2018-04-12 08:13:08] Removing 16 low quality variants with BQ < 25.
WARN [2018-04-12 08:13:08] MuTect stats file lacks contig and position columns.
INFO [2018-04-12 08:13:08] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:13:09] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:13:09] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:13:09] Removing 16 low quality variants with BQ < 25.
WARN [2018-04-12 08:13:09] MuTect stats file and VCF file do not align perfectly. Will remove 2325 unmatched variants.
WARN [2018-04-12 08:13:09] MuTect stats file lacks failure_reasons column. Keeping all variants listed in stats file.
INFO [2018-04-12 08:13:09] Removing 0 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:13:09] Removing 0 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:13:09] Removing 0 low quality variants with BQ < 25.
WARN [2018-04-12 08:13:11] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:13:11] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:13:11] Mean coverages: chrX: 108.58, chrY: 108.58, chr1-22: 104.87.
INFO [2018-04-12 08:13:11] Mean coverages: chrX: 108.58, chrY: 2.17, chr1-22: 104.87.
INFO [2018-04-12 08:13:11] Mean coverages: chrX: 108.58, chrY: 5.17, chr1-22: 104.87.
INFO [2018-04-12 08:13:12] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-04-12 08:13:12] No germline variants in VCF.
WARN [2018-04-12 08:13:15] Found 4 overlapping intervals, starting at line 2.
INFO [2018-04-12 08:13:16] Loading coverage data...
INFO [2018-04-12 08:13:17] Mean target coverages: 112X (tumor) 99X (normal).
INFO [2018-04-12 08:13:17] Mean target coverages: 112X (tumor) 43X (normal).
INFO [2018-04-12 08:13:18] ------------------------------------------------------------
INFO [2018-04-12 08:13:18] PureCN 1.8.1
INFO [2018-04-12 08:13:18] ------------------------------------------------------------
INFO [2018-04-12 08:13:18] Loading coverage files...
INFO [2018-04-12 08:13:19] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-04-12 08:13:19] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:13:19] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:13:19] Removing 22 small (< 5bp) targets.
INFO [2018-04-12 08:13:19] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-04-12 08:13:19] No normalDB provided. Provide one for better results.
INFO [2018-04-12 08:13:19] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-04-12 08:13:19] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-04-12 08:13:19] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-04-12 08:13:19] No gc.gene.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-04-12 08:13:19] Sample sex: ?
INFO [2018-04-12 08:13:19] Segmenting data...
INFO [2018-04-12 08:13:19] Target weights found, will use weighted CBS.
INFO [2018-04-12 08:13:19] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-04-12 08:13:22] Mean standard deviation of log-ratios: 0.37
INFO [2018-04-12 08:13:22] 2D-grid search of purity and ploidy...
INFO [2018-04-12 08:13:22] Local optima: 0.65/1.9, 0.33/1.9, 0.38/2.2, 0.27/2.1, 0.29/2, 0.78/2
INFO [2018-04-12 08:13:22] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2018-04-12 08:13:24] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-04-12 08:13:24] Recalibrating log-ratios...
INFO [2018-04-12 08:13:24] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-04-12 08:13:24] Recalibrating log-ratios...
INFO [2018-04-12 08:13:24] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-04-12 08:13:25] Recalibrating log-ratios...
INFO [2018-04-12 08:13:25] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-04-12 08:13:25] Testing local optimum 3/6 at purity 0.38 and total ploidy 2.20...
INFO [2018-04-12 08:13:27] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-04-12 08:13:27] Recalibrating log-ratios...
INFO [2018-04-12 08:13:27] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-04-12 08:13:27] Recalibrating log-ratios...
INFO [2018-04-12 08:13:27] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-04-12 08:13:28] Recalibrating log-ratios...
INFO [2018-04-12 08:13:28] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-04-12 08:13:28] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-04-12 08:13:29] Recalibrating log-ratios...
INFO [2018-04-12 08:13:29] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-04-12 08:13:29] Recalibrating log-ratios...
INFO [2018-04-12 08:13:29] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-04-12 08:13:30] Recalibrating log-ratios...
INFO [2018-04-12 08:13:30] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-04-12 08:13:30] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2018-04-12 08:13:31] Done.
INFO [2018-04-12 08:13:31] ------------------------------------------------------------
INFO [2018-04-12 08:13:31] Reading /tmp/RtmpmsF4Cj/file4ecb1627c83d.rds...
FATAL [2018-04-12 08:13:31] runAbsoluteCN was run without a VCF file.
FATAL [2018-04-12 08:13:31]
FATAL [2018-04-12 08:13:31] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:31] parameters (PureCN 1.8.1).
FATAL [2018-04-12 08:13:31] runAbsoluteCN was run without a VCF file.
FATAL [2018-04-12 08:13:31]
FATAL [2018-04-12 08:13:31] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:31] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:31] ------------------------------------------------------------
INFO [2018-04-12 08:13:31] PureCN 1.8.1
INFO [2018-04-12 08:13:31] ------------------------------------------------------------
INFO [2018-04-12 08:13:31] Loading coverage files...
FATAL [2018-04-12 08:13:32] Need a normal coverage file if log.ratio and seg.file are not provided.
FATAL [2018-04-12 08:13:32]
FATAL [2018-04-12 08:13:32] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:32] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
INFO [2018-04-12 08:13:32] PureCN 1.8.1
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
FATAL [2018-04-12 08:13:32] min.ploidy or max.ploidy not within expected range.
FATAL [2018-04-12 08:13:32]
FATAL [2018-04-12 08:13:32] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:32] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
INFO [2018-04-12 08:13:32] PureCN 1.8.1
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
FATAL [2018-04-12 08:13:32] min.ploidy or max.ploidy not within expected range.
FATAL [2018-04-12 08:13:32]
FATAL [2018-04-12 08:13:32] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:32] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
INFO [2018-04-12 08:13:32] PureCN 1.8.1
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
INFO [2018-04-12 08:13:32] PureCN 1.8.1
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
INFO [2018-04-12 08:13:32] PureCN 1.8.1
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
FATAL [2018-04-12 08:13:32] test.purity not within expected range.
FATAL [2018-04-12 08:13:32]
FATAL [2018-04-12 08:13:32] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:32] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
INFO [2018-04-12 08:13:32] PureCN 1.8.1
INFO [2018-04-12 08:13:32] ------------------------------------------------------------
INFO [2018-04-12 08:13:32] Loading coverage files...
FATAL [2018-04-12 08:13:33] Tumor and normal are identical. This won't give any meaningful results
FATAL [2018-04-12 08:13:33] and I'm stopping here.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] Loading coverage files...
FATAL [2018-04-12 08:13:33] Length of log.ratio different from tumor coverage.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
FATAL [2018-04-12 08:13:33] prior.purity must have the same length as test.purity.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
FATAL [2018-04-12 08:13:33] min.gof not within expected range or format.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
FATAL [2018-04-12 08:13:33] prior.purity not within expected range or format.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
FATAL [2018-04-12 08:13:33] prior.purity must add to 1. Sum is 1.5
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
FATAL [2018-04-12 08:13:33] max.homozygous.loss not within expected range or format.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
FATAL [2018-04-12 08:13:33] prior.K not within expected range or format.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
FATAL [2018-04-12 08:13:33] prior.contamination not within expected range or format.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
FATAL [2018-04-12 08:13:33] Iterations not in the expected range from 10 to 250.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
FATAL [2018-04-12 08:13:33] Iterations not in the expected range from 10 to 250.
FATAL [2018-04-12 08:13:33]
FATAL [2018-04-12 08:13:33] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:33] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:33] PureCN 1.8.1
INFO [2018-04-12 08:13:33] ------------------------------------------------------------
INFO [2018-04-12 08:13:34] ------------------------------------------------------------
INFO [2018-04-12 08:13:34] PureCN 1.8.1
INFO [2018-04-12 08:13:34] ------------------------------------------------------------
INFO [2018-04-12 08:13:34] Loading coverage files...
FATAL [2018-04-12 08:13:34] Interval files in normal and tumor different.
FATAL [2018-04-12 08:13:34]
FATAL [2018-04-12 08:13:34] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:13:34] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:13:34] ------------------------------------------------------------
INFO [2018-04-12 08:13:34] PureCN 1.8.1
INFO [2018-04-12 08:13:34] ------------------------------------------------------------
INFO [2018-04-12 08:13:34] Loading coverage files...
INFO [2018-04-12 08:13:35] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-04-12 08:13:35] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:13:35] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:13:35] Removing 22 small (< 5bp) targets.
INFO [2018-04-12 08:13:35] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-04-12 08:13:35] No normalDB provided. Provide one for better results.
INFO [2018-04-12 08:13:35] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-04-12 08:13:35] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-04-12 08:13:35] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-04-12 08:13:35] No gc.gene.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-04-12 08:13:35] Loading VCF...
INFO [2018-04-12 08:13:35] LIB-02240e4 is tumor in VCF file.
INFO [2018-04-12 08:13:35] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-04-12 08:13:35] Found 2331 variants in VCF file.
INFO [2018-04-12 08:13:35] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:13:36] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:13:36] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:13:36] Removing 16 low quality variants with BQ < 25.
INFO [2018-04-12 08:13:36] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-04-12 08:13:36] 1.0% of targets contain variants.
INFO [2018-04-12 08:13:36] Removing 2122 variants outside intervals.
INFO [2018-04-12 08:13:36] Found SOMATIC annotation in VCF.
INFO [2018-04-12 08:13:36] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-04-12 08:13:36] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-04-12 08:13:36] Sample sex: ?
INFO [2018-04-12 08:13:36] Segmenting data...
INFO [2018-04-12 08:13:36] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-04-12 08:13:40] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-04-12 08:13:40] Using 123 variants.
INFO [2018-04-12 08:13:40] Mean standard deviation of log-ratios: 0.37
INFO [2018-04-12 08:13:40] 2D-grid search of purity and ploidy...
INFO [2018-04-12 08:13:49] Local optima: 0.62/1.9, 0.38/2.2, 0.3/1.9, 0.48/2
INFO [2018-04-12 08:13:49] Testing local optimum 1/4 at purity 0.62 and total ploidy 1.90...
INFO [2018-04-12 08:13:50] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:13:51] Optimized purity: 0.65
INFO [2018-04-12 08:13:51] Testing local optimum 2/4 at purity 0.38 and total ploidy 2.20...
INFO [2018-04-12 08:13:53] Fitting variants for purity 0.40, tumor ploidy 2.43 and contamination 0.01.
INFO [2018-04-12 08:13:54] Optimized purity: 0.40
INFO [2018-04-12 08:13:54] Testing local optimum 3/4 at purity 0.30 and total ploidy 1.90...
INFO [2018-04-12 08:13:54] Recalibrating log-ratios...
INFO [2018-04-12 08:13:54] Testing local optimum 3/4 at purity 0.30 and total ploidy 1.90...
INFO [2018-04-12 08:13:54] Recalibrating log-ratios...
INFO [2018-04-12 08:13:54] Testing local optimum 3/4 at purity 0.30 and total ploidy 1.90...
INFO [2018-04-12 08:13:55] Recalibrating log-ratios...
INFO [2018-04-12 08:13:55] Testing local optimum 3/4 at purity 0.30 and total ploidy 1.90...
INFO [2018-04-12 08:13:55] Testing local optimum 4/4 at purity 0.48 and total ploidy 2.00...
INFO [2018-04-12 08:13:56] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:13:57] Optimized purity: 0.65
INFO [2018-04-12 08:13:57] Done.
INFO [2018-04-12 08:13:57] ------------------------------------------------------------
INFO [2018-04-12 08:13:57] ------------------------------------------------------------
INFO [2018-04-12 08:13:57] PureCN 1.8.1
INFO [2018-04-12 08:13:57] ------------------------------------------------------------
INFO [2018-04-12 08:13:57] Loading coverage files...
INFO [2018-04-12 08:13:58] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-04-12 08:13:58] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:13:58] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:13:58] Removing 22 small (< 5bp) targets.
INFO [2018-04-12 08:13:58] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-04-12 08:13:58] No normalDB provided. Provide one for better results.
INFO [2018-04-12 08:13:58] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-04-12 08:13:58] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-04-12 08:13:58] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-04-12 08:13:58] No gc.gene.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-04-12 08:13:58] Loading VCF...
INFO [2018-04-12 08:13:58] LIB-02240e4 is tumor in VCF file.
INFO [2018-04-12 08:13:58] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-04-12 08:13:58] Found 2331 variants in VCF file.
INFO [2018-04-12 08:13:58] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:13:58] Removing 63 variants with AF < 0.030 or AF >= 0.970 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:13:58] Removing 16 low quality variants with BQ < 25.
INFO [2018-04-12 08:13:59] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-04-12 08:13:59] 1.0% of targets contain variants.
INFO [2018-04-12 08:13:59] Setting somatic prior probabilities for dbSNP hits to 0.000500 or to 0.500000 otherwise.
INFO [2018-04-12 08:13:59] Sample sex: ?
INFO [2018-04-12 08:13:59] Segmenting data...
INFO [2018-04-12 08:13:59] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-04-12 08:14:03] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-04-12 08:14:03] Removing 282 variants outside segments.
INFO [2018-04-12 08:14:03] Using 1970 variants.
INFO [2018-04-12 08:14:03] Mean standard deviation of log-ratios: 0.37
INFO [2018-04-12 08:14:03] 2D-grid search of purity and ploidy...
INFO [2018-04-12 08:14:06] Local optima: 0.62/1.9, 0.38/2.2, 0.3/1.9, 0.48/2
INFO [2018-04-12 08:14:06] Testing local optimum 1/4 at purity 0.62 and total ploidy 1.90...
INFO [2018-04-12 08:14:08] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-04-12 08:14:13] Optimized purity: 0.65
INFO [2018-04-12 08:14:13] Testing local optimum 2/4 at purity 0.38 and total ploidy 2.20...
INFO [2018-04-12 08:14:15] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2018-04-12 08:14:20] Optimized purity: 0.40
INFO [2018-04-12 08:14:20] Testing local optimum 3/4 at purity 0.30 and total ploidy 1.90...
INFO [2018-04-12 08:14:21] Recalibrating log-ratios...
INFO [2018-04-12 08:14:21] Testing local optimum 3/4 at purity 0.30 and total ploidy 1.90...
INFO [2018-04-12 08:14:21] Recalibrating log-ratios...
INFO [2018-04-12 08:14:21] Testing local optimum 3/4 at purity 0.30 and total ploidy 1.90...
INFO [2018-04-12 08:14:22] Recalibrating log-ratios...
INFO [2018-04-12 08:14:22] Testing local optimum 3/4 at purity 0.30 and total ploidy 1.90...
INFO [2018-04-12 08:14:22] Testing local optimum 4/4 at purity 0.48 and total ploidy 2.00...
INFO [2018-04-12 08:14:24] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-04-12 08:14:29] Optimized purity: 0.65
INFO [2018-04-12 08:14:29] Done.
INFO [2018-04-12 08:14:29] ------------------------------------------------------------
INFO [2018-04-12 08:14:29] ------------------------------------------------------------
INFO [2018-04-12 08:14:29] PureCN 1.8.1
INFO [2018-04-12 08:14:29] ------------------------------------------------------------
INFO [2018-04-12 08:14:29] Loading coverage files...
INFO [2018-04-12 08:14:30] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-04-12 08:14:30] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:14:30] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:14:30] No Gene column in gc.gene.file. You won't get gene-level calls.
INFO [2018-04-12 08:14:30] No on_target column in gc.gene.file. Recreate this file with IntervalFile.R.
INFO [2018-04-12 08:14:30] No mappability column in gc.gene.file.
INFO [2018-04-12 08:14:30] Removing 15 low/high GC targets.
INFO [2018-04-12 08:14:30] Removing 21 small (< 5bp) targets.
INFO [2018-04-12 08:14:30] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-04-12 08:14:30] No normalDB provided. Provide one for better results.
INFO [2018-04-12 08:14:30] Removing 232 low coverage (< 15.0000X) targets.
INFO [2018-04-12 08:14:30] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2018-04-12 08:14:31] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-04-12 08:14:31] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2018-04-12 08:14:31] Loading VCF...
INFO [2018-04-12 08:14:31] LIB-02240e4 is tumor in VCF file.
INFO [2018-04-12 08:14:31] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-04-12 08:14:31] Found 2331 variants in VCF file.
INFO [2018-04-12 08:14:31] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:14:31] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:14:31] Removing 48 variants with AF < 0.030 or AF >= Inf or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:14:31] Removing 16 low quality variants with BQ < 25.
INFO [2018-04-12 08:14:31] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-04-12 08:14:31] 1.0% of targets contain variants.
INFO [2018-04-12 08:14:31] Removing 2122 variants outside intervals.
INFO [2018-04-12 08:14:31] Found SOMATIC annotation in VCF.
INFO [2018-04-12 08:14:31] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-04-12 08:14:31] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-04-12 08:14:31] Sample sex: ?
INFO [2018-04-12 08:14:31] Segmenting data...
INFO [2018-04-12 08:14:32] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-04-12 08:14:35] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-04-12 08:14:35] Using 123 variants.
INFO [2018-04-12 08:14:35] Mean standard deviation of log-ratios: 0.37
INFO [2018-04-12 08:14:35] 2D-grid search of purity and ploidy...
INFO [2018-04-12 08:14:35] Local optima: 0.65/1.9, 0.33/1.9, 0.38/2.2, 0.27/2.1, 0.29/2, 0.78/2
INFO [2018-04-12 08:14:36] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2018-04-12 08:14:37] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:14:38] Optimized purity: 0.65
INFO [2018-04-12 08:14:38] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-04-12 08:14:38] Recalibrating log-ratios...
INFO [2018-04-12 08:14:38] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-04-12 08:14:38] Recalibrating log-ratios...
INFO [2018-04-12 08:14:38] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-04-12 08:14:39] Recalibrating log-ratios...
INFO [2018-04-12 08:14:39] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-04-12 08:14:39] Testing local optimum 3/6 at purity 0.38 and total ploidy 2.20...
INFO [2018-04-12 08:14:41] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2018-04-12 08:14:42] Optimized purity: 0.40
INFO [2018-04-12 08:14:42] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-04-12 08:14:43] Recalibrating log-ratios...
INFO [2018-04-12 08:14:43] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-04-12 08:14:43] Recalibrating log-ratios...
INFO [2018-04-12 08:14:43] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-04-12 08:14:43] Recalibrating log-ratios...
INFO [2018-04-12 08:14:43] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-04-12 08:14:44] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-04-12 08:14:44] Recalibrating log-ratios...
INFO [2018-04-12 08:14:44] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-04-12 08:14:45] Recalibrating log-ratios...
INFO [2018-04-12 08:14:45] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-04-12 08:14:45] Recalibrating log-ratios...
INFO [2018-04-12 08:14:45] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-04-12 08:14:46] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2018-04-12 08:14:47] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:14:48] Optimized purity: 0.65
INFO [2018-04-12 08:14:48] Done.
INFO [2018-04-12 08:14:48] ------------------------------------------------------------
FATAL [2018-04-12 08:14:48] This function requires gene-level calls. Please add a column 'Gene'
FATAL [2018-04-12 08:14:48] containing gene symbols to the gc.gene.file.
FATAL [2018-04-12 08:14:48]
FATAL [2018-04-12 08:14:48] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:14:48] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:14:50] ------------------------------------------------------------
INFO [2018-04-12 08:14:50] PureCN 1.8.1
INFO [2018-04-12 08:14:50] ------------------------------------------------------------
INFO [2018-04-12 08:14:50] Loading coverage files...
INFO [2018-04-12 08:14:50] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-04-12 08:14:50] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:14:51] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:14:51] Removing 22 small (< 5bp) targets.
INFO [2018-04-12 08:14:51] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-04-12 08:14:51] No normalDB provided. Provide one for better results.
INFO [2018-04-12 08:14:51] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-04-12 08:14:51] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-04-12 08:14:51] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-04-12 08:14:51] No gc.gene.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-04-12 08:14:51] Loading VCF...
FATAL [2018-04-12 08:14:51] Different chromosome names in coverage and VCF.
FATAL [2018-04-12 08:14:51]
FATAL [2018-04-12 08:14:51] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:14:51] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:14:51] ------------------------------------------------------------
INFO [2018-04-12 08:14:51] PureCN 1.8.1
INFO [2018-04-12 08:14:51] ------------------------------------------------------------
INFO [2018-04-12 08:14:51] Loading coverage files...
INFO [2018-04-12 08:14:51] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-04-12 08:14:51] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:14:51] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:14:51] No on_target column in gc.gene.file. Recreate this file with IntervalFile.R.
INFO [2018-04-12 08:14:51] No mappability column in gc.gene.file.
FATAL [2018-04-12 08:14:51] tumor.coverage.file and gc.gene.file do not align.
FATAL [2018-04-12 08:14:51]
FATAL [2018-04-12 08:14:51] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:14:51] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:14:51] ------------------------------------------------------------
INFO [2018-04-12 08:14:51] PureCN 1.8.1
INFO [2018-04-12 08:14:51] ------------------------------------------------------------
INFO [2018-04-12 08:14:51] Loading coverage files...
INFO [2018-04-12 08:14:51] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-04-12 08:14:51] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:14:52] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:14:52] No on_target column in gc.gene.file. Recreate this file with IntervalFile.R.
INFO [2018-04-12 08:14:52] No mappability column in gc.gene.file.
INFO [2018-04-12 08:14:54] Removing 15 low/high GC targets.
INFO [2018-04-12 08:14:54] Removing 21 small (< 5bp) targets.
INFO [2018-04-12 08:14:54] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-04-12 08:14:54] No normalDB provided. Provide one for better results.
INFO [2018-04-12 08:14:54] Removing 232 low coverage (< 15.0000X) targets.
INFO [2018-04-12 08:14:54] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2018-04-12 08:14:54] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-04-12 08:14:54] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2018-04-12 08:14:54] Loading VCF...
INFO [2018-04-12 08:14:54] LIB-02240e4 is tumor in VCF file.
INFO [2018-04-12 08:14:54] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-04-12 08:14:54] Found 2331 variants in VCF file.
INFO [2018-04-12 08:14:54] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:14:54] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:14:54] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:14:54] Removing 16 low quality variants with BQ < 25.
INFO [2018-04-12 08:14:54] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-04-12 08:14:54] 1.0% of targets contain variants.
INFO [2018-04-12 08:14:54] Removing 2122 variants outside intervals.
INFO [2018-04-12 08:14:54] Found SOMATIC annotation in VCF.
INFO [2018-04-12 08:14:54] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-04-12 08:14:54] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-04-12 08:14:54] Sample sex: ?
INFO [2018-04-12 08:14:54] Segmenting data...
INFO [2018-04-12 08:14:54] Setting undo.SD parameter to 1.000000.
INFO [2018-04-12 08:14:58] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-04-12 08:14:58] Using 123 variants.
INFO [2018-04-12 08:14:58] Mean standard deviation of log-ratios: 0.37
INFO [2018-04-12 08:14:58] 2D-grid search of purity and ploidy...
INFO [2018-04-12 08:15:00] Local optima: 0.63/1.9, 0.4/2.2, 0.5/2
INFO [2018-04-12 08:15:00] Testing local optimum 1/3 at purity 0.63 and total ploidy 1.90...
INFO [2018-04-12 08:15:02] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:15:02] Optimized purity: 0.65
INFO [2018-04-12 08:15:02] Testing local optimum 2/3 at purity 0.40 and total ploidy 2.20...
INFO [2018-04-12 08:15:04] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2018-04-12 08:15:05] Optimized purity: 0.40
INFO [2018-04-12 08:15:05] Testing local optimum 3/3 at purity 0.50 and total ploidy 2.00...
INFO [2018-04-12 08:15:07] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:15:07] Optimized purity: 0.65
INFO [2018-04-12 08:15:07] Done.
INFO [2018-04-12 08:15:07] ------------------------------------------------------------
INFO [2018-04-12 08:15:07] ------------------------------------------------------------
INFO [2018-04-12 08:15:07] PureCN 1.8.1
INFO [2018-04-12 08:15:07] ------------------------------------------------------------
INFO [2018-04-12 08:15:07] Loading coverage files...
INFO [2018-04-12 08:15:08] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-04-12 08:15:08] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:15:08] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:15:08] No on_target column in gc.gene.file. Recreate this file with IntervalFile.R.
INFO [2018-04-12 08:15:08] No mappability column in gc.gene.file.
INFO [2018-04-12 08:15:09] Removing 15 low/high GC targets.
INFO [2018-04-12 08:15:09] Removing 21 small (< 5bp) targets.
INFO [2018-04-12 08:15:09] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-04-12 08:15:09] No normalDB provided. Provide one for better results.
INFO [2018-04-12 08:15:09] Removing 232 low coverage (< 15.0000X) targets.
INFO [2018-04-12 08:15:09] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2018-04-12 08:15:09] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-04-12 08:15:09] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2018-04-12 08:15:09] Loading VCF...
INFO [2018-04-12 08:15:09] LIB-02240e4 is tumor in VCF file.
INFO [2018-04-12 08:15:09] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-04-12 08:15:09] Found 2331 variants in VCF file.
INFO [2018-04-12 08:15:09] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:15:09] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:15:10] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:15:10] Removing 16 low quality variants with BQ < 25.
INFO [2018-04-12 08:15:10] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-04-12 08:15:10] 1.0% of targets contain variants.
INFO [2018-04-12 08:15:10] Removing 2122 variants outside intervals.
INFO [2018-04-12 08:15:10] Found SOMATIC annotation in VCF.
INFO [2018-04-12 08:15:10] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-04-12 08:15:10] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-04-12 08:15:10] Sample sex: ?
INFO [2018-04-12 08:15:10] Segmenting data...
INFO [2018-04-12 08:15:10] Setting undo.SD parameter to 1.000000.
INFO [2018-04-12 08:15:13] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-04-12 08:15:14] Using 123 variants.
INFO [2018-04-12 08:15:14] Mean standard deviation of log-ratios: 0.37
INFO [2018-04-12 08:15:14] 2D-grid search of purity and ploidy...
INFO [2018-04-12 08:15:15] Local optima: 0.63/1.9, 0.4/2.2, 0.5/2
INFO [2018-04-12 08:15:16] Testing local optimum 1/3 at purity 0.63 and total ploidy 1.90...
INFO [2018-04-12 08:15:17] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:15:18] Optimized purity: 0.65
INFO [2018-04-12 08:15:18] Testing local optimum 2/3 at purity 0.40 and total ploidy 2.20...
INFO [2018-04-12 08:15:19] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2018-04-12 08:15:20] Optimized purity: 0.40
INFO [2018-04-12 08:15:20] Testing local optimum 3/3 at purity 0.50 and total ploidy 2.00...
INFO [2018-04-12 08:15:22] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:15:23] Optimized purity: 0.65
INFO [2018-04-12 08:15:23] Done.
INFO [2018-04-12 08:15:23] ------------------------------------------------------------
FATAL [2018-04-12 08:15:23] chr1 not valid chromosome name(s). Valid names are:
FATAL [2018-04-12 08:15:23] 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22
FATAL [2018-04-12 08:15:23]
FATAL [2018-04-12 08:15:23] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:15:23] parameters (PureCN 1.8.1).
INFO [2018-04-12 08:15:23] ------------------------------------------------------------
INFO [2018-04-12 08:15:23] PureCN 1.8.1
INFO [2018-04-12 08:15:23] ------------------------------------------------------------
INFO [2018-04-12 08:15:23] Loading coverage files...
WARN [2018-04-12 08:15:24] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
WARN [2018-04-12 08:15:24] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:15:24] Allosome coverage missing, cannot determine sex.
INFO [2018-04-12 08:15:24] Using 10039 intervals (10039 on-target, 0 off-target).
INFO [2018-04-12 08:15:24] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-04-12 08:15:24] Loading VCF...
INFO [2018-04-12 08:15:24] LIB-02240e4 is tumor in VCF file.
INFO [2018-04-12 08:15:24] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-04-12 08:15:24] Found 2331 variants in VCF file.
INFO [2018-04-12 08:15:24] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-04-12 08:15:24] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-04-12 08:15:24] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-04-12 08:15:24] Removing 16 low quality variants with BQ < 25.
INFO [2018-04-12 08:15:24] Total size of targeted genomic region: 1.71Mb (2.70Mb with 50bp padding).
INFO [2018-04-12 08:15:25] 1.0% of targets contain variants.
INFO [2018-04-12 08:15:25] Removing 2120 variants outside intervals.
INFO [2018-04-12 08:15:25] Found SOMATIC annotation in VCF.
INFO [2018-04-12 08:15:25] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-04-12 08:15:25] Found SOMATIC annotation in VCF. Setting mapping bias to 0.976.
INFO [2018-04-12 08:15:25] Sample sex: ?
INFO [2018-04-12 08:15:25] Segmenting data...
WARN [2018-04-12 08:15:25] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
INFO [2018-04-12 08:15:25] Setting undo.SD parameter to 0.000000.
Setting multi-figure configuration
INFO [2018-04-12 08:15:28] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-04-12 08:15:28] Using 125 variants.
INFO [2018-04-12 08:15:28] Mean standard deviation of log-ratios: 0.40
INFO [2018-04-12 08:15:28] 2D-grid search of purity and ploidy...
INFO [2018-04-12 08:15:30] Local optima: 0.67/1.8, 0.4/2.2, 0.5/2
INFO [2018-04-12 08:15:30] Testing local optimum 1/3 at purity 0.67 and total ploidy 1.80...
INFO [2018-04-12 08:15:31] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:15:32] Optimized purity: 0.65
INFO [2018-04-12 08:15:32] Testing local optimum 2/3 at purity 0.40 and total ploidy 2.20...
INFO [2018-04-12 08:15:33] Fitting variants for purity 0.40, tumor ploidy 2.33 and contamination 0.01.
INFO [2018-04-12 08:15:34] Optimized purity: 0.40
INFO [2018-04-12 08:15:34] Testing local optimum 3/3 at purity 0.50 and total ploidy 2.00...
INFO [2018-04-12 08:15:35] Fitting variants for purity 0.65, tumor ploidy 1.73 and contamination 0.01.
INFO [2018-04-12 08:15:36] Optimized purity: 0.65
INFO [2018-04-12 08:15:36] Done.
INFO [2018-04-12 08:15:36] ------------------------------------------------------------
WARN [2018-04-12 08:15:36] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:15:36] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
FATAL [2018-04-12 08:15:57] Segmentation file expected with colnames ID, chrom, loc.start, loc.end,
FATAL [2018-04-12 08:15:57] num.mark, seg.mean
FATAL [2018-04-12 08:15:57]
FATAL [2018-04-12 08:15:57] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:15:57] parameters (PureCN 1.8.1).
FATAL [2018-04-12 08:15:57] seg.file contains multiple samples and sampleid missing.
FATAL [2018-04-12 08:15:57]
FATAL [2018-04-12 08:15:57] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:15:57] parameters (PureCN 1.8.1).
FATAL [2018-04-12 08:15:57] seg.file contains multiple samples and sampleid does not match any.
FATAL [2018-04-12 08:15:57]
FATAL [2018-04-12 08:15:57] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:15:57] parameters (PureCN 1.8.1).
WARN [2018-04-12 08:15:58] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:15:58] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
WARN [2018-04-12 08:16:16] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:16:16] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:16:17] No normalDB provided. Provide one for better results.
Setting multi-figure configuration
WARN [2018-04-12 08:16:42] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:16:42] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:16:44] Provided sampleid (Sample.1) does not match Sample1 found in segmentation.
Setting multi-figure configuration
WARN [2018-04-12 08:17:22] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:17:22] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:17:22] No normalDB provided. Provide one for better results.
WARN [2018-04-12 08:17:22] Sampleid looks like a normal in VCF, not like a tumor.
Setting multi-figure configuration
WARN [2018-04-12 08:17:52] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:17:52] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:17:52] No normalDB provided. Provide one for better results.
WARN [2018-04-12 08:18:05] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:18:05] Allosome coverage missing, cannot determine sex.
FATAL [2018-04-12 08:18:05] normalDB not a valid normalDB object. Use createNormalDatabase to
FATAL [2018-04-12 08:18:05] create one.
FATAL [2018-04-12 08:18:05]
FATAL [2018-04-12 08:18:05] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:18:05] parameters (PureCN 1.8.1).
WARN [2018-04-12 08:18:06] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:18:06] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:18:08] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:18:09] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
WARN [2018-04-12 08:18:41] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:18:41] Allosome coverage missing, cannot determine sex.
FATAL [2018-04-12 08:18:41] normalDB appears to be empty.
FATAL [2018-04-12 08:18:41]
FATAL [2018-04-12 08:18:41] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:18:41] parameters (PureCN 1.8.1).
WARN [2018-04-12 08:18:42] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:18:42] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:19:12] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2018-04-12 08:19:12] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:19:12] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:19:12] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2018-04-12 08:19:20] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:19:20] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
WARN [2018-04-12 08:19:46] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:19:46] Allosome coverage missing, cannot determine sex.
WARN [2018-04-12 08:19:46] Intervals in coverage and gc.gene.file have conflicting on/off-target annotation.
WARN [2018-04-12 08:19:47] No normalDB provided. Provide one for better results.
Setting multi-figure configuration
FATAL [2018-04-12 08:20:01] min.normals must be >=2.
FATAL [2018-04-12 08:20:01]
FATAL [2018-04-12 08:20:01] This is most likely a user error due to invalid input data or
FATAL [2018-04-12 08:20:01] parameters (PureCN 1.8.1).
RUNIT TEST PROTOCOL -- Thu Apr 12 08:20:02 2018
***********************************************
Number of test functions: 22
Number of errors: 0
Number of failures: 0
1 Test Suite :
PureCN RUnit Tests - 22 test functions, 0 errors, 0 failures
Number of test functions: 22
Number of errors: 0
Number of failures: 0
There were 18 warnings (use warnings() to see them)
>
> proc.time()
user system elapsed
522.095 5.917 534.799