Description Usage Arguments Details Value Note Author(s) References See Also Examples
Low-level function for multivariate fused Lars segmentation (GFLars)
1 2 | segmentByGFLars(Y, K, weights = defaultWeights(nrow(Y)), epsilon = 1e-09,
verbose = FALSE)
|
Y |
A |
K |
The number of change points to find |
weights |
A |
epsilon |
Values smaller than epsilon are considered null. Defaults to |
verbose |
A |
This function recrusively looks for the best candidate
change point according to group-fused LARS. This is a low-level
function. It is generally advised to use the wrapper
doGFLars
which also works on data frames, has a
convenient argument stat
, and includes a basic workaround
for handling missing values.
See also jointSeg
for combining group fused
LARS segmentation with pruning by dynamic programming
(pruneByDP
).
See PSSeg
for segmenting genomic signals
from SNP arrays.
The default weights √{n/(i*(n-i))} are calibrated as suggested by Bleakley and Vert (2011). Using this calibration, the first breakpoint maximizes the likelihood ratio test (LRT) statistic.
A list with elements:
bkp |
A vector of |
lambda |
The estimated lambda values for each change-point |
mean |
A vector of length |
value |
A |
c |
\hat{c}, a |
This implementation is derived from the MATLAB code by Vert and Bleakley: http://cbio.ensmp.fr/GFLseg.
Morgane Pierre-Jean and Pierre Neuvial
Bleakley, K., & Vert, J. P. (2011). The group fused lasso for multiple change-point detection. arXiv preprint arXiv:1106.4199.
Vert, J. P., & Bleakley, K. (2010). Fast detection of multiple change-points shared by many signals using group LARS. Advances in Neural Information Processing Systems, 23, 2343-2351.
PSSeg
, jointSeg
, doGFLars
, pruneByDP
1 2 3 4 5 6 7 8 9 | p <- 2
trueK <- 10
sim <- randomProfile(1e4, trueK, 1, p)
Y <- sim$profile
K <- 2*trueK
res <- segmentByGFLars(Y, K)
print(res$bkp)
print(sim$bkp)
plotSeg(Y, res$bkp)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.