This is a method for estimating a single-changepoint
which takes advantage of the special structure
of the Gaussian graphical model. It cannot take
arbitrary black-box models like simulated_annealing
or brute_force
. However, it can still be run within
binary_segmentation
.
rank_one(data, theta_init, buff = 10L, regularizer = 1, tau = -1L,
max_iter = 25L, update_w = 1, update_change = 0.9, mapping_iter = 1L,
tol = 1e-05)
N x P Matrix corresponding to the raw data.
Initial value for theta estimate.
Distance to maintain from edge of sample.
Regularizing constant, lambda.
Initial Estimate for change-point.
Maximum number of rank-one updates to be run.
Step size for prox-gradient.
Proportion of update_w
to keep when
the algorithm fails to successfully estimate theta.
Number of mapping iterations.
Tolerance at which the algorithm stops running.
List containing the estimated change-point and theta values.