Bonjour,
Renaud Lancelot a écrit :Peu de différences de rapidité entre les boucles et les *apply (qui sont des interfaces pour des boucles for: regarder par exemple le code de tapply). Il va peut-être falloir passer par du C. Difficile à dire sans un petit exemple complet de ce que vous chercher à faire.
Renaud
Renaud, est possible d'exécuter du code R sous C/C++ ?
Vos calculs, il n'y a pas moyen de
1/ les mettre en LUT (Look Up Table), (une fois pour toutes),
2/ puis ensuite les lire dans la table au lieu de recalculer ?
(c'est un truc courant pour des intégrales/fonctions un peu lourdes quand l'espace d'entrée est discret et pas trop gros).
J'utilise cette approche, mais il y a des contraintes mémoires...
Déja 2 matrices 37268*37268 et R commence à pleurnicher.
D'ailleurs, après plus de 24h ce calcul, j'ai eu une erreur du genre: impossible d'allouer 8Go !, et pourtant j'étais sur un calculateur performant.
Et de toute façon, je suis obligé de comparer mes modèles ou histogramme entre eux suivant un certain critère (c'est là qu'interviennent "les" intégrales) afin d'en sélectionner le 'meilleur'. Et comme j'aimerais calculer le risque (espérance mathématique), je suis obligé de faire du monté carlo !!!