在R中绘制ellipse3d？

Question

包rgl 包含一个非常有用的函数ellipse3d，它可以返回一个覆盖3D 中95% 点的椭球体。然后可以在rgl::plot3d 中使用这个对象来绘制它。我的问题是，是否可以将ellipse3d 的输出转换为可以通过 js 绘图包（如 plotly）绘制的东西？

library(rgl)
dt <- cbind(x = rnorm(100), y = rnorm(100), z = rnorm(100))
ellipse <- ellipse3d(cov(dt))
plot3d(dt)
plot3d(ellipse, add = T, color = "red", alpha = 0.5)

那我该怎么做才能通过 plotly 绘制椭球？

Answer 1

您可以从ellipse$vb 中提取椭圆的坐标。然后绘制这些。比如：

p <- plot_ly() %>% 
  add_trace(type = 'scatter3d', size = 1, 
     x = ellipse$vb[1,], y = ellipse$vb[2,], z = ellipse$vb[3,], 
     opacity=0.01) %>% 
  add_trace(data=dt, type = 'scatter3d', x=~x, y=~y, z=~z)