Built matrix of total out of sample interval square error costs for held-out means. One indexed.

const_costs(y, w, min_seg, indices)

Arguments

y

NumericVector, values to group in order.

w

NumericVector, weights.

min_seg

positive integer, minimum segment size (>=1).

indices

IntegerVector, order list of indices to pair.

Value

xcosts NumericMatix, for j>=i xcosts(i,j) is the cost of partition element [i,...,j] (inclusive).

Examples


const_costs(c(1, 1, 2, 2), c(1, 1, 1, 1), 1, 1:4)
#>               [,1]          [,2]          [,3]          [,4]
#> [1,] 1.797693e+308  0.000000e+00  1.500000e+00  1.777778e+00
#> [2,]  0.000000e+00 1.797693e+308  2.000000e+00  1.500000e+00
#> [3,]  1.500000e+00  2.000000e+00 1.797693e+308  0.000000e+00
#> [4,]  1.777778e+00  1.500000e+00  0.000000e+00 1.797693e+308