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]