r - 使用 R 中的 sn 包从偏态正态分布中绘制数据

Question

我需要表征数据的偏态正态分布的位置、形状和尺度，然后使用这些参数从具有相同参数的偏态正态分布中随机抽取值。在过去，我使用snR 中的包来执行此操作。因此，例如，如果我有一个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)