Переменные столбцы матрицы с заменой и кластеризацией

Question

Как перенести столбцы матрицы с заменой в R? Я нашел одну функцию, называемую rmperm {sna}, но она переставляет оба столбца и строки, тогда как я просто хочу переставить мои столбцы.

Редактировать: Я должен переставить матрицу 1000 раз, а затем выполнить иерархическую кластеризацию, чтобы у меня было окончательное дерево после 1000 рандомизаций.

Большое спасибо.

Answer 1

Попробуйте функцию sample().

> m <- matrix(as.integer(runif(9,0,9)),ncol=3) 
> m 
    [,1] [,2] [,3] 
[1,] 5 0 5 
[2,] 6 0 0 
[3,] 2 1 3 
> permuted <- m[,sample(ncol(m), 10, replace=TRUE)] 
> permuted 
    [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] 
[1,] 5 0 5 5 5 5 0 5 5  5 
[2,] 0 0 0 0 6 0 0 0 6  0 
[3,] 3 1 3 3 2 3 1 3 2  3 

Первый аргумент sample() определяет диапазон выборки (1:x), то второй аргумент определяет количество элементов на выбор (size), а параметр replace определяет, хотите ли мы использовать замену (требуется, если размер> x).