R语言画散点图-饼图-折线图-柱状图-箱线图-等高线图-曲线图-热力图-雷达图-韦恩图(三D)

散点图

可以通过几种不同的方式绘制3D散点图。最常用的工具包括 plotly 包和 scatterplot3d 包

使用 plotly 包

plotly 是一个非常强大的交互式图形包,支持多种3D绘图功能,包括3D散点图

示例

# 加载plotly包
library(plotly)

# 创建示例数据
set.seed(123)
data <- data.frame(
  x = rnorm(100),
  y = rnorm(100),
  z = rnorm(100)
)

# 绘制3D散点图
fig <- plot_ly(data, x = ~x, y = ~y, z = ~z, type = 'scatter3d', mode = 'markers',
                marker = list(size = 5, color = 'rgb(23, 190, 207)', opacity = 0.8))
fig <- fig %>% layout(title = "3D散点图示例",
                      scene = list(xaxis = list(title = 'X轴'),
                                   yaxis = list(title = 'Y轴'),
                                   zaxis = list(title = 'Z轴')))
fig

解析

  • plot_ly() 函数用于创建散点图对象,type = ‘scatter3d’ 指定图形类型为散点图
  • layout() 用于设置图表的标题和轴标签
  • 生成的图形可以在Web浏览器中交互查看

效果

在这里插入图片描述

使用 scatterplot3d 包

scatterplot3d 包提供了一种简单的方式来创建3D散点图,虽然功能不如 plotly 和 rgl 强大,但对于基本的3D绘图非常有效

示例

确保已安装 scatterplot3d 包:

install.packages("scatterplot3d")
# 加载scatterplot3d包
library(scatterplot3d)

# 创建示例数据
set.seed(123)
x <- rnorm(100)
y <- rnorm(100)
z <- rnorm(100)

# 绘制3D散点图
scatterplot3d(x, y, z, color = "red", pch = 16,
               xlab = "X轴", ylab = "Y轴", zlab = "Z轴",
               main = "3D散点图示例")

解析

  • scatterplot3d() 函数用于创建基本的3D散点图
  • color 设置点的颜色,pch 设置点的类型

效果

在这里插入图片描述

饼图

使用 plotly 包

plotly 是一个强大的交互式图形包,能够创建各种复杂的图形,包括类似于3D饼图的效果。

示例

# 加载plotly包
library(plotly)

# 创建数据
labels <- c("A", "B", "C", "D")
values <- c(10, 20, 30, 40)

# 绘制3D饼图
fig <- plot_ly(labels = labels, values = values, type = 'pie', 
                hole = 0.3,  # 控制饼图中心空洞的大小
                textinfo = 'label+percent',
                marker = list(colors = c('#FF9999','#66B2FF','#99FF99','#FFCC99')))

fig <- fig %>% layout(title = '3D饼图效果示例',
                      showlegend = TRUE,
                      scene = list(
                        xaxis = list(title = 'X轴'),
                        yaxis = list(title = 'Y轴'),
                        zaxis = list(title = 'Z轴')
                      ))

fig

解析

  • plot_ly() 用于创建饼图对象,type = ‘pie’ 指定图形类型为饼图
  • hole 参数用于创建一个环形图,这样的效果更接近3D饼图的外观

效果

在这里插入图片描述

使用 plotrix 包

示例

# 加载plotrix包
library(plotrix)

# 创建数据
values <- c(10, 20, 30, 40)
labels <- c("A", "B", "C", "D")
colors <- c("red", "green", "blue", "yellow")

# 绘制3D饼图
pie3D(values, labels = labels, col = colors, explode = 0.1,
      main = "3D饼图示例", radius = 0.9)

解析

  • pie3D() 函数用于创建类似3D的饼图效果,虽然它不是完全的3D图形,但通过阴影和爆炸效果模拟了3D外观

效果

在这里插入图片描述

折线图

