```
# NOT RUN {
resid_matrix <- matrix(rnorm(39 * 3),
nrow = 39,
ncol = 3,
dimnames = list(NULL, c('rumination', 'happiness', 'activity')))
resid_matrix[13, 2] <- 48
resid_matrix[23, 2] <- -62
resid_matrix[36, 2] <- 33
resid_matrix[27, 3] <- 75
resid_matrix
autovarCore:::residual_outliers(resid_matrix, 40)
# }
```

Run the code above in your browser using DataCamp Workspace