在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軟件包的作者。