我正在尝试在使用包中的gamlss函数拟合的模型上绘制蠕虫图残差gamlss。兴趣图如下所示：
最初，下面是使用包中函数的计算例程wormplot_gg，childsds但是，使用上述函数表示的结果与上面显示的示例不同，该示例应用于 R 中包含的数据集。
library(ggplot2)
library(gamlss)
library(childsds)
head(Orange)
Dados <- Orange
Model <- gamlss(circumference~age, family=NO,data=Dados); Model
wp(Model)
wormplot_gg(m = Model)
下面是通过包wp中的函数得到的传统结果gamlss。
最后，我们通过包中的wormplot_gg函数获得了结果childsds。然而，正如已经描述的那样，这个并没有以我感兴趣的方式呈现，即第一个图形的视觉结构。
使用带有选项的qqplotr https://aloy.github.io/qqplotr/index.htmldetrend=True
library(qqplotr)
set.seed(1)
df <- data.frame(z=rnorm(50))
ggplot(df, aes(sample=z)) +
stat_qq_point(detrend = T) +
stat_qq_band(detrend = T, color='black', fill=NA, size=0.5)
你也可以添加geom_hline(yintercept = 0)
编辑：
在将其与 gamlss 模型一起使用的情况下，首先必须从模型中提取随机残差，这对于 gamlss 只需使用 function 即可完成residuals，因此您可以只执行例如，df <- data.frame(z=residuals(Model)) 然后继续其余的的代码