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

Question

我有一组角度数据，我想要将两个von Mises分布混合在一起。如下所示，数据聚集在大约0和±π，因此在这种情况下需要具有周期性边界。

我已经尝试使用movMF包来适应这些数据的分布，但似乎它正在规范化每一行，因为这是一组1D数据，结果是±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))

资料来源：联系方式Kurt Hornik，movMF软件包的作者。