绘制 3D 折线图可以使用 plotly 或 rgl 等包。plotly 提供了直观的 3D 绘图功能,适合用于交互式图形,而 rgl 则适合用于需要更详细定制的 3D 图形

使用 plotly 包

示例

# 加载 plotly 包
library(plotly)

# 创建数据
x <- seq(1, 10, length.out = 100)
y <- sin(x)
z <- cos(x)

# 绘制 3D 折线图
fig <- plot_ly(x = x, y = y, z = z, type = 'scatter3d', mode = 'lines+markers')

# 添加标题和标签
fig <- fig %>% layout(
  title = '3D 折线图示例',
  scene = list(
    xaxis = list(title = 'X 轴'),
    yaxis = list(title = 'Y 轴'),
    zaxis = list(title = 'Z 轴')
  )
)

fig

解析

1. 创建数据:

  • x:自变量序列
  • y 和 z:因变量序列,分别对应 sin(x) 和 cos(x)

2. 绘制 3D 折线图:

  • plot_ly() 函数创建3D 折线图,type = ‘scatter3d’ 指定 3D 散点图(包含折线),mode = ‘lines+markers’ 指定显示折线和标记点

3. 设置图形属性:

  • layout() 函数设置图形的标题和坐标轴标签。

效果

在这里插入图片描述

使用 rgl 包

示例

# 加载 rgl 包
library(rgl)

# 创建数据
x <- seq(1, 10, length.out = 100)
y <- sin(x)
z <- cos(x)

# 绘制 3D 折线图
plot3d(x, y, z, type = 'l', col = 'blue', lwd = 2, xlab = 'X 轴', ylab = 'Y 轴', zlab = 'Z 轴')

# 添加标题
title3d(main = '3D 折线图示例')

# 设置视角
rgl.viewpoint(theta = 45, phi = 30)

解析

1. 创建数据:

  • x、y 和 z 为 3D 数据的坐标

2. 绘制 3D 折线图:

  • plot3d() 函数绘制 3D 折线图,type = ‘l’ 指定绘制线条,col 设置颜色,lwd 设置线宽

3. 设置图形属性:

  • title3d() 函数设置图形标题
  • rgl.viewpoint() 设置视角以获得更好的 3D 视图效果

效果

在这里插入图片描述

柱状图

示例

# 加载 rgl 包
library(rgl)

# 创建数据
x <- c(1, 2, 3, 4)
y <- c(1, 2, 3, 4)
heights <- c(10, 20, 15, 25)

# 绘制 3D 柱状图
for (i in seq_along(x)) {
  x_pos <- x[i]
  y_pos <- y[i]
  z_height <- heights[i]
  
  # 生成柱状体的顶点
  x_corners <- c(x_pos - 0.5, x_pos + 0.5, x_pos + 0.5, x_pos - 0.5)
  y_corners <- c(y_pos - 0.5, y_pos - 0.5, y_pos + 0.5, y_pos + 0.5)
  z_corners <- c(0, 0, 0, 0)
  
  # 绘制柱状体
  rgl::polygon3d(x_corners, y_corners, z_corners, col = 'blue', alpha = 0.6)
  rgl::polygon3d(x_corners, y_corners, z_corners + z_height, col = 'red', alpha = 0.6)
  
  # 绘制柱状体的侧面
  rgl::polygon3d(c(x_corners[1], x_corners[2], x_corners[2], x_corners[1]), 
                 c(y_corners[1], y_corners[1], y_corners[2], y_corners[2]),
                 c(0, 0, z_height, z_height), 
                 col = 'grey', alpha = 0.5)
}

# 设置视角和轴
rgl::rgl.viewpoint(theta = 45, phi = 30)
rgl::axes3d()
rgl::title3d(main = '3D 柱状图示例')

# 使图形窗口保持打开状态
rgl::rglwidget()

解析

1. 创建数据:

  • x 和 y 是柱状图的坐标位置
  • heights 是柱状图的高度

