图中丢失的GAM残差

Question

我正在将GAM模型应用于我的数据:随着时间的推移，细胞丰度。该模型工作得很好(虽然我知道我的大小中有一种模式，但这是一个与此无关的不同问题)。它只是无法在最后的绘图中显示部分残差，尽管我设置了residuals = TRUE。这是我的输出：https://i.stack.imgur.com/C1MlY.png

我还用了mgcv包。

以前，这段代码按照我的要求工作，但在不同的数据上工作。任何关于它为什么不起作用的想法都是受欢迎的！

GAM_EA <- mgcv::gam(EUB_FISH ~ s(Day, by = Heatwave), data = HnH, method = "REML")
gam.check(GAM_EA) #Checking the model
mgcv::anova.gam(GAM_EA) #Retrieving the statistical results. See ?anova.gam
summary.gam(GAM_EA)
plot(GAM_EA, shift = coef(GAM_EA)[1], residuals = TRUE)

Answer 1

参见参数by.resid in ?plot.gam。它们在plot.gam中使用的方式对于因子而言是毫无意义的，除非您要子集部分残差，并且只绘制by因子的特定级别上的观测残差。