多元偏态正态分布的参数估计
与多元正态不同,多元偏态正态分布的参数无法直接通过样本均值和样本方差-协方差矩阵得到,需要通过迭代过程估计。您可以利用 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$___