Бутстрапирования пар в Matlab

Question

Я хотел бы выразить «одномерную самонастройку» на «многомерный bootsrapping», то есть на первом этапе я рисую случайно с заменой из одномерного вектора, используя этот код:

s = RandStream.getGlobalStream(); 
reset(s) 

n = 100000;     % # of independent random trials 
h = 52;      % horizon 
T = size(Resid_standard, 1); 

Resid_bootstrapped = Resid_standard(unidrnd(T, h, n)); 

Теперь основной вектор Resid_standard не является одномерным вектором, а матрицей Tx2, и я хочу не только нарисовать случайные числа, но и случайные пары.

• Как мне изменить код для этого?
• Выход в одномерном корпусе составляет матрицу 100000x50. Выход для двумерного случая будет трехмерным. Как я могу сохранить результаты?

Answer 1

Одно из решений состоят в сохранении вектора индекса, использовать линейную индексацию и конкатенацию результаты:

r_ind = unidrnd(T, h, n); 
Resid_bootstrapped = cat(3, Resid_standard(r_ind), Resid_standard(r_ind + T)); 

Resid_bootstrapped будет тогда матрица h × n × 2.

Это даже может быть сокращен в однострочника:

Resid_bootstrapped = reshape(Resid_standard(unidrnd(T, h, n), [1,2]), h, n, 2);