Description Usage Arguments Value Author(s) References See Also Examples
View source: R/mutualInfoVAR1.r
Evaluate, within the VAR(1) model, the mutual information between each variate at the current time point and those at a future time point.
1  mutualInfoVAR1(A, SigmaE, T)

A 

SigmaE 
Covariance 
T 
Positive 
Object of class numeric
with elements corresponding to the mutual informations. The jth element represents the mutual information of the jth variate at the current time point with all variates at the T
th time point from now.
Wessel N. van Wieringen <w.vanwieringen@vumc.nl>
Miok, V., Wilting, S.M., Van Wieringen, W.N. (2017), "Ridge estimation of the VAR(1) model and its time series chain graph from multivariate timecourse omics data", Biometrical Journal, 59(1), 172191.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18  # set dimensions
p < 3
n < 4
T < 10
# set model parameters
SigmaE < diag(p)/4
A < matrix(c(0.1, 0.3, 0.6, 0.5, 0.4, 0, 0.3, 0.5, 0.2),
byrow=TRUE, ncol=3)
# generate data
Y < dataVAR1(n, T, A, SigmaE)
# fit VAR(1) model
VAR1hat < ridgeVAR1(Y, 1, 1)
# impulse response analysis
mutualInfoVAR1(VAR1hat$A, solve(symm(VAR1hat$P)), T=5)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs in the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.