phi <- matrix(
data = c(
-0.357, 0.771, -0.450,
0.0, -0.511, 0.729,
0.0, 0.0, -0.693
),
nrow = 3
)
sigma <- matrix(
data = c(
0.24455556, 0.02201587, -0.05004762,
0.02201587, 0.07067800, 0.01539456,
-0.05004762, 0.01539456, 0.07553061
),
nrow = 3
)
lambda <- diag(3)
theta <- diag(3)
cov_eta <- LinSDECovEta(
phi = phi,
sigma = sigma
)
LinSDECovY(
lambda = lambda,
theta = theta,
cov_eta = cov_eta
)
Run the code above in your browser using DataLab