Get the samples from the posterior distribution of the residual covariance matrix, together with the posterior mean and quantiles.

## Arguments

- m
a model fitted with `jtdm_fit`

## Value

A list containing:

- Ssamples
Sample from the posterior distribution of the residual covariance matrix. It is an array where the first two dimensions are the rows and columns of the matrix, and the third dimensions are the samples from the posterior distribution

- Smean
Posterior mean of the residual covariance matrix.

- Sq975,Sq025
97.5% and 0.25% posterior quantiles of the residual covariance matrix.

## Examples

```
data(Y)
data(X)
# Short MCMC to obtain a fast example: results are unreliable !
m = jtdm_fit(Y=Y, X=X, formula=as.formula("~GDD+FDD+forest"), sample = 1000)
# get the inferred residual covariance
Sigma =get_sigma(m)
```