This function extracts residuals from a fitted model estimated with `femlm`

, `feols`

or `feglm`

.

```
# S3 method for fixest
resid(
object,
type = c("response", "deviance", "pearson", "working"),
na.rm = TRUE,
...
)
# S3 method for fixest
residuals(
object,
type = c("response", "deviance", "pearson", "working"),
na.rm = TRUE,
...
)
```

## Arguments

- object
A `fixest`

object. Obtained using the functions `femlm`

, `feols`

or `feglm`

.

- type
A character scalar, either `"response"`

(default), `"deviance"`

, `"pearson"`

, or `"working"`

. Note that the `"working"`

corresponds to the residuals from the weighted least square and only applies to `feglm`

models.

- na.rm
Logical, default is `TRUE`

. Whether to remove the observations with NAs from the original data set. If `FALSE`

, then the vector returned is always of the same length as the original data set.

- ...
Not currently used.

## Value

It returns a numeric vector of the length the number of observations used for the estimation (if `na.rm = TRUE`

) or of the length of the original data set (if `na.rm = FALSE`

).

## Examples

```
# simple estimation on iris data, using "Species" fixed-effects
res_poisson = femlm(Sepal.Length ~ Sepal.Width + Petal.Length +
Petal.Width | Species, iris)
# we plot the residuals
plot(resid(res_poisson))
```