Función R para crear permutaciones monótonas (crecientes o decrecientes) [duplicar]
Estoy tratando de crear una función eficiente para generar todas las permutaciones que aumentan monótonamente de un vector grande. Obviamente, reducir las salidas de expand.grido gtools::permutationsfunciona, pero solo para vectores más pequeños.
Ejemplo:
x = 1:3
Salida deseada:
1, 1, 1
1, 1, 2
1, 1, 3
1, 2, 2
1, 2, 3
1, 3, 3
2, 2, 2
2, 2, 3
2, 3, 3
3, 3, 3
¿Alguna sugerencia sobre el uso de base R o paquetes existentes con esta capacidad?
EDITAR: Una solución ideal evitaría generar el conjunto completo de permutaciones en el subconjunto.
Respuestas
Usar data.tableesto es bastante fácil:
expand.monotonic <- function(x, len=length(x)){
do.call(CJ, lapply(integer(len), function(...) x ))[
eval(parse(text=paste0("V", 2:len, ">=", "V", 1:(len-1), collapse="&") )), ]
}
expand.monotonic(1:3)
V1 V2 V3
1: 1 1 1
2: 1 1 2
3: 1 1 3
4: 1 2 2
5: 1 2 3
6: 1 3 3
7: 2 2 2
8: 2 2 3
9: 2 3 3
10: 3 3 3
explicación:
Primero cree una lista que contenga los lentiempos de los vectores replicados . Use data.table::CJpara unir todos los vectores. Y aquí es donde ocurre la magia basada en la lencreación de una expresión básicamente V2>=V1&V3>=V2como V#es el nombre predeterminado para columnas sin nombre, y subconjunto por el resultado de evaluar dicha expresión.
parse(text=paste0("V", 2:len, ">=", "V", 1:(len-1), collapse="&") )
# expression(V2>=V1&V3>=V2)
Aquí hay un código que crea permutaciones con repeticiones permitidas como en su ejemplo, y detecta si cada permutación es monótona
x <- 1:3
# Generate permutations of length x
out <- gtools::permutations(length(x), length(x), v = x, repeats.allowed=TRUE)
# Detect if they're monotonic
mono <- apply(out, 1, function(x) { all(x == cummax(x)) })
output_with_monotonic_label <- cbind(out, mono)
# output_with_monotonic_label
# mono
# [1,] 1 1 1 1
# [2,] 1 1 2 1
# [3,] 1 1 3 1
# [4,] 1 2 1 0
# [5,] 1 2 2 1
# [6,] 1 2 3 1
# [7,] 1 3 1 0
# [8,] 1 3 2 0
# [9,] 1 3 3 1
# [10,] 2 1 1 0
# ....