Estimation des paramètres pour des lois skew-normales multivariées

Contrairement au cas normal multivarié, où l'on peut obtenir les estimations des paramètres à partir de la moyenne d'échantillon et de la matrice variance-covariance d'échantillon, les paramètres de la loi skew-normale doivent être estimés par un processus itératif. Vous pouvez utiliser les fonctions déjà offertes dans le paquet sn pour obtenir ces estimations.

Dans l'exercice précédent, nous avons examiné les variables Ht et Wt de l'objet ais.female et conclu qu'une loi skew-normale peut modéliser leur distribution. Nous allons maintenant estimer les paramètres de cette loi.