#simulate two time series
x <- zoo_simulate(seed = 1)
y <- zoo_simulate(seed = 2)
#distance matrix
dist_matrix <- distance_matrix_cpp(
x = x,
y = y,
distance = "euclidean"
)
#least cost matrix
cost_matrix <- cost_matrix_orthogonal_cpp(
dist_matrix = dist_matrix
)
#least cost path
cost_path <- cost_path_orthogonal_cpp(
dist_matrix = dist_matrix,
cost_matrix = cost_matrix
)
nrow(cost_path)
#remove blocks from least-cost path
cost_path_trimmed <- cost_path_trim_cpp(
path = cost_path
)
nrow(cost_path_trimmed)
#auto sum
auto_sum_cpp(
x = x,
y = y,
path = cost_path_trimmed,
distance = "euclidean",
ignore_blocks = FALSE
)
Run the code above in your browser using DataLab