r - 在 R 中拟合 von Mises 分布的混合

Question

我有一组角度数据，我想将它们拟合到两个冯·米塞斯分布的混合中。如下所示，数据聚集在大约 0 和 ±π 处，因此这种情况需要具有周期性边界。

我尝试使用 movMF 包来拟合这些数据的分布，但它似乎正在对每一行进行标准化，并且由于这是一组一维数据，因此结果是一个 ±1 的向量。其他人如何在 R 中拟合这样的分布混合？

Answer 1

问题在于使用角度向量作为 movMF 函数的输入。相反，角度必须转换为单位圆上的点

pts_on_unit_circle <- cbind(cos(angle_in_degrees * pi / 180), 
                            sin(angle_in_degrees * pi / 180))
d <- movMF(pts_on_unit_circle, number_of_mixed_vM_fxns)
mu <- atan2(d$theta[,2], d$theta[,1])
kappa <- sqrt(rowSums(d$theta^2))

来源:联系 movMF 包的作者 Kurt Hornik。