iterare su prodotto cartesiano di vettori

Question

Ho il seguente ciclo nidificato:

for (x in xs) { 
    for (y in ys) { 
     # Do something with x and y 
    } 
} 

Quale mi piacerebbe appiattire così ho pensato di costruire un prodotto cartesiano dei due vettori xs e ys e l'iterazione sul risultato . In Python, questo sarebbe banale:

for xy in product(xs, ys): 
    # x, y = xy[0], xy[1] 

Ma in R, l'equivalente più semplice che ho trovato sembra scoraggiante:

xys <- expand.grid(xs, ys) 
for (i in 1 : nrow(xys)) { 
    xy <- as.vector(xys[i, ]) 
    # x <- xy[1], y <- xy[2] 
} 

Sicuramente ci deve essere un modo migliore, no? (Per chiarire, io non voglio iterare su un indice di ... Penso che ci deve essere un modo per scorrere direttamente sopra le tuple nel prodotto.)

Answer 1

È possibile utilizzare la funzione apply per applicare una funzione a ciascuna riga del frame di dati. Basta sostituire "your function" con la tua funzione attuale.

# example data 
xs <- rnorm(10) 
ys <- rnorm(10)  

apply(expand.grid(xs, ys), 1, FUN = function(x) {"your function"}) 

Questo è un esempio molto semplice. Qui, la somma di entrambi i valori in una riga viene calcolato:

apply(expand.grid(xs, ys), 1, FUN = function(x) {x[1] + x[2]}) 

---

qui è una variante che utilizza argomenti denominati (xs, ys) invece di indici (x[1], x[2]):

myfun <- function(xs, ys) xs + ys 
arguments <- expand.grid(xs = rnorm(10), ys = rnorm(10)) 
apply(arguments, 1, function(x)do.call(myfun, as.list(x))) 

Answer 2

R ha un diverso paradigma di Python, così don' Mi aspetto che abbia generatori o tuple - abbiamo vettori e indici per questo.

In questo modo, per mappare una funzione su un prodotto cartesiano è sufficiente chiamare

outer(xs,ys,function(x,y) ...) 

e undim il risultato se lo si desidera.

EDIT: Nel caso xs o ys sono qualcosa di più complesso di vettori di base, una possibilità è quella di utilizzare indici, cioè

outer(seq(a=xs),seq(a=ys),function(xi,yi) ... xs[[xi]]/ys[xi,]/etc. ...) 

o mappare una funzione su un prodotto artigianale bit utilizzando mapply

mapply(function(x,y) ...,xs,rep(ys,each=length(xs)))