fitting de loi de pareto

Bonjour,

je travaille sur la modélisation de la sévérité en réassurance, j'utilise la fonction fitdist du package fitdistrplus qui marche plutôt bien sauf quand je veux estimer les paramètres d'une loi de pareto, j'ai aussi essayé avec la fonction mledist mais j'ai toujours des erreurs et pas de paramètres.

Je fais donc appel à votre aide, j'aimerais savoir si vous connaissez un moyen d'estimer les paramètres d'une loi de pareto et une loi de pareto généralisée suivant une distribution, je suis bloqué à cette étape de mon fitting parce que je ne teste que la lognormale, la weibull, la gamma et la loi exponentielle. Alors que, j'ai besoin de lois à queue épaisse.

Merci d'avance pour votre aide