将两个正态分布叠加在R中一个曲线图上的两个直方图上

install.packages("ggplot2")
library(ggplot2)

mu1 <- 0
sigma1 <- 1
mu2 <- 2
sigma2 <- 0.5

data1 <- rnorm(1000, mean = mu1, sd = sigma1)
data2 <- rnorm(1000, mean = mu2, sd = sigma2)

plot <- ggplot() +
  geom_histogram(data = data.frame(x = data1), aes(x = x, fill = "Distribution 1"), alpha = 0.5, binwidth = 0.2) +
  geom_histogram(data = data.frame(x = data2), aes(x = x, fill = "Distribution 2"), alpha = 0.5, binwidth = 0.2) +
  labs(title = "叠加的正态分布直方图", x = "值", y = "频数") +
  scale_fill_manual(values = c("Distribution 1" = "blue", "Distribution 2" = "red"))

print(plot)