# Sets up and estimates a set of filter probabilities
h <- 2
m <- 2
TT <- 100
u <- array(rnorm(TT*m*h), c(TT, m, h))
Sigma <- array(diag(m), c(m, m, h))
Q <- matrix(c(0.99, 0.2, 0.01, 0.8), h, h)
x <- BHLK.filter(u, Sigma, Q)Run the code above in your browser using DataLab