Obtain 2D kernel density estimates given sufficient statistics for lambdas and the initial data Y0
get_kernel(lambdas, sigma2, dens_grid, N, T, Y0)
sufficient statistics for the intercept term
variance of the shocks
grid over which the density is to be computed
cross-sectional dimension of the data
time dimension of the data
initial observations of the dependent variable