R语言做logistic回归分析限制性立方条图代码

以下是用R语言进行限制性立方条图的logistic回归分析代码示例：

# 导入数据
data <- read.csv("data.csv")

# 安装和加载必要的包
install.packages("ggplot2")
library(ggplot2)

# 限制性立方条图
ggplot(data, aes(x = factor(age), y = odds_ratio, fill = factor(sex))) +
  geom_bar(stat = "identity", position = "dodge", width = 0.5) +
  scale_fill_manual(values = c("#F8766D", "#00BFC4"), name = "性别",
                    labels = c("男", "女")) +
  scale_y_log10(name = "比值比") +
  labs(title = "年龄和性别对疾病风险的影响",
       x = "年龄组", y = "比值比") +
  theme_bw()

这段代码将数据集中的年龄和性别作为因子变量，并使用`ggplot2`包中的`geom_bar`函数生成限制性立方条图。其中，`odds_ratio`是数据集中的比值比变量。我们还使用`scale_fill_manual`函数设置性别的填充颜色，`scale_y_log10`函数将y轴刻度设为对数刻度，并使用`labs`函数添加图表标题和轴标签。最后，我们使用`theme_bw`函数将图表的背景设置为白色，以使图表更易读。