简介
小编方向所需的科研图形绘制,我已经分享了几期相关内容:
今天分享下,在模拟实验中如何实现以下图形:
注意:该图形主要展示不同(n,m)组合下,统计推断的估计性能(均方根误差,RMSE)比较。
该图不算很难,只需要把数据处理到合适格式,采用 ggplot2 中的 geom_point()
和 facet_wrap()
即可绘制得到。其中一些细节需要注意:
-
使用 latex2exp 包中的
TeX()
设置 Y 轴标签。 -
分面主题名称自定义,显示复杂的数学公式。
接下来,我们通过一个模拟数据来复现该图。本文数据和代码可见 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)