Hello all,

JAGS has some good matrix operations which other BUGS-esque programs lack, but

so far all programs (that I know of) cannot compute the Kronecker product.

Shouldn't it be theoretically possible, given JAGS's use of C++ matrix

libraries? (if someone has the time and sophistication to build a module).

My questions:

1) can the Kronecker product be computed "manually" through a series of

matrix-scalar multiplications? e.g.,

Omega <- Sigma1 * Sigma2

... or something along those lines. I guess it boils down to whether the

inheriting matrix Omega can be built-up from a series of sub-matrices.

2) is there a resource out there to refer to how to compile such a module from

C++ code?

Thanks,

Robert