R语言I绘制等高线图

等高线指的是地形图上高程相等的各点所连成的闭合曲线。而等高线图（contour map） 是可视化二维空间标量场的基本方法，可以将三维数据使用二维的方法可视化，同时用颜色视觉特征表示第三维数据

代码如下：

1.下载并调用MASS包

BiocManager::install("MASS")
library(MASS)

#matrix是创建矩阵
Sigma <- matrix(c(1,0.7,0.7,1),2,2)
Sigma

2.将画布调整为2*2的大小并生成随机数

par(mfrow = c(2, 2))

#基于MASS包的mvrnorm产生多元高斯分布的随机数
r <- mvrnorm(n=1000, c(0,3), Sigma)
head(r)

3.绘制等高线图
图一：

plot(density(r))

图二：

#二维核密度函数kde2d()来估计二维数据的密度函数
bivn.kde <- kde2d(r[,1], r[,2], n = 50)
#3D等高线图
persp(bivn.kde, phi = 40, theta = 35)