объединить гистограмму с диаграммой рассеяния в R

Question

Я пытаюсь создать график с гистограммой и диаграммой рассеяния только на одном графике с использованием вторичной оси. Более подробно, вот пример данных:

#generate example data 

set.seed(1) 
a <- rnorm(200,mean=500,sd=35) 
data <- data.frame(a = a, 
        b = rnorm(200, mean=10, sd=2), 
        c = c(rep(1,100), rep(0,100))) 

# produce a histogram of data$a 
hist(a, prob=TRUE, col="grey") 

#add a density line 
lines(density(a), col="blue", lwd=2) 

#scatter plot 
plot(data$a,data$b,col=ifelse(data$c==1,"red","black")) 

То, что я хочу сделать, это объединить гистограмму и график рассеяния вместе. Это означает, что моя ось x будет data$a, моя основная ось y равна frequency/density для гистограммы, а моя вторичная ось y равна data$b.

Answer 1

Может быть что-то вроде этого ...

# produce a histogram of data$a 
hist(a, prob=TRUE, col="grey") 

#add a density line 
lines(density(a), col="blue", lwd=2) 

par(new = TRUE) 

#scatter plot 
plot(data$a,data$b,col=ifelse(data$c==1,"red","black"), 
    axes = FALSE, ylab = "", xlab = "") 
axis(side = 4, at = seq(4, 14, by = 2)) 

Answer 2

Там хороший блог об этом здесь http://www.r-bloggers.com/r-single-plot-with-two-different-y-axes/.

В основном, как блог описывает, что нужно сделать:

par(new = TRUE) 
plot(data$a,data$b,col=ifelse(data$c==1,"red","black"), axes = F, xlab = NA, ylab = NA) 
axis(side = 4)