| Back to Multiple platform build/check report for BioC 3.22: simplified long |
|
This page was generated on 2025-11-20 12:04 -0500 (Thu, 20 Nov 2025).
| Hostname | OS | Arch (*) | R version | Installed pkgs |
|---|---|---|---|---|
| nebbiolo2 | Linux (Ubuntu 24.04.3 LTS) | x86_64 | 4.5.2 (2025-10-31) -- "[Not] Part in a Rumble" | 4615 |
| merida1 | macOS 12.7.6 Monterey | x86_64 | 4.5.2 Patched (2025-11-05 r88990) -- "[Not] Part in a Rumble" | 4610 |
| kjohnson1 | macOS 13.7.5 Ventura | arm64 | 4.5.2 Patched (2025-11-04 r88984) -- "[Not] Part in a Rumble" | 4598 |
| taishan | Linux (openEuler 24.03 LTS) | aarch64 | 4.5.0 (2025-04-11) -- "How About a Twenty-Six" | 4668 |
| Click on any hostname to see more info about the system (e.g. compilers) (*) as reported by 'uname -p', except on Windows and Mac OS X | ||||
| Package 662/2361 | Hostname | OS / Arch | INSTALL | BUILD | CHECK | BUILD BIN | ||||||||
| ELViS 1.2.0 (landing page) Jin-Young Lee
| nebbiolo2 | Linux (Ubuntu 24.04.3 LTS) / x86_64 | OK | OK | OK | |||||||||
| merida1 | macOS 12.7.6 Monterey / x86_64 | OK | OK | ERROR | OK | |||||||||
| kjohnson1 | macOS 13.7.5 Ventura / arm64 | OK | OK | ERROR | OK | |||||||||
| taishan | Linux (openEuler 24.03 LTS) / aarch64 | OK | OK | OK | ||||||||||
|
To the developers/maintainers of the ELViS package: - Allow up to 24 hours (and sometimes 48 hours) for your latest push to git@git.bioconductor.org:packages/ELViS.git to reflect on this report. See Troubleshooting Build Report for more information. - Use the following Renviron settings to reproduce errors and warnings. - If 'R CMD check' started to fail recently on the Linux builder(s) over a missing dependency, add the missing dependency to 'Suggests:' in your DESCRIPTION file. See Renviron.bioc for more information. |
| Package: ELViS |
| Version: 1.2.0 |
| Command: /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:ELViS.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings ELViS_1.2.0.tar.gz |
| StartedAt: 2025-11-18 07:05:30 -0500 (Tue, 18 Nov 2025) |
| EndedAt: 2025-11-18 07:26:01 -0500 (Tue, 18 Nov 2025) |
| EllapsedTime: 1231.1 seconds |
| RetCode: 1 |
| Status: ERROR |
| CheckDir: ELViS.Rcheck |
| Warnings: NA |
##############################################################################
##############################################################################
###
### Running command:
###
### /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:ELViS.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings ELViS_1.2.0.tar.gz
###
##############################################################################
##############################################################################
* using log directory ‘/Users/biocbuild/bbs-3.22-bioc/meat/ELViS.Rcheck’
* using R version 4.5.2 Patched (2025-11-05 r88990)
* using platform: x86_64-apple-darwin20
* R was compiled by
Apple clang version 14.0.0 (clang-1400.0.29.202)
GNU Fortran (GCC) 14.2.0
* running under: macOS Monterey 12.7.6
* using session charset: UTF-8
* using option ‘--no-vignettes’
* checking for file ‘ELViS/DESCRIPTION’ ... OK
* this is package ‘ELViS’ version ‘1.2.0’
* package encoding: UTF-8
* checking package namespace information ... OK
* checking package dependencies ... INFO
Imports includes 21 non-default packages.
Importing from so many packages makes the package vulnerable to any of
them becoming unavailable. Move as many as possible to Suggests and
use conditionally.
* checking if this is a source package ... OK
* checking if there is a namespace ... OK
* checking for hidden files and directories ... OK
* checking for portable file names ... OK
* checking for sufficient/correct file permissions ... OK
* checking whether package ‘ELViS’ can be installed ... OK
* checking installed package size ... OK
* checking package directory ... OK
* checking ‘build’ directory ... OK
* checking DESCRIPTION meta-information ... OK
* checking top-level files ... OK
* checking for left-over files ... OK
* checking index information ... OK
* checking package subdirectories ... OK
* checking code files for non-ASCII characters ... OK
* checking R files for syntax errors ... OK
* checking whether the package can be loaded ... OK
* checking whether the package can be loaded with stated dependencies ... OK
* checking whether the package can be unloaded cleanly ... OK
* checking whether the namespace can be loaded with stated dependencies ... OK
* checking whether the namespace can be unloaded cleanly ... OK
* checking whether startup messages can be suppressed ... OK
* checking dependencies in R code ... OK
* checking S3 generic/method consistency ... OK
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... OK
* checking Rd files ... OK
* checking Rd metadata ... OK
* checking Rd cross-references ... OK
* checking for missing documentation entries ... OK
* checking for code/documentation mismatches ... OK
* checking Rd \usage sections ... OK
* checking Rd contents ... OK
* checking for unstated dependencies in examples ... OK
* checking contents of ‘data’ directory ... OK
* checking data for non-ASCII characters ... OK
* checking data for ASCII and uncompressed saves ... OK
* checking files in ‘vignettes’ ... OK
* checking examples ... OK
Examples with CPU (user + system) or elapsed time > 5s
user system elapsed
run_ELViS 161.339 1.882 176.420
integrative_heatmap 85.999 2.579 93.473
gene_cn_heatmaps 29.126 0.412 31.182
plot_pileUp_multisample 5.371 0.277 6.454
get_depth_matrix 1.968 0.939 9.374
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
Running ‘testthat.R’
ERROR
Running the tests in ‘tests/testthat.R’ failed.
Last 13 lines of output:
[ FAIL 1 | WARN 2 | SKIP 0 | PASS 83 ]
══ Failed tests ════════════════════════════════════════════════════════════════
── Error ('test-Process_Bam_Test.R:210:9'): (code run outside of `test_that()`) ──
Error in `eval(code, test_env)`: object 'samtools_to_install' not found
Backtrace:
▆
1. └─reticulate::conda_create(...) at test-Process_Bam_Test.R:210:9
2. └─base::grepl("^python", packages)
3. └─base::is.factor(x)
[ FAIL 1 | WARN 2 | SKIP 0 | PASS 83 ]
Error:
! Test failures.
Execution halted
* checking for unstated dependencies in vignettes ... OK
* checking package vignettes ... OK
* checking running R code from vignettes ... SKIPPED
* checking re-building of vignette outputs ... SKIPPED
* checking PDF version of manual ... OK
* DONE
Status: 1 ERROR
See
‘/Users/biocbuild/bbs-3.22-bioc/meat/ELViS.Rcheck/00check.log’
for details.
ELViS.Rcheck/00install.out
############################################################################## ############################################################################## ### ### Running command: ### ### /Library/Frameworks/R.framework/Resources/bin/R CMD INSTALL ELViS ### ############################################################################## ############################################################################## * installing to library ‘/Library/Frameworks/R.framework/Versions/4.5-x86_64/Resources/library’ * installing *source* package ‘ELViS’ ... ** this is package ‘ELViS’ version ‘1.2.0’ ** using staged installation ** R ** data ** inst ** byte-compile and prepare package for lazy loading ** help *** installing help indices ** building package indices ** installing vignettes ** testing if installed package can be loaded from temporary location ** testing if installed package can be loaded from final location ** testing if installed package keeps a record of temporary installation path * DONE (ELViS)
ELViS.Rcheck/tests/testthat.Rout.fail
R version 4.5.2 Patched (2025-11-05 r88990) -- "[Not] Part in a Rumble"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20
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.
> # This file is part of the standard setup for testthat.
> # It is recommended that you do not modify it.
> #
> # Where should you do additional test configuration?
> # Learn more about the roles of various files in:
> # * https://r-pkgs.org/testing-design.html#sec-tests-files-overview
> # * https://testthat.r-lib.org/articles/special-files.html
>
> library(testthat)
> library(ELViS)
>
> test_check("ELViS")
ELViS run starts.
1
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
1| done
2
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
2| done
3
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
3| done
4
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
4| done
5
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
5| done
6
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
6| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
1| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
2| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
3| done
4| done
5| done
6| done
1
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
1| done
2
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
2| done
3
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
3| done
4
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
4| done
5
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
5| done
6
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
6| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
1| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
2| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
3| done
4| done
5| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
6| done
Normalization done.
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
1| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
2| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
3| done
4| done
5| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
6| done
1
1
2
2
3
3
4
4
5
5
6
6
Segmentation done.
ELViS run starts.
1
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
1| done
2
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
2| done
3
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
3| done
4
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
4| done
5
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
5| done
6
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
6| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
1| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
2| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
3| done
4| done
5| done
6| done
1
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
1| done
2
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
2| done
3
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
3| done
4
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
4| done
5
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
5| done
6
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
6| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
1| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
2| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
3| done
4| done
5| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
6| done
Normalization done.
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
1| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
2| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
3| done
4| done
5| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
6| done
1
1
2
2
3
3
4
4
5
5
6
6
Segmentation done.
1
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
1| done
2
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
2| done
3
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
3| done
4
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
4| done
5
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
5| done
6
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
6| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
1| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
2| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
! Argument ncluster was not provided. Selecting values with BIC
i BIC-selected number of class : ncluster = 2
BIC-selected number of segment : nseg = 2
3| done
4| done
5| done
6| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
1| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
2| done
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using ncluster = 2L
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
! Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running Segmentation/Clustering algorithm -----------------------------------
i Running Segmentation/Clustering with lmin = 5, Kmax = 3 and ncluster = 2L
> Calculating initial segmentation without clustering
v Initial segmentation with no cluster calculated.
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
-- Segmentation/Clustering with ncluster = 2
> Calculating initial segmentation without clustering
> Calculating initial segmentation without clustering
> Segmentation-Clustering for ncluster = 2 and nseg = 2/3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
v Segmentation-Clustering successful for ncluster = 2 and nseg = 2:3
> Segmentation-Clustering for ncluster = 2 and nseg = 3/3
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
v Smoothing successful for ncluster = 2
> Smoothing likelihood for ncluster = 2. This step can be lengthy.
> Calculating initial segmentation without clustering
v Segmentation/Clustering with ncluster = 2 successfully calculated.
BIC selected : nseg = 2
> Calculating initial segmentation without clustering
-- Segmentation/Clustering results ---------------------------------------------
v Best segmentation/clustering estimated with 2 clusters and 2 segments according to BIC
> Number of cluster should preferentially be selected according to biological
knowledge. Exploring the BIC plot with plot_BIC() can also provide advice to
select the number of clusters.
> Once number of clusters is selected, the number of segment cab be selected
according to BIC.
> Results of the segmentation/clustering may further be explored with plot()
and segmap()
3| done
4| done
5| done
6| done
1
1
2
2
3
3
4
4
5
5
6
6
-- Checking arguments ----------------------------------------------------------
v Segmentation with seg.var = c("z", "y")
v Using lmin = 300
v Using Kmax = 10
v Using scale.variable = FALSE
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("z", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "z"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
v Using subsample_by = 60
v subsampling by 60
v Adjusting lmin to subsampling.
Dividing lmin by 60, with a minimum of 5
> After subsampling, lmin = 5.
Corresponding to lmin = 300 on the original time scale
v Adjusting Kmax so that lmin*Kmax < nrow(x). Kmax = 3
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 3
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 3
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 2 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
n_cycle : 1
N_alt_ori
n_cycle : 1
N_alt_ori
The path to samtools not provided.
Default samtools is used : /Users/biocbuild/Library/r-miniconda/envs/env_samtools_auto/bin/samtools
The path to samtools not provided.
Default samtools is used : /Users/biocbuild/Library/r-miniconda/envs/env_samtools_auto/bin/samtools
+ /Users/biocbuild/Library/r-miniconda/bin/conda create --yes --name env_samtools_1.21 'python=3.12' 'samtools=1.21' --quiet -c conda-forge -c bioconda
Channels:
- conda-forge
- bioconda
- defaults
Platform: osx-64
Collecting package metadata (repodata.json): ...working... Retrying (Retry(total=2, connect=None, read=None, redirect=None, status=None)) after connection broken by 'SSLError(SSLError(1, '[SSL: TLSV1_ALERT_DECODE_ERROR] tlsv1 alert decode error (_ssl.c:1000)'))': /pkgs/main/noarch/repodata.json.zst
done
Solving environment: ...working... done
## Package Plan ##
environment location: /Users/biocbuild/Library/r-miniconda/envs/env_samtools_1.21
added / updated specs:
- python=3.12
- samtools=1.21
The following packages will be downloaded:
package | build
---------------------------|-----------------
samtools-1.21 | ha21ef43_1 471 KB bioconda
------------------------------------------------------------
Total: 471 KB
The following NEW packages will be INSTALLED:
bzip2 conda-forge/osx-64::bzip2-1.0.8-h500dc9f_8
c-ares conda-forge/osx-64::c-ares-1.34.5-hf13058a_0
ca-certificates conda-forge/noarch::ca-certificates-2025.11.12-hbd8a1cb_0
htslib bioconda/osx-64::htslib-1.22.1-h9f635df_0
krb5 conda-forge/osx-64::krb5-1.21.3-h37d8d59_0
libcurl conda-forge/osx-64::libcurl-8.17.0-h7dd4100_0
libcxx conda-forge/osx-64::libcxx-21.1.5-h3d58e20_0
libdeflate conda-forge/osx-64::libdeflate-1.25-h517ebb2_0
libedit conda-forge/osx-64::libedit-3.1.20250104-pl5321ha958ccf_0
libev conda-forge/osx-64::libev-4.33-h10d778d_2
libexpat conda-forge/osx-64::libexpat-2.7.1-h21dd04a_0
libffi conda-forge/osx-64::libffi-3.5.2-h750e83c_0
liblzma conda-forge/osx-64::liblzma-5.8.1-hd471939_2
libnghttp2 conda-forge/osx-64::libnghttp2-1.67.0-h3338091_0
libsqlite conda-forge/osx-64::libsqlite-3.51.0-h86bffb9_0
libssh2 conda-forge/osx-64::libssh2-1.11.1-hed3591d_0
libzlib conda-forge/osx-64::libzlib-1.3.1-hd23fc13_2
ncurses conda-forge/osx-64::ncurses-6.5-h0622a9a_3
openssl conda-forge/osx-64::openssl-3.6.0-h230baf5_0
pip conda-forge/noarch::pip-25.3-pyh8b19718_0
python conda-forge/osx-64::python-3.12.12-h74c2667_1_cpython
readline conda-forge/osx-64::readline-8.2-h7cca4af_2
samtools bioconda/osx-64::samtools-1.21-ha21ef43_1
setuptools conda-forge/noarch::setuptools-80.9.0-pyhff2d567_0
tk conda-forge/osx-64::tk-8.6.13-hf689a15_3
tzdata conda-forge/noarch::tzdata-2025b-h78e105d_0
wheel conda-forge/noarch::wheel-0.45.1-pyhd8ed1ab_1
zstd conda-forge/osx-64::zstd-1.5.7-h8210216_2
Preparing transaction: ...working... done
Verifying transaction: ...working... done
Executing transaction: ...working... done
The path to samtools not provided.
Default samtools is used : /Users/biocbuild/Library/r-miniconda/envs/env_samtools_1.21/bin/samtools
The path to samtools not provided.
Default samtools is used : /Users/biocbuild/Library/r-miniconda/envs/env_samtools_1.21/bin/samtools
The path to samtools not provided.
Default samtools is used : /Users/biocbuild/Library/r-miniconda/envs/env_samtools_1.21/bin/samtools
The path to samtools not provided.
Default samtools is used : /Users/biocbuild/Library/r-miniconda/envs/env_samtools_1.21/bin/samtools
Saving _problems/test-Process_Bam_Test-214.R
-- Checking arguments ----------------------------------------------------------
! Argument seg.var missing
taking default value seg.var = c("x","y")
v Segmentation with seg.var = c("x", "y")
v Using lmin = 5
v Using Kmax = 2
! Argument scale.variable missing
Taking default value scale.variable = FALSE for segmentation().
i Argument diag.var was not provided
Taking default seg.var as diagnostic variables diag.var.
Setting diag.var = c("x", "y")
i Argument order.var was not provided
Taking default diag.var[1] as ordering variable order.var.
Setting order.var = "x"
-- Preparing and checking data -------------------------------------------------
-- Subsampling --
! Subsampling automatically activated. To disable it, provide subsample = FALSE
i Argument subsample_over was not provided
Taking default value for segmentation()
Setting subsample_over = 10000
v nrow(x) < subsample_over, no subsample needed
-- Scaling and final data check --
v No variable rescaling.
To activate, use scale.variable = TRUE
v Data have no repetition of nearly-identical values larger than lmin
-- Running segmentation algorithm ----------------------------------------------
i Running segmentation with lmin = 5 and Kmax = 2
> Calculating cost matrix
v Cost matrix calculated
> Calculating cost matrix
> Dynamic Programming
v Optimal segmentation calculated for all number of segments <= 2
> Dynamic Programming
> Calculating segment statistics
v Best segmentation estimated with 1 segments, according to Lavielle's criterium
Other number of segments may be selected by looking for likelihood breaks with
plot_likelihood()
Results of the segmentation may be explored with plot() and segmap()
[ FAIL 1 | WARN 2 | SKIP 0 | PASS 83 ]
══ Failed tests ════════════════════════════════════════════════════════════════
── Error ('test-Process_Bam_Test.R:210:9'): (code run outside of `test_that()`) ──
Error in `eval(code, test_env)`: object 'samtools_to_install' not found
Backtrace:
▆
1. └─reticulate::conda_create(...) at test-Process_Bam_Test.R:210:9
2. └─base::grepl("^python", packages)
3. └─base::is.factor(x)
[ FAIL 1 | WARN 2 | SKIP 0 | PASS 83 ]
Error:
! Test failures.
Execution halted
ELViS.Rcheck/ELViS-Ex.timings
| name | user | system | elapsed | |
| coord_to_grng | 0.193 | 0.004 | 0.210 | |
| coord_to_lst | 0.003 | 0.001 | 0.004 | |
| depth_hist | 2.897 | 0.062 | 3.141 | |
| filt_samples | 1.370 | 0.025 | 1.474 | |
| gene_cn_heatmaps | 29.126 | 0.412 | 31.182 | |
| get_depth_matrix | 1.968 | 0.939 | 9.374 | |
| get_new_baseline | 0.426 | 0.015 | 0.443 | |
| integrative_heatmap | 85.999 | 2.579 | 93.473 | |
| norm_fun | 0.001 | 0.001 | 0.003 | |
| plot_pileUp_multisample | 5.371 | 0.277 | 6.454 | |
| run_ELViS | 161.339 | 1.882 | 176.420 | |