This page was generated on 2019-04-16 11:54:08 -0400 (Tue, 16 Apr 2019).
DNAcopy 1.56.0 Venkatraman E. Seshan
Snapshot Date: 2019-04-15 17:01:12 -0400 (Mon, 15 Apr 2019) |
URL: https://git.bioconductor.org/packages/DNAcopy |
Branch: RELEASE_3_8 |
Last Commit: e521826 |
Last Changed Date: 2018-10-30 11:41:41 -0400 (Tue, 30 Oct 2018) |
| malbec1 | Linux (Ubuntu 16.04.6 LTS) / x86_64 | OK | OK | OK | | |
merida1 | OS X 10.11.6 El Capitan / x86_64 | OK | OK | [ OK ] | OK | |
R version 3.5.3 (2019-03-11) -- "Great Truth"
Copyright (C) 2019 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.
> ######################################################################
> # Type: Redundancy test
> # Created by: Henrik Bengtsson <hb@stat.berkeley.edu>
> # Created on: 2009-06-10
> ######################################################################
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Startup
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> library("DNAcopy")
>
> # Record current random seed
> sample(1) # Assert that a random seed exists
[1] 1
> oldSeed <- .Random.seed
> # Alway use the same random seed
> set.seed(0xbeef)
>
> # Tolerance (maybe decrease?)
> tol <- .Machine$double.eps^0.5
>
> print(sessionInfo())
R version 3.5.3 (2019-03-11)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: OS X El Capitan 10.11.6
Matrix products: default
BLAS: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRlapack.dylib
locale:
[1] C/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] DNAcopy_1.56.0
loaded via a namespace (and not attached):
[1] compiler_3.5.3
>
>
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Simulating copy-number data
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Number of loci
> J <- 1000
>
> x <- sort(runif(J, min=0, max=1000))
> w <- runif(J)
> mu <- double(J)
> jj <- (200 <= x & x < 300)
> mu[jj] <- mu[jj] + 1
> jj <- (650 <= x & x < 800)
> mu[jj] <- mu[jj] - 1
> w[jj] <- 0.001
> eps <- rnorm(J, sd=1/2)
> y <- mu + eps
>
>
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Setting up a raw CNA object
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> cnR <- CNA(
+ genomdat = y,
+ chrom = rep(1, times=J),
+ maploc = x,
+ data.type = "logratio",
+ sampleid = "SampleA"
+ )
> print(cnR)
Number of Samples 1
Number of Probes 1000
Data Type logratio
>
>
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Test: Non-weighted segmentation
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> t <- system.time({
+ fitR <- segment(cnR, verbose=1)
+ })
Analyzing: SampleA
> cat("Processing time:\n")
Processing time:
> print(t)
user system elapsed
0.046 0.002 0.048
> print(fitR)
Call:
segment(x = cnR, verbose = 1)
ID chrom loc.start loc.end num.mark seg.mean
1 SampleA 1 1.368577 199.0840 209 0.0256
2 SampleA 1 201.604291 301.0669 105 1.0099
3 SampleA 1 303.775112 647.4270 337 -0.0084
4 SampleA 1 650.741212 798.9718 138 -0.9792
5 SampleA 1 800.302447 999.3290 211 -0.0289
>
> # Expected results
> # These were obtained by dput(fitR$output) using DNAcopy v1.19.0
> truth <- structure(list(ID = c("SampleA", "SampleA", "SampleA", "SampleA",
+ "SampleA"), chrom = c(1, 1, 1, 1, 1), loc.start = c(1.36857712641358,
+ 201.604291098192, 303.775111911818, 650.741211604327, 800.302447052673
+ ), loc.end = c(199.083976913244, 301.066882908344, 647.42697100155,
+ 798.971758922562, 999.329038895667), num.mark = c(209, 105, 337,
+ 138, 211), seg.mean = c(0.0256, 1.0099, -0.0084, -0.9792, -0.0289
+ )), .Names = c("ID", "chrom", "loc.start", "loc.end", "num.mark",
+ "seg.mean"), row.names = c(NA, -5L), class = "data.frame")
>
> stopifnot(all.equal(fitR$output, truth, tolerance=tol))
>
>
>
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Test: Weighted segmentation
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> t <- system.time({
+ fitR <- segment(cnR, weights=w, verbose=1)
+ })
Analyzing: SampleA
> cat("Processing time:\n")
Processing time:
> print(t)
user system elapsed
0.022 0.000 0.023
> print(fitR)
Call:
segment(x = cnR, weights = w, verbose = 1)
ID chrom loc.start loc.end num.mark seg.mean
1 SampleA 1 1.368577 199.0840 209 0.0259
2 SampleA 1 201.604291 301.0669 105 1.0004
3 SampleA 1 303.775112 999.3290 686 -0.0233
>
> # Expected results
> # These were obtained by dput(fitR$output) using DNAcopy v1.19.0
> truth <- structure(list(ID = c("SampleA", "SampleA", "SampleA"), chrom = c(1,
+ 1, 1), loc.start = c(1.36857712641358, 201.604291098192, 303.775111911818
+ ), loc.end = c(199.083976913244, 301.066882908344, 999.329038895667
+ ), num.mark = c(209, 105, 686), seg.mean = c(0.0259, 1.0004,
+ -0.0233)), .Names = c("ID", "chrom", "loc.start", "loc.end",
+ "num.mark", "seg.mean"), row.names = c(NA, -3L), class = "data.frame")
>
> stopifnot(all.equal(fitR$output, truth, tolerance=tol))
>
>
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Cleanup
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Reset to previous random seed
> .Random.seed <- oldSeed
>
> print(sessionInfo())
R version 3.5.3 (2019-03-11)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: OS X El Capitan 10.11.6
Matrix products: default
BLAS: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRlapack.dylib
locale:
[1] C/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] DNAcopy_1.56.0
loaded via a namespace (and not attached):
[1] compiler_3.5.3
>
>
> ######################################################################
> # HISTORY
> # 2009-06-10
> # o ROBUSTNESS: Added this test to assert that DNAcopy v1.19.2 and
> # newer will numerically give the same results as DNAcopy v1.19.0.
> # This test is ran each time with R CMD check.
> # o Created.
> ######################################################################
>
> proc.time()
user system elapsed
0.442 0.087 0.503