Parameterestimering för multivariata skev-normalfördelningar
Till skillnad från den multivariata normalfördelningen, där parameterestimaten kan beräknas med hjälp av stickprovsmedelvärdet och stickprovets varians-kovariansmatris, måste parametrarna för skev-normalfördelningen estimeras genom en iterativ process. Du kan använda befintliga funktioner i paketet sn för att få estimaten.
I föregående övning undersökte vi variablerna Ht och Wt från objektet ais.female och fastställde att en skev-normalfördelning kan modellera deras fördelning. Nu ska vi estimera parametrarna för denna fördelning.
Den här övningen är en del av kursen
Multivariata sannolikhetsfördelningar i R
Övningsinstruktioner
- Beräkna parametrarna för skev-normalfördelningen och tilldela resultatet till objektet
fit.ais. - Skriv ut enbart skevhetsparametrarna.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# Fit skew-normal parameters
fit.ais <- ___(y = cbind(ais.female$Ht, ais.female$Wt), opt.method = ___)
# Print the skewness parameters
fit.ais$dp$___