多變量偏態常態的參數估計
不像多變量常態可以用樣本平均與樣本變異—共變異數矩陣直接得到參數估計，偏態常態分配的參數需要透過反覆運算的程序來估計。你可以使用 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$___