Bonjour

je souhaite optimiser les paramètres du mon modèle hydrologique (un vecteur x à 4 dimensions) afin que les débits qu'il simule soient proches des débits observés. En d'autres termes, chercher les valeurs de x permettant de maximiser les valeurs de NASH entre les débits observés et simulés. Or j'ai des valeurs manquantes dans mes données que j'ai exprimées par NA et ceci pose un problème avec la méthode d'optimisation L-BFGS-B que j'utilise.
voici le code:

Code : Sélectionner tout - Visualiser dans une fenêtre à part
res_optim = optim(x,fonction_ABCD,method="L-BFGS-B",lower=c(-10, -10, -10, -10),upper=c(500,1,1,20),control=list(fnscale=-1))
Le message d'erreur est le suivant :
Error in optim(x, fonction_ABCD, method = "L-BFGS-B", lower = c(-10, -10, :
L-BFGS-B needs finite values of 'fn'
Comment je peux procéder ?
merci