`residuals.bamlss.Rd`

Function to compute quantile and response residuals.

- object
An object of class

`"bamlss"`

.- type
The type of residuals wanted, possible types are

`"quantile"`

residuals and`"response"`

residuals.- nsamps
If the fitted

`bamlss`

object contains samples of parameters, computing residuals may take quite some time. Therefore, to get a first feeling it can be useful to compute residuals only based on`nsamps`

samples, i.e.,`nsamps`

specifies the number of samples which are extracted on equidistant intervals.- x
Object returned from function

`residuals.bamlss()`

.- which
Should a histogram with kernel density estimates be plotted, a qq-plot or a worm plot?

- spar
Should graphical parameters be set by the plotting function?

- ...
For function

`residuals.bamlss()`

arguments passed to possible`$residuals()`

functions that may be part of a`bamlss.family`

. For function`plot.bamlss.residuals()`

arguments passed to function`hist.default`

and`qqnorm.default`

.

Response residuals are the raw residuals, i.e., the response data minus the fitted distributional
mean. If the `bamlss.family`

object contains a function `$mu(par, ...)`

, then
raw residuals are computed with `y - mu(par)`

where `par`

is the named list of fitted
values of distributional parameters. If `$mu(par, ...)`

is missing, then the fitted values
of the first distributional parameter are used.

Randomized quantile residuals are based on the cumulative distribution function of the
`bamlss.family`

object, i.e., the `$p(y, par, ...)`

function.

A vector of residuals.

Dunn P. K., and Smyth G. K. (1996). Randomized Quantile Residuals.
*Journal of Computational and Graphical Statistics* **5**, 236--244.

van Buuren S., and Fredriks M. (2001) Worm Plot: Simple Diagnostic Device for Modelling Growth
Reference Curves. *Statistics in Medicine*, **20**, 1259--1277

```
if (FALSE) ## Generate data.
d <- GAMart()
## Estimate models.
b1 <- bamlss(num ~ s(x1), data = d)
#> Error in bamlss.model.frame(formula, data, family, weights, subset, offset, na.action, specials, contrasts): object 'd' not found
b2 <- bamlss(num ~ s(x1) + s(x2) + s(x3), data = d)
#> Error in bamlss.model.frame(formula, data, family, weights, subset, offset, na.action, specials, contrasts): object 'd' not found
## Extract quantile residuals.
e1 <- residuals(b1, type = "quantile")
#> Error in residuals(b1, type = "quantile"): object 'b1' not found
e2 <- residuals(b2, type = "quantile")
#> Error in residuals(b2, type = "quantile"): object 'b2' not found
## Plots.
plot(e1)
#> Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'plot': object 'e1' not found
plot(e2)
#> Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'plot': object 'e2' not found
```