r - 如何在 R 中的 "y"轴的两侧绘制(几乎)相同的函数？

Question

我有一个函数，它取决于距离，并且根据您评估它的方向(东或西)表现不同。现在我有两个并排的图，但我需要将它们作为一个图，其中共享轴的标签并且 x 轴标签居中。这是代码现在的样子的示例:

x = (0:300)
par(mfrow=c(2,2))

Idriss70 = function(x){
exp(5.6315-0.4104*7-(2.9832-0.2339*7)*log(x+10)+0.00047*x+0.12)
}
plot(Idriss70(x), log = "x", type="l",xlim = c(300,1), xlab="Distancia [km]",ylab="PGA [g]", main="Aroma y Humayani extendida, Mw 7,0", col="green",   panel.first = grid(equilogs = TRUE))

Idriss70 = function(x){
ifelse (x >= 0 & x<=3, exp(5.6315-0.4104*7-(2.9832-0.2339*7)*log(0+10)+0.00047*0+0.12),
      exp(5.6315-0.4104*7-(2.9832-0.2339*7)*log(x+10)+0.00047*x+0.12))
}
plot(Idriss70(x), log = "x", type="l", xlab="Distancia [km]",ylab="PGA [g]", main="Aroma y Humayani extendida, Mw 7,0", col="green", panel.first = grid(equilogs = TRUE))

如您所见，绘图的对数比例不允许评估“负”值，因此我无法只使用一个绘图。

我如何在不使用 Illustrator 或其他图形软件的情况下将此图作为一个图，因为我必须为不同的区域创建很多这样的图？

我过去没有使用过 ggplot，但如果有必要，我愿意学习。

Answer 1

你基本上可以制作一个图并乱用 fig 这样你就可以将第一个图限制在设备的左半部分，第二个图限制在设备的右半部分

x <- 0:300

Idriss70 = function(x){
  exp(5.6315-0.4104*7-(2.9832-0.2339*7)*log(x+10)+0.00047*x+0.12)
}

Idriss71 = function(x){
  ifelse(x >= 0 & x<=3,
         exp(5.6315-0.4104*7-(2.9832-0.2339*7)*log(0+10)+0.00047*0+0.12),
         exp(5.6315-0.4104*7-(2.9832-0.2339*7)*log(x+10)+0.00047*x+0.12))
}

par(fig = c(0, .5, 0, 1), mar = c(5, 4, 4, 0), xaxs = 'i', ann = FALSE)
plot(Idriss70(x), log = "x", type="l", xlim = c(300,1),
     col="green", axes = FALSE, panel.first = grid(equilogs = TRUE))
xx <- axis(1)
axis(2)

par(fig = c(.5, 1, 0, 1), new = TRUE, mar = c(5, 0, 4, 2))
plot(Idriss71(x), log = "x", type="l", col="green",
     panel.first = grid(equilogs = TRUE), axes = FALSE)
axis(1, at = xx, labels = c('', xx[-1]))

title(xlab = "Distancia [km]", main = 'Aroma y Humayani extendida, Mw 7,0',
      ylab = 'PGA [g]', outer = TRUE, line = -1.5)

祝 ggplot 好运。你可能需要召唤@baptiste

我想你仍然可以使用 mfrow

graphics.off()
par(mfrow = c(1, 2), mar = c(5, 4, 4, 0), xaxs = 'i', ann = FALSE)
plot(Idriss70(x), log = "x", type="l", xlim = c(300,1),
     col="green", axes = FALSE, panel.first = grid(equilogs = TRUE))
xx <- axis(1)
axis(2)

par(mar = c(5, 0, 4, 2))
plot(Idriss71(x), log = "x", type="l", col="green",
     panel.first = grid(equilogs = TRUE), axes = FALSE)
axis(1, at = xx, labels = c('', xx[-1]))

title(xlab = "Distancia [km]", main = 'Aroma y Humayani extendida, Mw 7,0',
      ylab = 'PGA [g]', outer = TRUE, line = -1.5)