复现 sci 顶刊中的分面散点图

简介

小编方向所需的科研图形绘制，我已经分享了几期相关内容：

1. R问题｜数值模拟流程记录和分享

2. R 语言数值实验中常见技巧整理

3. ggplot 绘制分面条形图并标记数字

4. 分面中添加不同表格

5. 分面中添加不同的直线

今天分享下，在模拟实验中如何实现以下图形：

注意：该图形主要展示不同（n,m）组合下，统计推断的估计性能（均方根误差，RMSE）比较。

该图不算很难，只需要把数据处理到合适格式，采用 ggplot2 中的 geom_point() 和 facet_wrap() 即可绘制得到。其中一些细节需要注意：

1. 使用 latex2exp 包中的 TeX() 设置 Y 轴标签。

2. 分面主题名称自定义，显示复杂的数学公式。

接下来，我们通过一个模拟数据来复现该图。本文数据和代码可见 GitHub。

教程

# 导入包
library(openxlsx)
library(ggplot2)
library(tidyverse)
library(ggsci)
library(latex2exp)

数据导入

原始数据结构如下所示：

sum_dat <- read.xlsx("数据汇总.xlsx", sheet = 1)

注意：可以导入多个 sheet，sum_dat <- sapply(1:5, function(i) read.xlsx("数据汇总.xlsx", sheet = i))

数据处理

首先，提取出 RMSE 的行。由于数据小数点后位数过多，我们进行尺度变化，乘以100。此时，数据格式如下所示

data.frame(sum_dat)