ggplot2:如何在单个图形中组合直方图、地毯图和逻辑回归预测

# 加载必要的库
library(ggplot2)

# 创建一个示例数据集
set.seed(123)
data <- data.frame(
  x = rnorm(100),
  y = rbinom(100, 1, 0.5)
)

# 拟合逻辑回归模型
model <- glm(y ~ x, data = data, family = binomial(link = "logit"))

# 使用ggplot2创建图形
ggplot(data, aes(x = x)) +
  # 添加直方图
  geom_histogram(aes(y = ..density..), colour = 'black', fill = 'white') +
  # 添加地毯图
  geom_rug() +
  # 添加逻辑回归预测线
  stat_function(fun = function(x) predict(model, newdata = data.frame(x = x), type = "response"),
                geom = "line", color = "blue") +
  # 添加标题和标签
  labs(title = "Combining Histogram, Rug Plot, and Logistic Regression Prediction",
       x = "X-axis",
       y = "Density")