将第 50 个百分位数曲线 gamlss 合并在一个图中

Question

我是一位经验丰富的 R 程序员，试图帮助我的朋友使用 gamlss 可视化百分位数增长曲线。我已经成功地为每组绘制了曲线。现在的问题是她要求我将两组的第 50 个百分位合并到一个图中。我怎样才能实现这一目标？

这是我每组 50 个百分位数的代码：

男孩身高：

centiles(mboys, dataa$Age, cent=50, ylab="Berat Badan (kg)", xlab="Usia (tahun)", lwd.centiles = 2.5, col.centiles = "red", main = "Anak Laki-laki", lty.centiles = 1)

女生身高：

centiles(girlsmodel, girlsbhh$Age, cent = 50, ylab = "Berat Badan (kg)", xlab = "Usia (tahun)", lwd.centiles = 2.5, col.centiles = "black", main = "Anak Perempuan", lty.centiles = 1)

我想实现这个：

我尝试阅读文档，但仍然不知道该怎么做。我尝试过使用 chatGPT 但失败了。

Answer 1

尝试使用

centiles.pred

函数提取每个模型的百分位数。

centiles.pred(mboys,xname="Age",xvalues=seq(0,19,by=1))

xvalues

参数是进行预测的解释变量的新值（我假设男孩的年龄为 0 到 19 岁）。结果是

data.frame

，您可以将其与

girlsmodel

结合并绘图。