cghseg (version 1.0.5)

uniseg: Segmentation of univariate CGH profiles

Description

Segmentation of univariate CGH profiles

Usage

uniseg(.Object,CGHo,uniKmax=NULL)

Arguments

.Object

an object of the CGHdata class

CGHo

an object of the CGHoption class

uniKmax

list containing the maximum number of segments per profile. If !is.null(uniKmax) in the entry of the method, the method checks its validity.

Value

CGHr

object of class CGHresult

References

Joint segmentation, calling and normalization of multiple CGH profiles, Picard et al. Biostatistics (2011), Pruned dynamic programming for optimal multiple change-point detection, G. Rigaill, arXiv preprint arXiv:1004.0887 (2010)

Examples

# NOT RUN {
set.seed(1)
simul = simulprofiles(M=5,n=100,k.mean=2,SNR=5,lambda=1)
CGHd  = new("CGHdata",Y=simul$Y)
CGHo  = new("CGHoptions")
CGHr  = uniseg(CGHd,CGHo)
# }