2. 使用 rgl:

  • polygon3d() 函数绘制柱状体的底面、顶部和侧面
  • rgl.viewpoint() 设置视角以获得更好的 3D 视图效果
  • rglwidget() 保持图形窗口打开状态,以便进行交互查看

效果

在这里插入图片描述

箱线图

示例

library(rgl)

# 生成一些示例数据
set.seed(123)
data <- matrix(rnorm(300), ncol = 3)

# 定义一个函数来绘制一个箱线图的长方体
draw_box <- function(x, y, z, width, height, depth, col = "blue") {
  # 绘制箱子的底面
  quads3d(x + c(-width/2, width/2, width/2, -width/2),
          y + c(-height/2, -height/2, height/2, height/2),
          z + rep(-depth/2, 4),
          col = col, alpha = 0.5)
  # 绘制箱子的顶面
  quads3d(x + c(-width/2, width/2, width/2, -width/2),
          y + c(-height/2, -height/2, height/2, height/2),
          z + rep(depth/2, 4),
          col = col, alpha = 0.5)
  # 绘制箱子的侧面
  quads3d(x + c(-width/2, width/2, width/2, -width/2),
          y + c(-height/2, -height/2, -height/2, -height/2),
          z + c(-depth/2, -depth/2, depth/2, depth/2),
          col = col, alpha = 0.5)
  quads3d(x + c(-width/2, width/2, width/2, -width/2),
          y + c(height/2, height/2, height/2, height/2),
          z + c(-depth/2, -depth/2, depth/2, depth/2),
          col = col, alpha = 0.5)
  quads3d(x + c(-width/2, width/2, width/2, -width/2),
          y + c(-height/2, height/2, height/2, -height/2),
          z + c(-depth/2, -depth/2, -depth/2, -depth/2),
          col = col, alpha = 0.5)
  quads3d(x + c(-width/2, width/2, width/2, -width/2),
          y + c(-height/2, height/2, height/2, -height/2),
          z + c(depth/2, depth/2, depth/2, depth/2),
          col = col, alpha = 0.5)
}

# 创建一个 3D 窗口
open3d()

# 绘制几个 3D 箱线图长方体
for (i in 1:ncol(data)) {
  # 计算长方体的边界
  q1 <- quantile(data[, i], 0.25)
  q3 <- quantile(data[, i], 0.75)
  median <- median(data[, i])
  lower <- min(data[, i])
  upper <- max(data[, i])
  
  # 绘制长方体
  draw_box(x = i, y = median, z = 0, width = 0.5, height = q3 - q1, depth = 1,
           col = "blue")
}

# 添加坐标轴标签
axes3d()
title3d(main = "3D Boxplot", xlab = "Sample", ylab = "Value", zlab = "Depth")

解析

  • draw_box 函数创建了一个 3D 长方体,模拟箱线图的外观
  • 使用 quads3d 函数绘制长方体的各个面
  • 对于每个数据列,计算了四分位数和中位数,并在 3D 空间中绘制了长方体

效果

在这里插入图片描述

直方图

示例

# 安装和加载必要的包
install.packages("rgl")
library(rgl)

# 生成示例数据
set.seed(123)
data <- matrix(rnorm(1000), ncol = 10)  # 10组数据,每组100个数据点

# 计算每组数据的频数
hist_data <- apply(data, 2, function(x) hist(x, breaks = 10, plot = FALSE))

# 清空之前的图形
open3d()

# 设定颜色
colors <- rainbow(length(hist_data))

