Discussion :

[boubouchon]

Bonjour à tous,

Sans boucle et sans package,

Je voudrais retourner N produits de M lignes choisis aléatoirement d'une matrice A :

Pour l'instant, j'ai réussi avec une boucle :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
B <- matrix(NA, ncol = ncol(A), nrow = 0) 
for (i in 1 : N) B <- rbind(B, apply(A[sample(1 : nrow(A), M, replace = T), ], 2, prod))

Avez-vous une solution ?

Merci d'avance !

[Sengar]

Re salut,

Honnetement tu ne pourras pas beaucoup optimiser ce code la.
Tout simplement parce que si tu enleves la boucle, ce sera pour ajouter une fonction de la gamme des apply( aggregate, tapply, sapply...) qui s'appliquera sur les marges 1:N. (donc revient un peu au même que la boucle...).

Exemple:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
l<-tapply(rep(M,N),1:N,function(x){A[sample(1 : nrow(A), M, replace = T), ]})
t(sapply(l,apply,2,prod))

Qui fait la même chose que ton code.
La boucle est implicite.

Il doit quand même y avoir un petit gain de temps car je ne fais pas de rbind à chaque étape. (le rbind obligeant à recharger le tableau à chaque incrémentation, contrairement à la liste qui se contente d'ajouter un element sans recharger tout ce qu'il y a avant)