segselect {segclust}R Documentation

segselect

Description

model selection for segmentation

Usage

        out <- segselect(J,Kmax,S=0.5)

Arguments

J Residual Sum of Squares for segmentation models up to Kmax segments
Kmax maximum number of segments
S threshold for the adaptive method, set to 0.5 (default)

Value

Kselect Selected number of segments

References

Picard, F., Robin, S., Lavielle, M., Vaisse, C., & Daudin, J. -J. (2005). A statistical approach for array CGH data analysis. BMC Bioinformatics, 6(1), 1-14.

Examples

        x1      <- rnorm(20,0,1)
        x2      <- rnorm(30,2,1)
        x       <- c(x1,x2)
        Kmax    <- 20
        vh      <- TRUE
        out     <- segmean(x,Kmax)
        Kselect <- segselect(out$J.est,Kmax)
        output  <- segout(x,K=Kselect,th = out$t.est,draw=TRUE)

[Package segclust version 0.74 Index]