# 绘制3D直方图
for (i in 1:length(hist_data)) {
  for (j in 1:(length(hist_data[[i]]$breaks) - 1)) {
    x0 <- hist_data[[i]]$breaks[j]    # 柱体底边左端
    x1 <- hist_data[[i]]$breaks[j + 1] # 柱体底边右端
    y0 <- i - 0.5                       # 柱体底边在Y轴上的位置
    y1 <- i + 0.5                       # 柱体顶边在Y轴上的位置
    z0 <- 0                            # 柱体底边的Z坐标
    z1 <- hist_data[[i]]$counts[j]    # 柱体顶边的Z坐标

    # 绘制底面
    quads3d(
      x = c(x0, x0, x1, x1), 
      y = c(y0, y1, y1, y0), 
      z = rep(z0, 4), 
      col = colors[i], 
      alpha = 0.5
    )

    # 绘制顶面
    quads3d(
      x = c(x0, x0, x1, x1), 
      y = c(y0, y1, y1, y0), 
      z = rep(z1, 4), 
      col = colors[i], 
      alpha = 0.5
    )

    # 绘制侧面
    quads3d(
      x = c(x0, x0, x0, x0), 
      y = c(y0, y1, y1, y0), 
      z = c(z0, z0, z1, z1), 
      col = colors[i], 
      alpha = 0.5
    )
    
    quads3d(
      x = c(x1, x1, x1, x1), 
      y = c(y0, y1, y1, y0), 
      z = c(z0, z0, z1, z1), 
      col = colors[i], 
      alpha = 0.5
    )
  }
}

# 添加坐标轴和标题
axes3d(edges = "bbox", xlab = "Value", ylab = "Group", zlab = "Frequency", box = TRUE)
title3d("3D Histogram", xlab = "Value", ylab = "Group", zlab = "Frequency")

解析

1. 坐标轴: 使用 axes3d 函数绘制坐标轴和标签,edges = “bbox” 表示显示边界框的坐标轴
2. 柱体: 每个柱体的底面、顶面和侧面都用 quads3d 绘制
3. 颜色: 使用 rainbow 函数生成颜色,并将其应用于每组数据的柱体

效果

在这里插入图片描述

等高线图

使用 plotly 包(方式1)

示例

# 安装和加载必要的包
install.packages("plotly")
library(plotly)

# 生成示例数据
n <- 50
x <- seq(-10, 10, length.out = n)
y <- seq(-10, 10, length.out = n)
z <- outer(x, y, function(x, y) sin(sqrt(x^2 + y^2)))

# 绘制3D等高线图
plot_ly(z = ~z, x = ~x, y = ~y, type = "surface") %>%
  layout(title = "3D Terrain Plot",
         scene = list(
           xaxis = list(title = "X"),
           yaxis = list(title = "Y"),
           zaxis = list(title = "Z")
         ))

解析

  • plot_ly 函数用于绘制图形,type = “surface” 指定绘制3D等高线图
  • layout 函数用于设置图形的标题和坐标轴标签

效果

在这里插入图片描述

使用 plotly 包(方式2)

示例

library(plotly)

# 生成示例数据
set.seed(123)
x <- rnorm(100)
y <- rnorm(100)
z <- rnorm(100)

# 创建等高线图数据
x_bins <- seq(min(x), max(x), length.out = 10)
y_bins <- seq(min(y), max(y), length.out = 10)
x_cut <- cut(x, breaks = x_bins)
y_cut <- cut(y, breaks = y_bins)

# 计算每个箱子中的数据量
hist_data <- table(x_cut, y_cut)

# 转换为矩阵格式
z_data <- as.matrix(hist_data)

# 绘制3D直方图
plot_ly(
  z = z_data,
  type = 'surface',
  colorscale = list(c(0, 1), c('blue', 'red')),
  colorbar = list(title = 'Count')
) %>%
  layout(
    title = "3D Histogram",
    scene = list(
      xaxis = list(title = 'X'),
      yaxis = list(title = 'Y'),
      zaxis = list(title = 'Count')
    )
  )

解析

1. 生成数据:

  • 生成 x, y, z 随机数据

2. 计算直方图数据:

  • 使用 cut 将 x 和 y 数据分箱
  • 使用 table 计算每个箱子的频数

3. 绘制等高线图:

  • 使用 plot_ly 和 type = ‘surface’ 绘制3D表面图
  • 传递 z_data 矩阵表示每个箱子的频数
  • colorscale 和 colorbar 用于设置颜色和图例

