Hola Ignacio
existe alguna función en r, que realice lo de a función proj_l1 de Python?
saludos.
Ah... no. Esa función la hicimos nosotros. Pero podés mirar el código y traducirla. Si querés te paso el pseudocódigo. Es muy sencillo.
dale
El mié, 23 nov 2022 a la(s) 19:55, Ignacio Ramirez (vía FING) (
El mié, 23 nov 2022 a la(s) 19:55, Ignacio Ramirez (vía FING) (
proj_l
Entrada: vector x de largo n
Salida: vector proyectado y
1. Si ||x||_1 <= 1, devolver x así como está
2. a <- vector de valores absolutos de vector x (a_i = |x_i|)
3. s <- vector de signos del vector x (s_i = -1 si z_i < 0, sino 1)
4. b <- a en orden decreciente
5. Repetir para j=1,...,n:
w <- ( sum_{i=1}^{i=j} - 1) / j
y <- max(b - w,0) % elemento a elemento
si ||y||_1 == 1: % aproximadamente
y <- s * y
reordenar elementos de y en su posicion original
devolver y como resultado
Entrada: vector x de largo n
Salida: vector proyectado y
1. Si ||x||_1 <= 1, devolver x así como está
2. a <- vector de valores absolutos de vector x (a_i = |x_i|)
3. s <- vector de signos del vector x (s_i = -1 si z_i < 0, sino 1)
4. b <- a en orden decreciente
5. Repetir para j=1,...,n:
w <- ( sum_{i=1}^{i=j} - 1) / j
y <- max(b - w,0) % elemento a elemento
si ||y||_1 == 1: % aproximadamente
y <- s * y
reordenar elementos de y en su posicion original
devolver y como resultado
Subí un util.R con proj_l1 en R.
No soy un experto en R así que probala con unos cuantos valores a ver si da cosas razonables.
Saludos.
No soy un experto en R así que probala con unos cuantos valores a ver si da cosas razonables.
Saludos.