在R中获取标准差椭圆几何形状

Question

我正在处理由点组成的数据框。将下面的数据作为我对个人 i 观察到的点的虚假示例：

D <- tibble(trip_long = c(-87.7, -87.6, -87.6, -87.7, -87.7, -87.8, -87.5, -87.6, -87.6, -87.6, -87.7, -87.6, -87.6, -87.6, -87.6, -87.6, -87.6, -87.7, -87.7, -87.6), trip_lat = c(41.7, 41.8, 41.7, 41.7, 41.7, 41.7, 41.6, 41.7, 41.8, 41.8, 41.7, 41.7, 41.7, 41.7, 41.7, 41.7, 41.8, 41.7, 41.7, 41.7) )

我想计算这个人的点的标准差椭圆，为此我使用了 aspace 库。

library(aspace)
sde <- calc_sde(id = 1, points = D))

sde 包含我的标准差椭圆：

plot_sde(sde)

我想提取该椭圆作为我所有观察结果的几何列，但我还没有弄清楚如何做到这一点。

我不完全确定这是否有可能。

Answer 1

尝试

suppressPackageStartupMessages({
  library(aspace) 
  library(dplyr)}
  )

D <- tibble(trip_long = c(-87.7, -87.6, -87.6, -87.7, -87.7, -87.8, -87.5, -87.6, -87.6, -87.6, -87.7, -87.6, -87.6, -87.6, -87.6, -87.6, -87.6, -87.7, -87.7, -87.6), trip_lat = c(41.7, 41.8, 41.7, 41.7, 41.7, 41.7, 41.6, 41.7, 41.8, 41.8, 41.7, 41.7, 41.7, 41.7, 41.7, 41.7, 41.8, 41.7, 41.7, 41.7) )

sde <- calc_sde(id = 1, points = D)
plot_sde(sde)

^{创建于 2023-10-18，使用 reprex v2.0.2}

并查看

calc_sde

和

plot_sde

的帮助页面以获取更多信息。输入

?calc_sde

和

?plot_sde

。