==48707== Memcheck, a memory error detector ==48707== Copyright (C) 2002-2017, and GNU GPL'd, by Julian Seward et al. ==48707== Using Valgrind-3.15.0 and LibVEX; rerun with -h for copyright info ==48707== Command: /data/blackswan/ripley/R/R-devel-vg/bin/exec/R --vanilla ==48707== R Under development (unstable) (2020-02-08 r77783) -- "Unsuffered Consequences" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale 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. > pkgname <- "cghseg" > source(file.path(R.home("share"), "R", "examples-header.R")) > options(warn = 1) > library('cghseg') Loading required package: parallel > > base::assign(".oldSearch", base::search(), pos = 'CheckExEnv') > base::assign(".old_wd", base::getwd(), pos = 'CheckExEnv') > cleanEx() > nameEx("CGHdata-class") > ### * CGHdata-class > > flush(stderr()); flush(stdout()) > > ### Name: CGHdata-class > ### Title: Class "CGHdata" > ### Aliases: CGHdata-class [,CGHdata-method initialize,CGHdata-method > ### getuniKmax,CGHdata-method getmultiKmax,CGHdata-method > ### multiseg,CGHdata-method uniseg,CGHdata-method print,CGHdata-method > ### show,CGHdata-method summary,CGHdata-method > ### Keywords: classes > > ### ** Examples > > #set.seed(1) > #simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > #CGHd = new("CGHdata",Y=simul$Y) > #show(CGHd) > #summary(CGHd) > > > > cleanEx() > nameEx("CGHoptions-class") > ### * CGHoptions-class > > flush(stderr()); flush(stdout()) > > ### Name: CGHoptions-class > ### Title: Class "CGHoptions" > ### Aliases: CGHoptions-class [,CGHoptions-method > ### calling<-,CGHoptions-method alpha<-,CGHoptions-method > ### beta<-,CGHoptions-method nblevels<-,CGHoptions-method > ### wavenorm<-,CGHoptions-method GCnorm<-,CGHoptions-method > ### itermax<-,CGHoptions-method print,CGHoptions-method > ### select<-,CGHoptions-method nbprocs<-,CGHoptions-method > ### show,CGHoptions-method > ### Keywords: classes > > ### ** Examples > > #CGHo = new("CGHoptions") > #showClass("CGHoptions") > > > > cleanEx() > nameEx("CGHresults-class") > ### * CGHresults-class > > flush(stderr()); flush(stdout()) > > ### Name: CGHresults-class > ### Title: Class "CGHresults" > ### Aliases: CGHresults-class [,CGHresults-method getbp,CGHresults-method > ### initialize,CGHresults-method getprobebias,CGHresults-method > ### getsegprofiles,CGHresults-method getlevels,CGHresults-method > ### getbgoutliers,CGHresults-method summary,CGHresults-method > ### Keywords: classes > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > CGHr = multiseg(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening ==48707== Mismatched free() / delete / delete [] ==48707== at 0x4839EFC: operator delete(void*) (/builddir/build/BUILD/valgrind-3.15.0/coregrind/m_replacemalloc/vg_replace_malloc.c:586) ==48707== by 0x165A9A29: colibri_c(double*, int*, int*, double*, double*, int*, double*) (packages/tests-vg/cghseg/src/colibri.cc:132) ==48707== by 0x49DEE6: do_dotCode (svn/R-devel/src/main/dotcode.c:1799) ==48707== by 0x4D193C: bcEval (svn/R-devel/src/main/eval.c:7054) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== Address 0x15eeca10 is 0 bytes inside a block of size 800 alloc'd ==48707== at 0x4839593: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.15.0/coregrind/m_replacemalloc/vg_replace_malloc.c:433) ==48707== by 0x165A956C: colibri_c(double*, int*, int*, double*, double*, int*, double*) (packages/tests-vg/cghseg/src/colibri.cc:61) ==48707== by 0x49DEE6: do_dotCode (svn/R-devel/src/main/dotcode.c:1799) ==48707== by 0x4D193C: bcEval (svn/R-devel/src/main/eval.c:7054) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== ==48707== Mismatched free() / delete / delete [] ==48707== at 0x4839EFC: operator delete(void*) (/builddir/build/BUILD/valgrind-3.15.0/coregrind/m_replacemalloc/vg_replace_malloc.c:586) ==48707== by 0x165A9A31: colibri_c(double*, int*, int*, double*, double*, int*, double*) (packages/tests-vg/cghseg/src/colibri.cc:133) ==48707== by 0x49DEE6: do_dotCode (svn/R-devel/src/main/dotcode.c:1799) ==48707== by 0x4D193C: bcEval (svn/R-devel/src/main/eval.c:7054) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== Address 0x15eec6b0 is 0 bytes inside a block of size 800 alloc'd ==48707== at 0x4839593: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.15.0/coregrind/m_replacemalloc/vg_replace_malloc.c:433) ==48707== by 0x165A94E8: colibri_c(double*, int*, int*, double*, double*, int*, double*) (packages/tests-vg/cghseg/src/colibri.cc:33) ==48707== by 0x49DEE6: do_dotCode (svn/R-devel/src/main/dotcode.c:1799) ==48707== by 0x4D193C: bcEval (svn/R-devel/src/main/eval.c:7054) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== ==48707== Mismatched free() / delete / delete [] ==48707== at 0x4839EFC: operator delete(void*) (/builddir/build/BUILD/valgrind-3.15.0/coregrind/m_replacemalloc/vg_replace_malloc.c:586) ==48707== by 0x49DEE6: do_dotCode (svn/R-devel/src/main/dotcode.c:1799) ==48707== by 0x4D193C: bcEval (svn/R-devel/src/main/eval.c:7054) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== Address 0x15eec350 is 0 bytes inside a block of size 800 alloc'd ==48707== at 0x4839593: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.15.0/coregrind/m_replacemalloc/vg_replace_malloc.c:433) ==48707== by 0x165A94DA: colibri_c(double*, int*, int*, double*, double*, int*, double*) (packages/tests-vg/cghseg/src/colibri.cc:32) ==48707== by 0x49DEE6: do_dotCode (svn/R-devel/src/main/dotcode.c:1799) ==48707== by 0x4D193C: bcEval (svn/R-devel/src/main/eval.c:7054) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== by 0x4E82B7: Rf_eval (svn/R-devel/src/main/eval.c:688) ==48707== by 0x4E9E76: R_execClosure (svn/R-devel/src/main/eval.c:1853) ==48707== by 0x4EAC53: Rf_applyClosure (svn/R-devel/src/main/eval.c:1779) ==48707== by 0x4DB76D: bcEval (svn/R-devel/src/main/eval.c:7022) ==48707== [multiseg] multisegmean running > bplist = getbp(CGHr) Ind1 Ind2 Ind3 Ind4 Ind5 > segprofiles = getsegprofiles(CGHr) > > > > cleanEx() > nameEx("alpha") > ### * alpha > > flush(stderr()); flush(stdout()) > > ### Name: alpha<- > ### Title: Replace slot alpha > ### Aliases: alpha<- alpha<--methods > ### Keywords: ts methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > alpha(CGHo) = 0.5 > > > > cleanEx() > nameEx("beta") > ### * beta > > flush(stderr()); flush(stdout()) > > ### Name: beta<- > ### Title: Replace slot beta > ### Aliases: beta<- beta<--methods > ### Keywords: ts methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > beta(CGHo) = 0.5 > > > > cleanEx() > nameEx("calling") > ### * calling > > flush(stderr()); flush(stdout()) > > ### Name: calling<- > ### Title: Replace slot calling > ### Aliases: calling<- calling<--methods > ### Keywords: ts methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > calling(CGHo) = TRUE > > > > cleanEx() > nameEx("getbgoutliers") > ### * getbgoutliers > > flush(stderr()); flush(stdout()) > > ### Name: getbgoutliers > ### Title: getbgoutliers > ### Aliases: getbgoutliers getbgoutliers-methods > ### Keywords: methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > calling(CGHo) = FALSE > wavenorm(CGHo) = "position" > CGHr.pos = multiseg(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening [multiseg] ILS running [multiseg] ILS running > wavenorm(CGHo) = "spline" > CGHr.spl = multiseg(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening [multiseg] ILS running [multiseg] ILS running > getbgoutliers(CGHr.pos,CGHr.spl,fdr=0.05) [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 [38] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 [75] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 > > > > cleanEx() > nameEx("getbp") > ### * getbp > > flush(stderr()); flush(stdout()) > > ### Name: getbp > ### Title: Breakpoints frequencies > ### Aliases: getbp getbp-methods > ### Keywords: ts methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > CGHr = multiseg(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening [multiseg] multisegmean running > bplist = getbp(CGHr) Ind1 Ind2 Ind3 Ind4 Ind5 > > > > cleanEx() > nameEx("getmultiKmax") > ### * getmultiKmax > > flush(stderr()); flush(stdout()) > > ### Name: getmultiKmax > ### Title: Determining the maximum number of segments for joint > ### segmentation > ### Aliases: getmultiKmax getmultiKmax-methods > ### Keywords: ts methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > > ## Default initialization > > multiKmax = getmultiKmax(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening > > ## Modification of all fields: > ## to set multiKmax to 50% of the sum of Kmax of each profile: > > beta(CGHo) = 0.5 > multiKmax = getmultiKmax(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening > > ## Modification of uniKmax > > uniKmax = getuniKmax(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening > uniKmax[[1]] = 80 > multiKmax = getmultiKmax(CGHd,CGHo,uniKmax) > > ## Modification of uniKmax and multiKmax > ## and check that the value is correct wrt options > > uniKmax[[1]] = 80 > multiKmax = getmultiKmax(CGHd,CGHo,uniKmax,multiKmax = 10) > > > > cleanEx() > nameEx("getsegprofiles") > ### * getsegprofiles > > flush(stderr()); flush(stdout()) > > ### Name: getsegprofiles > ### Title: getsegprofiles > ### Aliases: getsegprofiles getsegprofiles-methods > ### Keywords: methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > calling(CGHo) = FALSE > CGHr = multiseg(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening [multiseg] multisegmean running > segprofiles = getsegprofiles(CGHr) > > > > cleanEx() > nameEx("getuniKmax") > ### * getuniKmax > > flush(stderr()); flush(stdout()) > > ### Name: getuniKmax > ### Title: Maximum number of segments per profile > ### Aliases: getuniKmax getuniKmax-methods > ### Keywords: ts methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > > ## Default initialization > > uniKmax = getuniKmax(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening > > ## Modification of all fields: > ## to tune uniKmax using to 50% of the length of each profile : > ## this may provide the same result since a model selection is performed > ## for tuning. > alpha(CGHo) = 0.5 > uniKmax = getuniKmax(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening > > ## Modification of one field for one profile > ## and check that the value is correct wrt options > > uniKmax[[1]] = 80 > uniKmax = getuniKmax(CGHd,CGHo,uniKmax) > > > > cleanEx() > nameEx("itermax") > ### * itermax > > flush(stderr()); flush(stdout()) > > ### Name: itermax<- > ### Title: Replace slot itermax > ### Aliases: itermax<- itermax<--methods > ### Keywords: ts methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > itermax(CGHo) = 10 > > > > cleanEx() > nameEx("multiseg") > ### * multiseg > > flush(stderr()); flush(stdout()) > > ### Name: multiseg > ### Title: Joint segmentation of multiple CGH profiles > ### Aliases: multiseg multiseg-methods > ### Keywords: methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > CGHr = multiseg(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening [multiseg] multisegmean running > > > > cleanEx() > nameEx("nbprocs") > ### * nbprocs > > flush(stderr()); flush(stdout()) > > ### Name: nbprocs<- > ### Title: Replace slot nbprocs > ### Aliases: nbprocs<- nbprocs<--methods > ### Keywords: ts methods > > ### ** Examples > > CGHo = new("CGHoptions") > nbprocs(CGHo) = 2 # if 2 processes are available > > > > cleanEx() > nameEx("select") > ### * select > > flush(stderr()); flush(stdout()) > > ### Name: select<- > ### Title: Replace slot select > ### Aliases: select<- select<--methods > ### Keywords: ts methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > select(CGHo) = "mBIC" > CGHr = multiseg(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening [multiseg] multisegmean running > > select(CGHo) = "none" > CGHr = multiseg(CGHd,CGHo) [getuniKmax] if no selection is performed while uniKmax is not specified [getuniKmax] uniKmax initialized by pre-screening [check getmultiKmax] if no selection is performed while multiKmax is not specified [check getmultiKmax] multiKmax is initialized by default ********************************************************** [multiseg] multisegmean running > > > > cleanEx() > nameEx("simulprofiles") > ### * simulprofiles > > flush(stderr()); flush(stdout()) > > ### Name: simulprofiles > ### Title: simulprofiles > ### Aliases: simulprofiles > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > > > > cleanEx() > nameEx("uniseg") > ### * uniseg > > flush(stderr()); flush(stdout()) > > ### Name: uniseg > ### Title: Segmentation of univariate CGH profiles > ### Aliases: uniseg uniseg-methods > ### Keywords: methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > CGHr = uniseg(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening > > > > cleanEx() > nameEx("wavenorm") > ### * wavenorm > > flush(stderr()); flush(stdout()) > > ### Name: wavenorm<- > ### Title: Replace slot wavenorm > ### Aliases: wavenorm<- wavenorm<--methods > ### Keywords: ts methods > > ### ** Examples > > set.seed(1) > simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1) > CGHd = new("CGHdata",Y=simul$Y) [raw data: check] minimum number of positions per signal [raw data: check] number of records per position [raw data: format] changing dataframe to compact format (list) > CGHo = new("CGHoptions") > wavenorm(CGHo) = "spline" > CGHr = multiseg(CGHd,CGHo) [getuniKmax] uniKmax initialized by pre-screening [multiseg] ILS running [multiseg] ILS running > > > > ### *