# NOT RUN {
## Generate data.
d <- GAMart()
## Estimate models.
b1 <- bamlss(num ~ s(x1), data = d)
b2 <- bamlss(num ~ s(x1) + s(x2) + s(x3), data = d)
## Extract quantile residuals.
e1 <- residuals(b1, type = "quantile")
e2 <- residuals(b2, type = "quantile")
## Plots.
plot(e1)
plot(e2)
# }
Run the code above in your browser using DataCamp Workspace