segmentation.OU.specialCase.lasso performs the segmentation using a
lasso regression to select for the edges where the shifts are added.
segmentation.OU.specialCase.lasso(phylo, nbr_of_shifts, D, Xp,
penscale = rep(1, (nrow(phylo$edge) + 1)), ...)a phylogenetic tree
Number of shifts on the phylogeny allowed
moments of the conditional law of X given Y, result
of function compute_M.OU.specialCase
the selection strength
List containing : beta_0 : the optimal value at the root shifts : list containing the computed tau and delta costs : vector of costs
This function re-write the sum of costs to be minimized as a least squares
regression problem, and uses a lasso regression to solve it. It uses
functions incidence.matrix.full to express the problem as a
linear model.