phi <- matrix(
data = c(
-0.357, 0.771, -0.450,
0.0, -0.511, 0.729,
0, 0, -0.693
),
nrow = 3
)
colnames(phi) <- rownames(phi) <- c("x", "m", "y")
sigma <- matrix(
data = c(
0.24455556, 0.02201587, -0.05004762,
0.02201587, 0.07067800, 0.01539456,
-0.05004762, 0.01539456, 0.07553061
),
nrow = 3
)
# Specific time interval ----------------------------------------------------
TotalCentralStd(
phi = phi,
sigma = sigma,
delta_t = 1
)
# Range of time intervals ---------------------------------------------------
total_central_std <- TotalCentralStd(
phi = phi,
sigma = sigma,
delta_t = 1:30
)
plot(total_central_std)
# Methods -------------------------------------------------------------------
# TotalCentralStd has a number of methods including
# print, summary, and plot
total_central_std <- TotalCentralStd(
phi = phi,
sigma = sigma,
delta_t = 1:5
)
print(total_central_std)
summary(total_central_std)
plot(total_central_std)
Run the code above in your browser using DataLab