多変量歪正規分布のパラメータ推定
多変量正規分布では、標本平均と標本分散共分散行列を使ってパラメータを推定できます。しかし、歪正規分布のパラメータは反復的なプロセスによって推定する必要があります。sn パッケージの既存の関数を活用することで、推定値を求められます。
前の演習では、ais.female オブジェクトの Ht 変数と Wt 変数を調べ、歪正規分布がそのデータ分布をモデル化できると判断しました。ここでは、この分布のパラメータを推定しましょう。
この演習はコースの一部です
Rによる多変量確率分布
演習の手順
- 歪正規分布のパラメータを計算し、オブジェクト
fit.aisに代入してください。 - 歪度パラメータのみを出力してください。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Fit skew-normal parameters
fit.ais <- ___(y = cbind(ais.female$Ht, ais.female$Wt), opt.method = ___)
# Print the skewness parameters
fit.ais$dp$___