小编wei*_*ner的帖子

我有一个包含三列的数据框:ref,target,distance.每个ref都有一个到同一组目标的测量距离,我想得到每个ref的最小距离矢量.现在我用for循环做这个,但似乎应该有一种方法来矢量化它.

这是我的代码:

refs <- levels(data$ref)

result <- c()
for (ref in refs) {
    # Find the minimum distance for observations with the current ref
    # but be sure to protect against ref == target!
    best_dist <- min(data[data$ref == ref & data$target != ref,]$distance)
    result <- c(result, best_dist)
}

我注定要通过这种方式设置数据框,还是有一种很好的方法来对此进行矢量化？谢谢您的帮助!

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

我已经尝试使用movMF包来适应这些数据的分布,但似乎它正在规范每一行,因为这是一组1D数据,结果是±1的向量.其他人如何在R中拟合这样的分布混合？