Learn R Programming

networkTomography (version 0.3)

mle_filter: Filtering & smoothing at MLE for calibration SSM

Description

Run Kalman filtering and smoothing at calculated MLE for parameters of calibration SSM. This is used to obtain point and covariance estimates for the actual OD flows X following estimation of other parameters.

Usage

mle_filter(mle, Ft, yt, Zt, Rt, k = ncol(Ft), tau = 2, initScale = 1/(1 - diag(Ft)^2), nugget = sqrt(.Machine$double.eps))

Arguments

mle
numeric vector (length k+1) of parameters. theta[-1] = log(lambda), and theta[1] = log(phi)
Ft
evolution matrix (k x k) for OD flows; include fixed
yt
matrix (k x n) of observed link loads, one observation per column
Zt
observation matrix for system; should be routing matrix A
Rt
covariance matrix for observation equation; typically small and fixed
k
integer number of OD flows to infer
tau
numeric power parameter for mean-variance relationship
initScale
numeric inflation factor for time-zero state covariance; defaults to steady-state variance setting
nugget
small positive value to add to diagonal of state evolution covariance matrix to ensure numerical stability

Value

numeric marginal log-likelihood obtained via Kalman smoothinglist containing result of Kalman smoothing; see SSModel and KFS for details

References

A.W. Blocker and E.M. Airoldi. Deconvolution of mixing time series on a graph. Proceedings of the Twenty-Seventh Conference Annual Conference on Uncertainty in Artificial Intelligence (UAI-11) 51-60, 2011.

See Also

Other calibrationModel: calibration_ssm; llCalibration