beta <- matrix(
data = c(
0.7, 0.5, -0.1,
0.0, 0.6, 0.4,
0.0, 0.0, 0.5
),
nrow = 3
)
alpha <- rep(x = 1, times = 3)
lambda <- diag(3)
nu <- rep(x = 1, times = 3)
mean_eta <- SSMMeanEta(
beta = beta,
alpha = alpha
)
mean_y <- SSMMeanY(
nu = nu,
lambda = lambda,
mean_eta = mean_eta
)
SSMInterceptY(
mean_y = mean_y,
mean_eta = mean_eta,
lambda = lambda
)
Run the code above in your browser using DataLab