Back to Multiple platform build/check report for BioC 3.22:   simplified   long
ABCD[E]FGHIJKLMNOPQRSTUVWXYZ

This page was generated on 2025-11-20 12:04 -0500 (Thu, 20 Nov 2025).

HostnameOSArch (*)R versionInstalled pkgs
nebbiolo2Linux (Ubuntu 24.04.3 LTS)x86_644.5.2 (2025-10-31) -- "[Not] Part in a Rumble" 4615
merida1macOS 12.7.6 Montereyx86_644.5.2 Patched (2025-11-05 r88990) -- "[Not] Part in a Rumble" 4610
kjohnson1macOS 13.7.5 Venturaarm644.5.2 Patched (2025-11-04 r88984) -- "[Not] Part in a Rumble" 4598
taishanLinux (openEuler 24.03 LTS)aarch644.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/2361HostnameOS / ArchINSTALLBUILDCHECKBUILD BIN
ELViS 1.2.0  (landing page)
Jin-Young Lee
Snapshot Date: 2025-11-17 13:45 -0500 (Mon, 17 Nov 2025)
git_url: https://git.bioconductor.org/packages/ELViS
git_branch: RELEASE_3_22
git_last_commit: c2a03f6
git_last_commit_date: 2025-10-29 11:34:26 -0500 (Wed, 29 Oct 2025)
nebbiolo2Linux (Ubuntu 24.04.3 LTS) / x86_64  OK    OK    OK  UNNEEDED, same version is already published
merida1macOS 12.7.6 Monterey / x86_64  OK    OK    ERROR    OK  
kjohnson1macOS 13.7.5 Ventura / arm64  OK    OK    ERROR    OK  
taishanLinux (openEuler 24.03 LTS) / aarch64  OK    OK    OK  


CHECK results for ELViS on merida1

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.

raw results


Summary

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

Command output

##############################################################################
##############################################################################
###
### 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.


Installation output

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)

Tests output

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

Example timings

ELViS.Rcheck/ELViS-Ex.timings

nameusersystemelapsed
coord_to_grng0.1930.0040.210
coord_to_lst0.0030.0010.004
depth_hist2.8970.0623.141
filt_samples1.3700.0251.474
gene_cn_heatmaps29.126 0.41231.182
get_depth_matrix1.9680.9399.374
get_new_baseline0.4260.0150.443
integrative_heatmap85.999 2.57993.473
norm_fun0.0010.0010.003
plot_pileUp_multisample5.3710.2776.454
run_ELViS161.339 1.882176.420