效果

在这里插入图片描述

曲线图

示例

# 加载必要的包
library(plotly)

# 生成 3D 螺旋曲线数据
t <- seq(0, 10 * pi, length.out = 100)
x <- sin(t)
y <- cos(t)
z <- t

# 创建数据框
data <- data.frame(
  t = t,
  x = x,
  y = y,
  z = z
)

# 绘制 3D 曲线图
fig <- plot_ly(
  data = data,
  x = ~x,
  y = ~y,
  z = ~z,
  type = 'scatter3d',
  mode = 'lines',
  line = list(width = 4, color = 'blue')
) %>%
  layout(
    scene = list(
      xaxis = list(title = 'X Axis'),
      yaxis = list(title = 'Y Axis'),
      zaxis = list(title = 'Z Axis')
    )
  )

# 显示图形
fig

解析

1. 生成 3D 螺旋曲线数据:

  • t 是参数,用于生成曲线的点
  • x 和 y 是通过 sin(t) 和 cos(t) 生成的曲线的横纵坐标
  • z 是随 t 增长的直线分量,用于使曲线在 z 轴上上升

2. 创建数据框:

  • 将生成的数据存储在数据框 data 中,以便用于绘图

3. 绘制 3D 曲线图:

  • plot_ly() 函数用于创建图形对象,type = ‘scatter3d’ 表示绘制 3D 散点图,而 mode = ‘lines’ 指定将点连接成线形成曲线
  • line 参数用于设置线条的宽度和颜色

4. 设置布局:

  • layout() 函数用于设置坐标轴的标题以及其他图形布局选项

效果

在这里插入图片描述

热力图

示例

# 加载必要的包
library(plotly)

# 生成模拟数据
set.seed(123)
x <- seq(-10, 10, length.out = 30)
y <- seq(-10, 10, length.out = 30)
z <- seq(-10, 10, length.out = 30)
data <- expand.grid(x = x, y = y, z = z)

# 计算热力图值,这里使用一个简单的函数作为例子
data$values <- with(data, sin(sqrt(x^2 + y^2 + z^2)))

# 绘制 3D 热力图
fig <- plot_ly(
  data,
  x = ~x,
  y = ~y,
  z = ~z,
  color = ~values,
  colors = colorRamp(c('blue', 'green', 'yellow', 'red')),
  type = 'mesh3d',
  intensity = ~values
) %>%
  layout(
    scene = list(
      xaxis = list(title = 'X Axis'),
      yaxis = list(title = 'Y Axis'),
      zaxis = list(title = 'Z Axis')
    ),
    coloraxis = list(colorbar = list(title = 'Intensity'))
  )

# 显示图形
fig

解析

1. 生成模拟数据:

  • x、y 和 z 是三维网格的坐标轴,使用 seq() 函数生成从 -10 到 10 的数据
  • expand.grid() 函数生成所有可能的 (x, y, z) 组合

2. 计算热力图值:

  • data$values 是我们用来表示热力图强度的数据,这里使用一个简单的 sin(sqrt(x^2 + y^2 + z^2)) 函数计算值,可以根据需要替换为实际数据或其他函数

3. 绘制 3D 热力图:

  • plot_ly() 函数创建图形对象,type = ‘mesh3d’ 指定绘制 3D 网格图,intensity 参数设置为热力图的值
  • colorRamp() 函数定义颜色渐变,展示不同强度的区域

4. 设置布局:

  • layout() 函数用于设置坐标轴标题和颜色条

效果

在这里插入图片描述

雷达图

示例

library(plotly)
library(dplyr)

# 生成示例数据
data <- data.frame(
  axis = c("A", "B", "C", "D", "E"),
  value = c(3, 4, 2, 5, 4)
)

# 转换数据为适合绘图的格式
n <- nrow(data)
theta <- seq(0, 2 * pi, length.out = n + 1)
r <- c(data$value, data$value[1])
x <- r * cos(theta)
y <- r * sin(theta)

