@Dirk Eddelbuettel: Основная идея отличная, но код, как показано, может быть улучшен. [Занимает много времени, чтобы объяснить, следовательно, отдельный ответ, а не комментарий.]
Функция hist()
по умолчанию рисует графики, поэтому вам нужно добавить параметр plot=FALSE
. Более того, яснее установить область графика с помощью вызова plot(0,0,type="n",...)
, в котором вы можете добавить метки оси, название сюжета и т. Д. Наконец, я хотел бы упомянуть, что можно также использовать затенение для различения двух гистограмм. Вот код:
set.seed(42)
p1 <- hist(rnorm(500,4),plot=FALSE)
p2 <- hist(rnorm(500,6),plot=FALSE)
plot(0,0,type="n",xlim=c(0,10),ylim=c(0,100),xlab="x",ylab="freq",main="Two histograms")
plot(p1,col="green",density=10,angle=135,add=TRUE)
plot(p2,col="blue",density=10,angle=45,add=TRUE)
И вот результат (слишком большой из-за RStudio :-)):
[/g0]