Elements 1-4 of x specify the sigmas for each of the four levels in
the toy example. Elements 5-7 specify the correlation scales for
level 1, elements 8-10 the scales for level 2, and so on.
Internal function pdm.maker()
shows how the B matrix is
obtained from the various elements of input argument x
. Note
how, in this simple example, the B matrices are diagonal, but
generalizing to non-diagonal matrices should be straightforward (if
you can guarantee that they remain positive definite).