Bonjour,
En tant que débutant sur R, je cherche à maximiser une vraisemblance, dépendante de 5 paramètres. Le code ci-dessous risque de vous faire comprendre :
ML = function(x1,x2,x3,x4,x5) (-sum((((sla-x1*(dbh^x2)-x3*exp(-x4*age))^2)/((sigi*dbh^x5)^2)+log(2*pi)+log((sigi*dbh^x5)^2))))
Donc ça, c'est la fonction que je veux minimiser (ou maximiser si j'enlève le signe - devant le sum).
Une fois avoir défini cette fonction, j'utilise la fonction nlm, dans le but d'obtenir mes 5 paramètres (x1, x2, ...) qui minimise ma fonction.
donc je tape le code suivant, suivi de la réponse de R :
> nlm(ML,p=c(12,0.2,6,0.5,-1))
Erreur dans f(x, ...) :
l'élément 2 est vide ;
la partie de la liste d'arguments de '^' en cours d'évaluation était :
(dbh, x2)
J'ai du mal à comprendre ce qui ne va pas dans mon code...
Quelqu'un peut-il m'expliquer comment fonctionne cette fonction nlm (je n'ai pas trouvé de réponse compréhensible lors de mes recherches), si elle peut calculer ce que je lui demande et si non, quelle fonction le peut ?
Merci !