Joint segmentation of multiple CGH profiles

Share:

Description

Joint segmentation of multiple CGH profiles

Usage

1
multiseg(.Object,CGHo,uniKmax=NULL,multiKmax=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.

multiKmax

total number of segments for joint segmentation

Value

CGHr

object of class CGHresult

Author(s)

F. Picard, E. Lebarbier, M. Hoebeke, G. Rigaill, B. Thiam, S. Robin

References

Joint segmentation, calling and normalization of multiple CGH profiles, Biostatistics (2011)

Examples

1
2
3
4
5
#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           = multiseg(CGHd,CGHo)