使用R中的sn包從偏正態分布中繪制數據

Question

我需要表征數據的偏態正態分布的位置，形狀和比例，然后使用這些參數從具有相同參數的偏態正態分布中隨機抽取值。 過去，我是使用R中的sn包來完成此操作的。因此，例如，如果我有一個數據向量v ，它遵循偏正態分布，那么我將使用sn.em估計位置，形狀和比例功能：

require(sn)
data(ais)
v <- ais$Fe
sn.em(,v)

然后，我將使用rsnorm函數從具有相同參數的偏態正常dsitrubtion中隨機抽取：

rsnorm(100, shape = x, location = y, scale = z)

這兩個功能在sn軟件包中不再存在。 我該如何使用sn軟件包中的不同功能或一起使用不同的功能來做到這一點？

Answer 1

在“ sn”包中使用sn.mple（）函數。 例如，您可以使用

cp.est <- sn.mple(y=v,opt.method = "nlminb")$cp 
dp.est <- cp2dp(cp.est,family="SN")
dp.est
       xi     omega     alpha 
20.244158 73.840301  9.142412 

要從SN分布中抽取樣本，請在同一軟件包中使用rsn（）函數。 例如，

rsn(n=100, xi=20.24, omega=73.84, alpha=9.14)