# 创建3D雷达图
plot_ly(
  x = ~x,
  y = ~y,
  z = ~rep(0, length(x)),
  type = 'scatter3d',
  mode = 'lines+markers',
  line = list(color = 'blue'),
  marker = list(size = 5, color = 'blue')
) %>%
  layout(
    scene = list(
      xaxis = list(title = 'X'),
      yaxis = list(title = 'Y'),
      zaxis = list(title = 'Z', visible = FALSE)
    ),
    title = '3D Radar Plot'
  )

效果

在这里插入图片描述

韦恩图

示例

# 加载必要的包
library(plotly)

# 定义 3D 圆的参数
circle_radius <- 1
n_points <- 100
theta <- seq(0, 2 * pi, length.out = n_points)

# 计算圆的 x 和 y 坐标
circle_x <- circle_radius * cos(theta)
circle_y <- circle_radius * sin(theta)

# 生成 3 个圆的坐标
circle1 <- data.frame(x = circle_x, y = circle_y, z = 0)
circle2 <- data.frame(x = circle_x - 0.5, y = circle_y - 0.5, z = 0)
circle3 <- data.frame(x = circle_x + 0.5, y = circle_y - 0.5, z = 0)

# 创建一个 3D 图形
fig <- plot_ly() %>%
  # 添加圆1
  add_trace(
    x = circle1$x,
    y = circle1$y,
    z = circle1$z,
    type = 'scatter3d',
    mode = 'lines',
    line = list(width = 2, color = 'blue'),
    name = 'Circle 1'
  ) %>%
  # 添加圆2
  add_trace(
    x = circle2$x,
    y = circle2$y,
    z = circle2$z,
    type = 'scatter3d',
    mode = 'lines',
    line = list(width = 2, color = 'green'),
    name = 'Circle 2'
  ) %>%
  # 添加圆3
  add_trace(
    x = circle3$x,
    y = circle3$y,
    z = circle3$z,
    type = 'scatter3d',
    mode = 'lines',
    line = list(width = 2, color = 'red'),
    name = 'Circle 3'
  ) %>%
  layout(
    scene = list(
      xaxis = list(title = 'X Axis'),
      yaxis = list(title = 'Y Axis'),
      zaxis = list(title = 'Z Axis')
    ),
    title = '3D Venn Diagram'
  )

# 显示图形
fig

解析

1. 定义圆的参数:

  • circle_radius 为圆的半径
  • n_points 设置圆上点的数量以确保圆的平滑

2. 计算圆的坐标:

  • 使用 cos 和 sin 函数计算圆的 x 和 y 坐标
  • 创建三个不同位置的圆(通过调整 x 和 y 的位置)

3. 创建 3D 图形:

  • 使用 plot_ly() 创建图形对象
  • add_trace() 函数用于添加每个圆的坐标
  • layout() 设置坐标轴标题和图形标题

效果

在这里插入图片描述

最近更新

  1. docker php8.1+nginx base 镜像 dockerfile 配置

    2024-07-22 06:22:03       52 阅读
  2. Could not load dynamic library ‘cudart64_100.dll‘

    2024-07-22 06:22:03       54 阅读
  3. 在Django里面运行非项目文件

    2024-07-22 06:22:03       45 阅读
  4. Python语言-面向对象

    2024-07-22 06:22:03       55 阅读

热门阅读

  1. redis+spring面试题

    2024-07-22 06:22:03       17 阅读
  2. 线段树分治+可撤销并查集 学习笔记

    2024-07-22 06:22:03       16 阅读
  3. 含有罗马字母的txt转换为csv文件读取-报错

    2024-07-22 06:22:03       12 阅读
  4. go标准库---net/http服务端

    2024-07-22 06:22:03       13 阅读
  5. python多进程库(multiprocessing)

    2024-07-22 06:22:03       17 阅读