我有一个使用样条拟合的模型:
ssfit.3 <- smooth.spline(anage$lifespan ~ log(anage$Metabolic.by.mass),
df = 3)
我正在尝试获得模型诊断,如残差图和该模型的QQ图。我知道对于线性模型,你可以做
plot(lm)
其输出所有不同的绘图。由于plot(ssfit.3)
的输出不相同,我如何使用样条曲线模型进行此操作?
提取残差并使用qqnorm()
/qqline()
。
example(smooth.spline) ## to get a model to work with
qqnorm(residuals(s2m))
qqline(residuals(s2m))