Stima dei parametri per skew-normal multivariate

Diversamente dal caso normale multivariato, in cui le stime dei parametri si ottengono con la media campionaria e la matrice di varianza-covarianza, i parametri della distribuzione skew-normal devono essere stimati con un procedimento iterativo. Puoi usare le funzioni già disponibili nel pacchetto sn per ottenere le stime.

Nell’esercizio precedente abbiamo esaminato le variabili Ht e Wt dell’oggetto ais.female e abbiamo verificato che una distribuzione skew-normal può modellarne la distribuzione. Ora stimiamo i parametri di questa distribuzione.