两个变量正态分布

Question

有人会帮助我使用persp函数获取R中两个变量正态分布的3d图吗？我正在使用mvtnorm包，这有点令人困惑......

任何帮助都将受到高度赞赏。

Answer 1

根据您的评论，您似乎想要这样的内容：

show.2dnorm <- function(n, mean = rep(0, nrow(sigma)), sigma = diag(length(mean))) {
  require(mvtnorm)
  require(ggplot2)
  norm2d <- as.data.frame(rmvnorm(n, mean, sigma))
  colnames(norm2d) <- c('x', 'y')
  ggplot(norm2d, aes(x,y)) + geom_point()
}

# standard normal
show.2dnorm(1e4, c(0, 0))
# 0.6 correlation
show.2dnorm(1e4, sigma = matrix(c(1, 0.6, 0.6, 1), 2))