Filter the seasonal DLM
Run Particle Gibbs Sampling on the first order DLM
Use Particle Gibbs to determine the parameters of the poisson DGLM
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