R: Как сгенерировать все возможные векторы формы c (a, b, c), где a, b, c - целые числа, a <b <c <= 100?

Question

Теперь я могу получить эту работу сделали с помощью

x <- expand.grid(1:98, 1:99, 1:100); 
x.subset <- subset(x, Var1 < Var2 & Var2 < Var3) 

Но есть более элегантный способ? Я думаю, что шаг, на котором создается x, можно отбросить. Представьте, если я хочу сделать это для < b < c < 10,000. Распределяющая сетка потерпит неудачу на меньших машинах (я знаю, что x.subset тоже будет огромным!), Но у меня есть больше шансов создать его, если я могу пропустить шаг, на котором создается x.).

Answer 1

Невозможно довести до 10 000 единиц памяти (9998 * 9999 * 10000 == 9.997e + 11 элементов). Даже 1000 не подходит. Я думаю, у вас есть довольно хорошая версия прямо здесь, которая быстро работает примерно до 100, но намного выше ее просто невозможно.

Однако, возможно, пересмотреть свою «необходимость» для этого подхода и посмотреть, есть ли другой способ его достижения.

Answer 2

combn должен быть полезен здесь. Пример:

combn(1:4, 3) 
#  [,1] [,2] [,3] [,4] 
#[1,] 1 1 1 2 
#[2,] 2 2 3 3 
#[3,] 3 4 4 4 

Однако до 10.000 все равно будет проблематично и есть, вероятно, некоторые лучше алгоритм, чтобы добраться до конечной цели.

Answer 3

Публикация как ответ для удобочитаемости. Мое время показывает выигрыш cgrid. R302, Windows7, процессор i7 (конечно, одноядерный).

Rgames> cgrid <-function(x,y,z) { 
+ x <- expand.grid(1:98, 1:99, 1:100); 
+ x.subset <- subset(x, Var1 < Var2 & Var2 < Var3) 
+ } 
Rgames> ccomb<-function(z) combn(1:z,3) 
Rgames> microbenchmark(cgrid(1:98,1:99,1:100),ccomb(100),times=10) 
Unit: milliseconds 
        expr  min   lq  median   uq 
cgrid(1:98, 1:99, 1:100) 217.330657 220.113192 253.5495825 325.167895 
       ccomb(100) 336.545629 346.839341 358.8957100 380.154500 
     max neval 
346.787171 10 
452.278904 10