R的ggplot2自定义绘图代码如何转换为Julia的等价实现？

你已经通过 RCall 在 Julia 环境中调用 ggplot2 库，只需对原有 R 语法做少量适配即可实现目标效果，完整可运行代码如下：

依赖安装（如未安装对应包可执行） #

using Pkg
Pkg.add(["CairoMakie", "DataFrames", "Effects", "GLM", "StatsModels", "StableRNGs", "RCall", "DataFramesMeta", "Chain"])

完整实现代码 #

using CairoMakie, DataFrames, DataFramesMeta, Chain, Effects, GLM, StatsModels, StableRNGs, RCall
@rlibrary ggplot2

# 原有数据生成逻辑
rng = StableRNG(42)
growthdata = DataFrame(; age=[13:20; 13:20],
                       sex=repeat(["male", "female"], inner=8),
                       weight=[range(100, 155; length=8); range(100, 125; length=8)] .+ randn(rng, 16))

mod_uncentered = lm(@formula(weight ~ 1 + sex * age), growthdata)

refgrid = copy(growthdata)
filter!(refgrid) do row
    return mod(row.age, 2) == (row.sex == "male")
end
effects!(refgrid, mod_uncentered)

refgrid[!, :lower] = @. refgrid.weight - 1.96 * refgrid.err
refgrid[!, :upper] = @. refgrid.weight + 1.96 * refgrid.err

df = refgrid

# 筛选男性最大体重数据，等价R代码的 t1 生成逻辑
t1 = @chain df begin
    @rsubset :sex == "male"
    @orderby -:weight
    first(1)
end

# 绘图代码，和你提供的R版ggplot参数完全对应
ggplot(df, aes(x=:age, y=:weight, group=:sex, shape=:sex, linetype=:sex, fill=:sex, colour=:sex)) +
  geom_line(position=position_dodge(width=0.15)) +
  geom_point(position=position_dodge(width=0.15)) +
  geom_errorbar(aes(ymin=:lower, ymax=:upper), width=0.1,
                linetype="solid", position=position_dodge(width=0.15)) +
  geom_ribbon(aes(ymin=:lower, ymax=:upper, fill=:sex, colour=:sex), alpha=0.2) +
  geom_text(data=t1, aes(x=:age, y=:weight, label=round(:weight, 1)), hjust=-0.25, size=7, show_guide=false) +
  scale_y_continuous(limits=[70, 150], breaks=collect(80:20:140)) +
  theme_classic() +
  scale_colour_manual(values=["orange", "blue"]) +
  guides(color=guide_legend(override.aes=Dict("linetype" => ["dotted", "dashed"])),
         linetype="none") +
  xlab("Age") + ylab("Average marginal effects") + ggtitle("Title") +
  theme(
    axis.title.y = element_text(color="Black", size=28, face="bold", hjust=0.9),
    axis.text.y = element_text(face="bold", color="black", size=16),
    plot.title = element_text(hjust=0.5, color="Black", size=28, face="bold"),
    legend.title = element_text(color="Black", size=13),
    legend.text = element_text(color="Black", size=16),
    legend.position="bottom",
    axis.text.x = element_text(face="bold", color="black", size=11),
    strip.text = element_text(face="bold", size=15)
  )

注意说明 #

• 代码完全保留了你提供的R版ggplot的所有样式配置，输出效果和目标示意图一致
• 如果不想依赖RCall调用ggplot2，也可以改用纯Julia的Gadfly或CairoMakie库实现相同效果，语法会有少量差异

内容的提问来源于stack exchange，提问作者Jack