создание нескольких точечных диаграмм с одинаковыми осями в R

С помощью lattice и ggplot2 необходимо изменить форму данных. Например:

1. создайте 4 данных.рамка (x=x1, y=y1)...
2. добавьте столбец группы для каждого типа данных.рамка, группа=1,2,...
3. rbind данные 4.кадр в один раз

Вот пример использования lattice

dat <- data.frame(x = rep(sample(1:100,size=10),4),
                  y = rep(rnorm(40)),
                  group = rep(1:4,each =10))

xyplot(y~x|group,       ## conditional formula to get 4 panels
       data =dat,       ## data
       type='l',        ## line type for plot
       groups=group,     ## group ti get differents colors
       layout=c(2,2))   ## equivalent to par or layout

PS: нет необходимости устанавливать сакли. В xyplot настройки sacles по умолчанию - same (одинаковые sacles для всех панелей). Вы можете изменить его, например :

xyplot(y~x|group, data =dat, type='l',groups=group,
       layout=c(2,2), scales =list(y = list(relation='free')))

EDIT

Существует большое количество аргументов для построения решетки функции, позволяющие контролировать многие детали сюжета, вот например я настраиваю:

1. Текст, используемый для надписей и заголовков для полос
2. размер и размещение меток галочек оси,

3. Размер зазоров между столбцами и рядами панелей.

   xyplot(y~x|group, data =dat, type='l',groups=group,
         between =list(y=2,x=2),
         layout=c(2,2), 
         strip = myStrip,
         scales =list(y = list(relation='same',alternating= c(3,3))))
   

Где

myStrip <- function(var.name,which.panel, which.given,...) {
  var.name <- paste(var.name ,which.panel)
  strip.default(which.given,which.panel,var.name,...)
  }

EDIT In для того чтобы получить решетку сюжета base-graphics plots, вы можете попробовать следующее:

xyplot(y~x|group, data =dat, type='l',groups=group,
       between=list(y=2,x=2),
       layout=c(2,2), 
       strip =FALSE,
       xlab=c('a','a'),
       xlab.top=c('a','a'),
       ylab=c('b','b'),
       ylab.right = c('b','b'),
       main=c('plot1','plot2'),
       sub=c('plot3','plot4'),
       scales =list(y = list(alternating= c(3,3)),
                    x = list(alternating= c(3,3))))

Хотя ответ уже выбран, этот ответ использует ggplot, а не базу R, что и требовалось OP. Хотя ggplot действительно хорош для быстрого построения графиков, для публикации часто требуется более тонкий контроль над сюжетами, чем предлагает ggplot. Вот где базовый сюжет выделяется.

Я бы предложил прочитать объяснение Шона Андерсона магии, которую можно использовать с умным использованием par, а также несколько других хороших трюков, таких как использование layout() и split.screen(). Используя его объяснение, я пришел к следующему выводу:

# Assume that you are starting with some data, 
# rather than generating it on the fly
data_mat <- matrix(rnorm(600), nrow=4, ncol=150)
x_val <- iris$Petal.Width

Ylim <- c(-3, 3)
Xlim <- c(0, 2.5)

# You'll need to make the ylimits the same if you want to share axes


par(mfrow=c(2,2))
par(mar=c(0,0,0,0), oma=c(4,4,0.5,0.5))
par(mgp=c(1, 0.6, 0.5))
for (n in 1:4) { 
  plot(x_val, data_mat[n,], "p", asp=1, axes=FALSE, ylim=Ylim, xlim=Xlim)
  box()
  if(n %in% c(1,3)){
    axis(2, at=seq(Ylim[1]+0.5, Ylim[2]-0.5, by=0.5))
  }
  if(n %in% c(3,4)){
    axis(1, at=seq(min(x_val), max(x_val), by=0.1))
  }
}

Здесь еще есть над чем поработать. Так же, как и в ОП, данные выглядят раздавленными посередине. Конечно, было бы неплохо настроить все так, чтобы использовалась вся область построения.