r - 如何使用更少的包绘制二元正态分布的表面和轮廓

Question

我将绘制二元正态分布的 3D 曲面及其轮廓(可以是任何二元正态分布)。我想在绘图中使用 persp 和 contour。我在网上搜索了一下，但发现了很多方法。他们中的大多数人都使用过一些软件包。但我想以使用更少的软件包甚至不安装任何软件包的方式来执行此操作(但是，我不知道这是可能的)。您能否通过二元正态分布的示例向我展示您自己的方式？[更少的包，更好的解决方案]。事实上，我希望看到绘制二元正态分布的表面和轮廓的最简单方法。

Answer 1

只要您准备好编写自己的代码来执行此操作，您就可以在 R 中执行任何您喜欢的操作，而无需使用包。使用persp和contour只需几行R代码即 cocoa 视化二维表面，但是生成表面(即创建二维密度)情节，并非微不足道。

以下是如何仅使用 mnormt 包生成二维二元正态分布曲面(按照评论中的要求)。对于这个例子，我创建了一个虚拟方差-协方差矩阵。

library(mnormt)

x     <- seq(-5, 5, 0.25) 
y     <- seq(-5, 5, 0.25)
mu    <- c(0, 0)
sigma <- matrix(c(2, -1, -1, 2), nrow = 2)
f     <- function(x, y) dmnorm(cbind(x, y), mu, sigma)
z     <- outer(x, y, f)

contour(x, y, z)

persp(x, y, z, theta = -30, phi = 25, 
      shade = 0.75, col = "gold", expand = 0.5, r = 2, 
      ltheta = 25, ticktype = "detailed")

^{由reprex package于2020年6月17日创建(v0.3.0)}