Parameterinschatting voor multivariate skew-normalverdelingen
In tegenstelling tot de multivariate normaalverdeling, waar je de parameters kunt schatten met het steekproefgemiddelde en de steekproefvariantie-covariantiematrix, moeten de parameters van de skew-normalverdeling via een iteratief proces worden geschat. Je kunt bestaande functies uit het sn-pakket gebruiken om de schattingen te krijgen.
In de vorige oefening bekeken we de variabelen Ht en Wt uit het object ais.female en stelden we vast dat een skew-normalverdeling deze data goed kan modelleren. Nu gaan we de parameters van deze verdeling schatten.
Deze oefening maakt deel uit van de cursus
Multivariate kansverdelingen in R
Oefeninstructies
- Bereken de parameters van de skew-normalverdeling en wijs ze toe aan het object
fit.ais. - Print alleen de scheefheidsparameters.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Fit skew-normal parameters
fit.ais <- ___(y = cbind(ais.female$Ht, ais.female$Wt), opt.method = ___)
# Print the skewness parameters
fit.ais$dp$___