Как построить в цикле, используя par (new = T)?

Question

Я хотел бы построить два набора значений (из разных наборов данных, но для одних и тех же стран) с двумя осями x для нескольких стран в двух циклах, используя par (new = T). Однако R перекрывает второй набор данных, начиная с последней страны. Здесь вы можете увидеть результат https://www.dropbox.com/s/9x0dn7oom5czpwm/Rplot02.jpg?dl=0

Вот мой код в качестве примера только для 2 стран (у меня 28).

par(mfrow=c(2,2))
par(mar = c(5,5,5,5), cex.axis=1)
for (i in colnames(DPreachdifft[,2:3])) {
  plot(DPreachdifft$Diff, DPreachdifft[[i]],
       type="p",col="red", xaxt="n", yaxt="n", 
       ylim = range(c(-25, 25)),
       xlim = rev(range(DPreachdifft$Diff)),
       ylab="delay(+) or advance(-) in years", xlab=NA, pch=16,
       main = i, adj=0)
  axis(side = 1, at=ticks,labels=ticks) 
  axis(side = 2)
  abline(h=0)
}
par(new=T)
for (i in colnames(NPreachdifft[,2:3])) {
  plot(NPreachdifft$Diff1, NPreachdifft[[i]],
  type="p",col="blue", xaxt="n", yaxt="n",
  ylim = range(c(-40, 40)),
  xlab=NA, ylab=NA, pch=16)
axis(side = 3, at=ticks1)
mtext(side = 1, text="difference in DP", line = 2, adj = 0.5)
mtext(side = 3, text="difference in NP", line = 2, adj = 0.5)
}