Bonjour, je cherche à implémenter l'algorithme de newton-raphson à la main pour les statistiques, pour le maximum de vraisemblance. On dispose d'une fonction unidimensionnelle avec un échantillon d'observation.

C'est en fait l'équivalent de la fonction nlm sous R, mais à programmer à la main.
Je ne connais pas la formule de mise à jour dans le cas unidimensionnel.