Skip to contents

Intercept from Steady-State Mean Vector for the Observed Variables in the State Space Model

Source: R/simStateSpace-ssm-intercept-y.R
SSMInterceptY.Rd

The intercept vector for the observed variables in the state space model \(\boldsymbol{\nu}\) is given by $$ \boldsymbol{\nu} = \mathrm{Mean} \left( \mathbf{y} \right) - \boldsymbol{\Lambda} \mathrm{Mean} \left( \boldsymbol{\eta} \right) $$ where \(\boldsymbol{\Lambda}\) is the matrix of factor loadings, \(\mathrm{Mean} \left( \mathbf{y} \right)\) is the steady-state mean vector for the observed variables, and \(\mathrm{Mean} \left( \boldsymbol{\eta} \right)\) is the steady-state mean vector for the latent variables.

Usage

SSMInterceptY(mean_y, mean_eta, lambda)

Arguments

mean_y

Numeric vector. Steady-state mean vector of the observed variables \(\mathrm{Mean} \left( \mathbf{y} \right)\).

mean_eta

Numeric vector. Steady-state mean vector of the latent variables \(\mathrm{Mean} \left( \boldsymbol{\eta} \right)\).

lambda

Numeric matrix. Factor loading matrix linking the latent variables to the observed variables (\(\boldsymbol{\Lambda}\)).

See also

Other Simulation of State Space Models Data Functions: LinSDE2SSM(), LinSDECovEta(), LinSDECovY(), LinSDEMeanEta(), LinSDEMeanY(), ProjectToHurwitz(), ProjectToStability(), SSMCovEta(), SSMCovY(), SSMInterceptEta(), SSMMeanEta(), SSMMeanY(), SimAlphaN(), SimBetaN(), SimBetaN2(), SimBetaNCovariate(), SimCovDiagN(), SimCovN(), SimIotaN(), SimNuN(), SimPhiN(), SimPhiN2(), SimPhiNCovariate(), SimSSMFixed(), SimSSMIVary(), SimSSMLinGrowth(), SimSSMLinGrowthIVary(), SimSSMLinSDEFixed(), SimSSMLinSDEIVary(), SimSSMOUFixed(), SimSSMOUIVary(), SimSSMVARFixed(), SimSSMVARIVary(), SpectralRadius(), TestPhi(), TestPhiHurwitz(), TestStability(), TestStationarity()

Author

Ivan Jacob Agaloos Pesigan

Examples

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
)
#> [1]  7.666667 14.333333 14.333333