# getmultiKmax: Determining the maximum number of segments for joint... In cghseg: Segmentation Methods for Array CGH Analysis

## Description

Determining the total maximum number of segments according to the maximum number of segments for each profile

## Usage

 1 getmultiKmax(.Object,CGHo,uniKmax=NULL,multiKmax=NULL) 

## Arguments

 .Object an object of the CGHdata class CGHo an object of the CGHoptions class uniKmax NULL if no value is proposed, list of Kmax for each profile multiKmax NULL or a proposed value for multiKmax, to check validity

## Details

CGHo["beta"] sets the proportion of each Kmax in multiKmax. If uniKmax[[i]] is the Kmax of profile i then

\code{multiKmax} = β \times ∑_i \code{uniKmax[[i]]}

MultiKmax must be greater than the number of profiles, but should not be too small since it should be greater than the selected number of segments (unknown).

## Value

 multiKmax total number of segments for joint segmentation

## 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 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 #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") ## Default initialization #multiKmax = getmultiKmax(CGHd,CGHo) ## 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) ## Modification of uniKmax #uniKmax = getuniKmax(CGHd,CGHo) #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) 

cghseg documentation built on May 30, 2017, 6:56 a.m.

Search within the cghseg package
Search all R packages, documentation and source code