R中的二元偏态正态分布？

Question

有谁能告诉我一个很好的方法来适应R中的双变量偏差正态分布？我已经看过像“sn”和“fMultivar”这样的包，但这些包都没有被证实有用。谢谢。

Answer 1

您可以在“ sn”包中使用msn.mple（）函数。假设您的数据是y，即n×2数据帧/矩阵。然后只需使用：

msn.mple(y)或msn.mple(y, penalty='Qpenalty')

当样本量较小时，建议使用刑罚=“ Qpenalty”选项，以避免形状参数出现爆炸估计。