更像是“中等数据”，因为仅仅16k行就足以复制它。

工作状态：

library(ggplot2)
outputdir = "/tmp"
mtdf = NULL
for(i in 1:100){
  mti = mtcars[,c("disp","hp","wt","gear","carb","mpg","qsec","vs","am")]
  mti$disp = mti$disp + 10*rnorm(length(mti$disp))
  mti$hp = mti$hp + 10*rnorm(length(mti$hp))
  mti$wt = mti$wt + 10*rnorm(length(mti$wt))
  mtdf = rbind(mtdf,mti)
}
p = ggplot(mtdf,aes(disp,hp,color=wt)) + 
      geom_point()+
      facet_nested(gear ~ am,scales="free")
ggsave(plot=p,file=file.path(outputdir,"facet_scatter.pdf"),width=10,height=10)

不起作用（各种替代方法）：

下面是sessionInfo()的输出，显示了加载的包版本。这个MWE在Linux上运行。

library(ggplot2)
library(ggpubr)
library(ggh4x)

outputdir = "/tmp"

mtdf = NULL
for(i in 1:500){
  mti = mtcars[,c("disp","hp","wt","gear","carb","mpg","qsec","vs","am")]
  mti$disp = mti$disp + 10*rnorm(length(mti$disp))
  mti$hp = mti$hp + 10*rnorm(length(mti$hp))
  mti$wt = mti$wt + 10*rnorm(length(mti$wt))
  mtdf = rbind(mtdf,mti)
}
p = ggplot(mtdf,aes(disp,hp,color=wt)) + 
      geom_point()+
      facet_nested(gear ~ am,scales="free")
ggsave(plot=p,file=file.path(outputdir,"facet_scatter.png"),width=10,height=10)

p = ggplot(mtdf,aes(disp,hp,color=wt)) + 
      geom_point()+
      facet_grid(gear ~ am,scales="free")
ggsave(plot=p,file=file.path(outputdir,"facet_scatter.png"),width=10,height=10)

p = ggscatter(mtdf,x="disp",y="hp",color="wt")
p = facet(p,facet.by=c("gear","am"))
ggsave(plot=p,file=file.path(outputdir,"facet_scatter.png"),width=10,height=10)

三次尝试均会产生以下错误：

R version 4.1.1 (2021-08-10) -- "Kick Things"
Copyright (C) 2021 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

source("/mnt/md0/insshare/mkarikom/Active_Project_Backup/DURIAN/DURIAN/slurm/test.R", encoding = "UTF-8")
Error: `scale_id` must not be `NA`
Run `rlang::last_error()` to see where the error occurred.
> sessionInfo()
R version 4.1.1 (2021-08-10)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 20.04.3 LTS

Matrix products: default
BLAS/LAPACK: /usr/lib/x86_64-linux-gnu/libmkl_rt.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] ggh4x_0.2.1   ggpubr_0.4.0  ggplot2_3.3.5

loaded via a namespace (and not attached):
 [1] pillar_1.6.4     compiler_4.1.1   tools_4.1.1      jsonlite_1.7.2  
 [5] lifecycle_1.0.1  tibble_3.1.6     gtable_0.3.0     pkgconfig_2.0.3 
 [9] rlang_0.4.12     DBI_1.1.2        withr_2.4.3      dplyr_1.0.7     
[13] generics_0.1.1   vctrs_0.3.8      grid_4.1.1       tidyselect_1.1.1
[17] glue_1.6.0       R6_2.5.1         rstatix_0.7.0    fansi_1.0.0     
[21] carData_3.0-5    purrr_0.3.4      tidyr_1.1.4      farver_2.1.0    
[25] car_3.0-12       magrittr_2.0.1   scales_1.1.1     backports_1.4.1 
[29] ellipsis_0.3.2   assertthat_0.2.1 abind_1.4-5      colorspace_2.0-2
[33] ggsignif_0.6.3   utf8_1.2.2       munsell_0.5.0    broom_0.7.11    
[37] crayon_1.4.2

rang::last_error()

如错误消息所示：

> rlang::last_error()
<error/rlang_error>
`scale_id` must not be `NA`
Backtrace:
  1. ggplot2::ggsave(...)
  3. ggplot2:::grid.draw.ggplot(plot)
  5. ggplot2:::print.ggplot(x)
  7. ggplot2:::ggplot_build.ggplot(x)
  8. layout$train_position(data, scale_x(), scale_y())
  9. ggplot2:::f(..., self = self)
 10. self$facet$train_scales(...)
 11. ggplot2:::f(...)
 12. ggplot2:::scale_apply(layer_data, x_vars, "train", SCALE_X, x_scales)
Run `rlang::last_trace()` to see the full context.

rang::last_trace()

如slang::last_error()的输出所示

> rlang::last_trace()
<error/rlang_error>
`scale_id` must not be `NA`
Backtrace:
     █
  1. └─ggplot2::ggsave(...)
  2.   ├─grid::grid.draw(plot)
  3.   └─ggplot2:::grid.draw.ggplot(plot)
  4.     ├─base::print(x)
  5.     └─ggplot2:::print.ggplot(x)
  6.       ├─ggplot2::ggplot_build(x)
  7.       └─ggplot2:::ggplot_build.ggplot(x)
  8.         └─layout$train_position(data, scale_x(), scale_y())
  9.           └─ggplot2:::f(..., self = self)
 10.             └─self$facet$train_scales(...)
 11.               └─ggplot2:::f(...)
 12.                 └─ggplot2:::scale_apply(layer_data, x_vars, "train", SCALE_X, x_scales)

更新：

这可能是由于scales1.1.1包造成的