Filter the seasonal DLM
Sample the state using FFBS algorithm
Use Gibbs sampling with Inverse Gamma priors on the observation variance and diagonal system covariance
Simulate data from a Seasonal DLM
Run backward smoothing on the seasonal DLM
Use Kalman Filtering to determine the parameters of the Student's t-distribution DGLM