ggplot2中的蠕虫图残差图答案

【问题标题】：Worm plot residuals graph in ggplot2ggplot2中的蠕虫图残差图
【发布时间】：2021-08-23 23:23:05
【问题描述】：

我正在尝试在使用 gamlss 包中的 gamlss 函数拟合的模型上绘制蠕虫图残差。兴趣图如下所示：

最初，下面是使用childsds包中的wormplot_gg函数的计算例程，但是，使用上述函数表示的结果与上面所示的示例不同，正在应用到 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)

以下是通过gamlss 包中的wp 函数获得的传统结果。

最后，我们通过childsds 包中的wormplot_gg 函数获得了结果。然而，正如已经描述的，这个并没有以我感兴趣的方式呈现自己，即第一个图形的视觉结构。

【问题讨论】：

你的意思是你想让情节看起来不同吗？喜欢主题和颜色？
我想获得一个具有第一张图像外观和感觉的图形，即具有ggplot2 包中的功能。
childsds::wormplot_gg 中的源代码看起来很简单（可能要复杂得多）。您可以删除 theme_minimal() 和 theme() 行以撤消自定义主题并删除 geom_() 行中的自定义
但是，该函数似乎不接受 ggplot2 包中的某些功能。我在ggplot2 中看到您有一个包含蠕虫图构建的站点，我什至尝试重现它以解决我的问题，但我做不到。考虑到通过gamlss 包进行调整的示例，您能否帮助我通过 ggplot2 构建此图？？

标签： r ggplot2 statistics regression gamlss

【解决方案1】：

使用 qqplotr https://aloy.github.io/qqplotr/index.html 和 detrend=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 而言，这只需使用函数 residuals 即可完成，因此您可以只执行例如 df <- data.frame(z=residuals(Model)) 然后只需继续剩下的代码

【讨论】：

喂！我想知道如何通过考虑 gamlss 模型的调整来完成...考虑到上述情况，您能否编辑您的答案？
对不起，我不明白你的意思。什么模板？
我想根据通过 gamlss 函数调整的模型结果来执行蠕虫图。对不起！我表达错误。
请注意，上面执行的蠕虫图是基于 gamls 模型的调整..