Nous travaillons actuellement sur un modèle et utilisons la fonction nlm() pour estimer les paramètres minimisant une fonction "vraisemblance", qui est en fait l'équivalent à une constante près de -log(V).
Les temps de calcul nous semblant plutôt longs et la fonction nlm() "explosant" régulièrement, je voulais savoir si quelqu'un avait utilisé une autre fonction de minimisation, peut-être plus performante.
Mon code est celui-ci :
Code : Tout sélectionner
param <- nlm(f = vraisemblance, p = param, sexe=sexe,
iterlim = 1000, stepmax = 1000, gradtol = 1e-6, steptol = 1e-6)
où vraisemblance() est une fonction de param et de sexe, qui fait intervenir dans son déroulement d'autres fonctions et des objets de l'environnement global.
Bonne journée,
Jacques