微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何将大圆距离编码为 glmmtmb 混合效果模型?

如何解决如何将大圆距离编码为 glmmtmb 混合效果模型?

我正在尝试使用具有空间协方差结构的“glmmtmb”包运行混合效应模型,该结构考虑球体上点之间的距离。我已经深入研究了源代码并确定了我认为他们在哪里计算空间协方差结构的欧几里得距离。我知道基于此网站使用欧几里得距离: https://cran.r-project.org/web/packages/glmmTMB/vignettes/covstruct.html

通过调出源代码

trace(getReStruc,edit = T)

第 44 行是他们对距离矩阵使用 dist(coords) 的地方。

我想更改那个 = 代码,以便它计算大圆距离而不是欧几里得距离。但是,“geosphere”包中的 disthaversine() 等函数需要 4 个参数(x1 的纬度、x1 的经度、x2 的经度、x2 的经度),所以我不能直接插入:

geosphere::disthaversine(coords)

有没有人可以解决这个问题?任何帮助将不胜感激!

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。