將直方圖與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)