Малые кратные распределения выборки в одном графике

Question

В R у меня 9 образцов асимметрии и эксцесса в разных условиях. Я хотел бы построить распределение выборки двух статистик на соответствующих графиках, а затем нарисовать линию соответствия. Я думаю, что-то похожее на этот рисунок: ошибка дисперсии вокруг линии соответствия.

Хотя моделирование моих данных довольно длинное, следующий код даст аналогичные данные. Однако я хочу разместить свои данные на одном графике. Может быть, ящики внутри ячеек на оси X и общей оси Y.

sample.and.plot.function<-function(draws,mean,sd)
{
  sample<-rnorm(draws,mean,sd) #sample
  A<-hist(sample,plot=FALSE) #histogram

  plot(NULL, type = "n", xlim = c(0, max(A$counts)), ylim = c(range(A$breaks))) #false plot
  rect(0, A$breaks[1:(length(A$breaks) - 1)], A$counts, A$breaks[2:length(A$breaks)]) #plot sideways
}

par(mfrow=c(3,3)) #par
draws<-1000 #var
mean<-(-0.4) #var
sd<-0.1 #var

while(mean<=0.4)
{
  sample.and.plot.function(draws,mean,sd)#plot
  mean<-mean+0.1 #iterate
}

Я полагаю, что нет "правильного" способа решения этой проблемы, но у меня в голове образ, который я просто не знаю, как кодировать.

Я ценю ваше рассмотрение этой проблемы.