r - R のフォン ミーゼス分布の混合をあてはめる

Question

2 つのフォン ミーゼス分布の混合を当てはめたい角度データのセットがあります。以下に示すように、データは約 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))

出典: movMF パッケージの作成者である Kurt Hornik に問い合わせました。